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CONFIDENCE INTERVALS FOR VARIANCE 
RATIOS SPECIFYING GENETIC HERITABILITY* 


Franxkuin A. GRAYBILL, FRANK Martin, AND GEorGE GopFREY** 


Oklahoma Agricultural and Mechanical College 


1. Summary. The purpose of this paper is to present a method for 
setting confidence intervals on the ratio of variances in the twofold 
analysis of variance model. Since the ratio of variance components is 
of some importance in the study of genetic heritability, a portion of the 
paper is couched in genetic terminology. 


2. Introduction. Phenotypic differences between individuals in 
most traits are partly due to differences in heredity and partly due to 
the differences in individuals’ environments. Each developed trait is 
the result of the action of genes, the action of the environment, and the 
interaction of the genes and the environment. Heritability is a quantita- 
tive description of the amount of hereditary variation in a trait. 

It is important for the livestock breeder to know which traits have 
some degree of heritability if he wants to make any permanent improve- 
ment in his livestock. The only permanent changes in livestock quality 
are genetic changes brought about by a breeding program that will 
bring together favorable gene combinations. 

It is thus important for the breeder to be able to estimate herita- 
bility of a certain trait. One method of estimating heritability is by 
the technique of the analysis of variance. 

If r sires are each mated at random to s dams and ¢ offspring result 
from each mating, the analysis of variance takes the form as outlined 
in Table 1, where o2 , 0% , and o are respectively the variance com- 
ponents between sires, dams, and individuals. 

Using the variance components from Table 1, there are three ways of 
estimating heritability, h?: 


2 2 2 
*This work was sponsored by the Oklahoma Agricultural Experiment Station Project Number 850, 


**Respectively Associate Professor of Mathematics, Research Assistant in Mathematics, and 
Professor of Poultry, 
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TABLE 1 
Analysis of Variance of Inherited Trait 
Source of Degrees of Mean Expected 
variation freedom square mean square 
Between sires n=r—-1 A; o3 = + to} + sto? 
Between dams within sires| mn, = r(s — 1) A; o3 = o2 + lo} 
Between full siblings m = rs(t — 1) A, =o 


If in Table 1 the expected and observed mean squares are equated, the 
resulting equations can be solved for o2 , 03 , and. o2. If these values are 
substituted in the formulas above, the point estimates of h’ are obtained. 
While point estimates give valuable information, it is also desirable to 
have a confidence interval estimate of h’. Reference [1] presents a 
well known method for setting confidence limits on h’ for a simpler case 
than the one presented here. Reference [2] gives a method for finding 
the approximate standard error of h’ for the case presented above. 

The purpose of this paper is to present a method for setting confidence 
limits on the ratio 


+ 
(1) 


for the method of mating explained above. 


3. The approximate distribution of a linear combination of chi-square 
functions. Consider the twofold classification model 


(2) Yin 


where = 1,2, --- ,r;7 = 1,2, ---, 8; andk = 1,2,---,¢# Itis 
assumed that the a; are distributed normally with mean zero and vari- 
ance a2. Similarly the b,; and the c,;, are distributed normally with 
zero means and variances o; and o% , respectively. It is further assumed 
that all the terms are uncorrelated. The analysis of variance takes the 
form given in Table 1. 

It is well known that 


are independently distributed as chi-square with n, , n, , and n; degrees 
of freedom, respectively. While it is known that the linear combination 
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of independent chi-square variates is distributed as a chi-square variate 
if and only if the coefficients of the linear combination are equal to 
unity, it seems reasonable to assume that the distribution of a linear 
combination of independent chi-square variates can be reasonably well 
approximated by a chi-square distribution even if the coefficients of 
the linear sum are not all equal to one. 

The method presented in this paper is somewhat patterned after the 
method of attack used by Satterthwaite for testing a hypothesis on 
variance components (reference [3]). That is to say, the method pro- 
posed in this paper consists of equating the moments of 


(3) where of = tata, 


to the moments of the function x7v)/N, where xin) represents a chi- 
square variate with N degrees of freedom. We will determine the con- 
stants, N, a, , a3 , and y, so that the first two moments of Y/N will be 
equal to the first two moments of x{v)/N. We have chosen o( so that o2 
and o} have equal coefficients. The importance of this will be seen later. 
If this linear combination of chi-square variates, Y/N, is closely approxi- 
mated by xiw)/N, then the ratio (Y/N)/(A,/o2) will be approximately 
distributed as Snedecor’s F, and this will be used to obtain approximate 
confidence limits on h’. : 
The moment generating function of Y/N is 


My,yx(6) = (1 — — 


where 


Expanding M y,y(@) we get 
= 1 + (mB, + 
+ [n.(nz + 2)Bz + + + 2)B5]6°/2! 
+ + 2)(n2 + + + 2)nsB2Bs + 3nyns(ns + 2)B2Bs 


+ 2)(ng + 4)B3]0°/3! [> B,(n, + 23) 


k-1 k-i-1 


+2 (*) + 2j) IT B,(n3 + 2m) 


t=1 


iS 

2 
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The moment generating function of xiv) /N is (1 — 20/N)~*”. Expand- 
ing into an infinite series we get 


N(N + 2)6°/2! | N(N + 2)(N + 4)6°/3! 
N? 


M = 1 6 + N3 


(5) 


k-1 ak 

(N+ 

If we equate the first moments of My,y(6) and M,.,y(0), we find 
that n.B, + n;B,; = 1. Substituting for the B,’s we find that 


(6) 203 _ 1. 


It is now possible to determine the values of a, , a; , and y. It follows 
from (6) that y = t, a. = (s — 1)/s, and a; = 1/s. 

N will be determined by equating the coefficients of 6°/2! in equations 
(4) and (5), i.e., by equating the second moments of Y/N and xiy)/N. 
We obtain 


1+ = + 2n.n;B.B; + N3(Ns + 


Substituting for the B,’s, we get 


(7) 
2 
N30202 + 


We see that the number of degrees of freedom, N, is a function of the 
unknown parameters , , and . 
If we let K = o2/(to} + 02), then 


Ltr tte +o.) LI + 
It is then easily verified that (7) simplifies to 


+ tK)’ 
+ + stK)° 


Substituting for n. , n; , a, , and a; , the formula for N becomes 


rs'(r — 1)(1 + tK)’ q 
~ — 1s — 1) + + stk)’ 


and 


N= 


(8) N 
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It has been suggested (reference [3]) that the parameter, K, be 
estimated by the analysis of variance, and the value of N be obtained 
by substituting the estimated value of K into (8). If this procedure 
is followed, the estimated value of K, (see Table I), which is 
(A; — A,)/tsA, , is easily computed. After a value for N has been 
determined, the method explained in section 5 can be used to set con- 
fidence limits on h’. 

While this method for determining N may be satisfactory for many 
purposes, we will present an alternative method which will be applicable 
in many important situations. For this alternative method, we will be 
able to set upper and lower bounds on N. 

Let 


? 2~ 


and let f and g be two positive constants such that f < h? < g. After 
substituting D and W into h’ and performing some algebraic manipu- 
lations, we arrive at the following inequality for K, the only unknown 
in determining N from (8): 


f g rr 
If we can show that N in (8) is a monotonic function of K, then we can 
substitute the two bounds of K into (8) and get bounds for N. 

To show that N in (8) is a monotonic function of K we proceed as 
follows: 

Since N is a rational function of K with non-vanishing denominator 
(when r > 1, s > 1, and éK > 0), N is a continuous function of K and 
has continuous derivatives of all orders. Taking 6N/5K we get 


_ 2rs’t(1+tK)[(r—1)(s— 1) 1)(14+ tK)*(1+8tK) 

By straightforward manipulations it can be shown (for r > 1, s > 1, 
t > 1,andtK > 0) that 6N/d5K < 0, and thus N is a monotonic function 


of K for the ranges of r, s, and ¢ specified above. Therefore the inequality 
N, < N < N, holds where N, and N, are the values of (8) when 


K 


g 
(2 — g(D + 1) + gtD 
and when 


fi 
@—fiD+) +f 


respectively. 
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If the model (2) represents a genetic model, and if for example in 
@ particular problem it seems reasonable to assume that .1 < h? < .6 
(say), and if the inheritance is additive and autosomal, then D = 1. 
Therefore the inequality 


is obtained. Using these values for K, bounds can be put on N that 
are applicable to the particular problem under investigation. On the 
other hand if the experimenter does not have enough information at 
his disposal to make any particular statement about the bounds on h’, 
it is known that in genetic models 0 < h? < 1. Also in many important 


genetic cases D > 1. Using these inequalities we get 


1 1 
and the relationship 
— 1)(t + 1)? rs'(r — 1) 
9) 4rs(t + 1)? + (s — 1)[rst’ 


provides bounds on the quantity N. 


4. The third moment. In this section we will see how the third 
moment of Y/N compares with the third moment of xiy)/N when N 
is given by (8) and when the inequalities D > 1 and 0 < h® < 1 hold. 
If we let d = C, — C, where C, is the coefficient of ¢*/3! in equation (4), 
and C, is the coefficient of #°/3! in equation (5), we get 


d = n,(n, + 2)(n. + + + 2)B3B, + + 2)B.B; 


+ +0n 22 +4) 


Using equations (6) and (7), which were obtained by equating the first 
and second moments of Y/N and x%v)/N, we get 


nines 


Substituting for the constants and simplifying, we obtain 


8(s — 1)(1 + stK)(1 + rstK)’ 
r(r — 1)*s‘(1 + tK)* 


| 
< K < 
6 
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If we proceed as we did with the second moment, we can obtain the 
following inequalities which are valid if D > 1 and if 0 < h’ < 1: 


. a< (s — 1)(t + 2)(st + t + 2)(rst + t + 2)” 
= 2r°(r — + 


(10) 


2 


Thus we see that the first two moments of Y/N and x‘w)/N are equal 
when N is given by (8) and the third moments differ by d where the 
bounds of d are given by (10). 


5. Confidence interval on h?. To obtain an approximate (1 — 2a) 


confidence interval on h? as defined by (1), we proceed as follows: Y is 
approximately distributed as xin), 2:A,/o7 is distributed as x{,,) , and 
the two are independent. Therefore 


‘A 1 

w/o * 
is approximately distributed as Snedecor’s F distribution with N degrees 
of freedom in the numerator and n, degrees of freedom in the denomina- 
tor. Let F, and F, be two values such that 


[50 du = du = a. 
Using these we get 
P(F; <u < F,) = 1 — 2a 


where P(f, < q < f2) = V means the probability that f, < q < f. is 
equal to V. Substituting for uw and making some manipulations, we get 


2(C — F,) 2C — F,) 


where 


(s — 1)A, + As. 
8A, 


Thus the quantity in the brackets determines an approximate (1 — 2a) 
confidence interval on h’. 


6. Empirical investigation of accuracy. If the value of N in (8) is 
used, we have shown that the first two moments of Y/N and xiw)/N are 
equal and, if D > 1 and0O < h’ < 1, the third moments are almost equal 
when r is of appreciable size. We could examine higher moments to see 


C= 


= 

| = 

te 
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how well they fit. The computations, however, become very cumber- 
some when this course is pursued; therefore, it was decided to use 
empirical sampling to investigate the accuracy of the confidence limits 
given in (11) when N is calculated by using the two bounds in (9). 

One hundred thousand random normal deviates (reference [4]) were 
combined to form the observations Y;;, where 


4=1,2,---,7; j=1,2,-::,8; k=1,2,---,t. 


Thus the a; , 6;; , and c;;, are independent observations drawn from 

normal populations whose means are zero and whose variances o2 , 93 , 

and o; are each equal to unity. The value of » was taken as 15 in order 

to eliminate negative values of Y;;, (the range of the random normal 

deviates was from —4.417 to +4.417). This gave a value for h’ of 4/3. 
Eight combinations of r, s, and ¢ were used as follows: 


Total number of 
Case r 8 t observations 
1 2 2 2 8 
2 3 2 2 12 
3 5 2 2 20 
4 10 2 2 40 
5 2 4 2 16 
6 3 4 2 24 
7 5 4 2 40 
8 10 4 2 80 


For each case, 500 sets of observations were drawn from the random 
normal deviates, 500 analyses of variance were made and from each of 
the 500 analyses of variance confidence limits were computed for h’. 
For example in case 4, which consisted of a total of 40 observations, 
500 sets of 40 observations were drawn, and for each set of 40 an analysis 
of variance was run and confidence limits computed for h” by using the 
formula given in (11). Actually two sets of confidence limits were 
calculated for each analysis of variance, one using the upper bound of 
N (when tK = 0) and the other using the lower bound of N (when 
tK = t/(t + 2), given in (9)). Confidence limits were also placed on o: 
as a check to insure that our sample was not too divergent from what 
was expected. 95% confidence intervals were used throughout, and 
since an exact method of setting confidence limits on o- is available, 


re 
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we would expect them to contain o2 95% of the time. A chi-square 
test was run for each case on the number of times out of 500 that the 
calculated limits did not contain 02. In none of the eight cases was the 
chi-square value significant at the 5% level. Thus there was no reason 
to believe that the samples were “bad samples.” The results of the 
sampling are reported below. 


Percentage of 
Percentage of confidence limits containing confidence limits 
Case h? = 4/3 containing o? = 1 
number 
Maximum N Minimum N 

it 94.8 97.6 95.6 
2 91.6 93.6 94.0 
3 96.2 96.8 94.4 
4 95.8 96.6 95.0 
5 91.6 96.8 94.6 
6 91.0 95.6 95.6 
7 89.6 93.0 93.4 
8 92.6 94.6 96.0 


The accuracy of the confidence limits on h’ is satisfactory for the 
cases cited. In fact, the results were exceptionally good considering 
that such small values of r, s, and ¢ were used. Since h? = 4/3 is not in 
the range which was used to calculate N in (9), it appears that the 
confidence limits would be even more accurate if a population had been 
selected in which h” was between zero and one. 


7. Numerical example. We will present an example using data 
which were obtained from the Poultry Department of Oklahoma 
A and M College. Varying numbers of dams were mated at random to 
sires and the 12-week body weight of male offspring was recorded. If 
the data on all the offspring had been used, it would have demanded an 
unequal subclass analysis. In order to avoid this, and since the above 
method of setting confidence limits on h? was discussed only for equal 
subclass numbers, the data were sampled to provide subclass numbers 
which were equal.* The data as finally analyzed consisted of 22 sires 
each mated at random to 6 dams, and the 12-week body weight of 8 
offspring from each mating was used. An analysis of variance was 
made, and the results are given in Table 2. 


*Other minor adjustments were also made on these data such as filling in some missing items, etc. 
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TABLE 2 
Analysis of Variance of Poultry Weights 


Source of Degrees of Mean Expected 
variation freedom square Mean square 
Total 1055 
Between sires 21 0.5629 = A; o2 + toZ + sto? 
Between dams 
within sires 110 0.2055 = A; o2 + toy 
Between full 
siblings 924 0.0924 = A, o 


To calculate confidence limits on 


h? + o;) 
oto 
by using (11), we must calculate C, F, , and F,. By substituting the 
values from the table we get 


C = 2.87. 


To obtain F, and F, we must first evaluate N and n,. By definition, 
nm, = 924, and N will be calculated by using (8). Since we know that 
0 < h® < 1, and since, for these data, we are willing to assume that 
D > 1, it follows that N satisfies the inequality 64 < N < 131. Since 
n, and N are quite large it will make very little difference which of the 
extreme values of N we use. We will use the upper bound to illustrate 
the method and then we will also compute the confidence limits for the 
lower bound of N and compare the two results. If we set 90% confidence 
limits on h’, then a = .05, and F, is the value in Snedecor’s F table which 
corresponds to N degrees of freedom for the numerator (larger mean 
square) and n, degrees of freedom for the denominator (smaller mean 
square). This gives F, = 1.23. To obtain F, , let F* be the tabulated 
F value with n, degrees of freedom for the numerator and N degrees of 
freedom for the denominator. Then F, = 1/F*. This gives F, = 
1/1.26 = 0.79. 
Upon substituting the values into (11) we get 


29 < h’ < .50 


as the 90% confidence limits on h’. If we use the lower bound of N 
instead of the upper, the limits are .26 < h’ < .54. 
We could also calculate confidence limits on h’ by estimating K 
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from the analysis of variance and using this to obtain N. We get 0.0363 
for the estimate of K and 102 as the corresponding value of N. This 
gives .28 < < .51. 
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A PLAN FOR PROGRAMMING ANALYSIS OF VARIANCE 
FOR GENERAL PURPOSE COMPUTERS 
H. O. Harrier 
Towa State College 


1. The necessity for standardization 


The analysis of variance of data arising from an experiment is a 
numerical procedure which is very frequently performed in numerous 
departments at universities and colleges as well as in other research 
centers. In spite of the considerable volume of computation expended 
on this activity most of the work is still carried out on desk computers, 
and this is even true of many centers at which the services of a high 
speed general purpose computer are available. The reason for this is 
undoubtedly the great variety of experimental designs, each of which 
gives rise to a different type of analysis of variance each applied to a 
small body of data. There is no difficulty in setting up and testing 
suitable programs every time data from a new design are ready for 
analysis, but in so far as the time and effort of doing this is usually 
much greater than the effort of completing the analysis of variance on a 
desk computer, there is clearly no point in enlisting the high speed 
machines.* The question, however, arises whether it is not possible to 
so adapt the analyses of variance for the diverse designs that they can 
all be covered by a standard computing program. Such a standard 
program would be set up and tested once and for all, and would then be 
available for the analysis of variance of data from any design, it being 
only necessary to convey minor program instructions pertaining to each 
particular design. Whilst such a standard program is clearly not possible 
for the analysis of any set of data an attempt is here made to at least 
cover the majority of designs arising in experimental work. Because 
of the above reasons standardization of the program is here regarded 
as of prime importance. It is fully realized that when dealing with a 
particular design on a particular machine there may well be alternative 
programs resulting in shorter computing times. A number of such 
particular programs have been published in the past. (See List of 
References.) The basic analysis to which that of other designs will 
here be reduced is that of a general factorial experiment which will be 
discussed in the next section. 


*The comparative clerical labor of preparing the data for input into the respective machines, 
although by no means negligible, is not discussed here as this depends on details of the organization 
of the computing center. 


110 


‘ 
\ 
| 
4 
= 
Af 
oH 
qi 
{4 
j 


PROGRAMMING FOR COMPUTERS 111 


2. The program for the k-factor experiment based on a special operator 
calculus 


For convenience we confine ourselves to a k = 3 factor experiment. 
Let x,;; denote the experimental result from ith level of factor ‘T’, 
ith level of factor ‘I’ and jth level of factor ‘J’. The symbols T, I and 
J will also denote the number of levels for each factor so that 
t=1,2,---,T,i=1,2,---,Zandj = 1,2,---,J. The complete 
analysis of variance of the T X I X J results into its 2* — 1 = 7 com- 
ponents is shown in Table 1. 


TABLE 1 
Analysis of Variance for 3-Factor Experiment 
Component Degrees of freedom 
(T 1) 
I (I — 1) 
J (J — 1) 
(T - 
TXJ (T 1) 
IxJ (I — 1)V 1) 


For the corresponding sums of squares we shall require the familiar 
notation for group totals, viz. 


Xai = and likewise X,.; X,;. 
(1) 


T 


X5= and likewise X.;, X,... 


t=1 


The sums of squares can be obtained by repeated application of the 
following operators which will be explained in terms of factor T’. 


Operator 2, = sum over all levels of ¢ = 1, 2, --- 7 whilst 
keeping the other subscripts constant 
Operator D, = multiply all items by 7 and subtract the 


result of 2, from all items 

(2) Operator ( )? = form the sum of the squares of the items in- 
side the brackets and divide by the number 
of items. 
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For example if we apply the first two operators to the original set of 
results z,,; we have 


= = the total for the 7, 7 combination of 

(3) factors J and J 
= T tu; — J the deviate of z,;; from the 7, 7 mean 
= T(z.;; — (multiplied by T. 


The above simple operations represent the first two lines in the schedule 
of operations shown in Table 2 which gives complete formulas for the 
totals and deviates resulting from the sequence of operations 2,D,2; 
D;2,;D; applied to the data z,;; . The seven sets of deviates finally 
reached in lines 9 to 15 are finally subjected to the Mean Square Opera- 
tion (_ )? and the results are the ‘Sums of Squares of Deviations’ (all 
multiplied by TIJ) for the seven components of variance shown in 
the fifth column of Table 2. Table 3 (below) illustrates these operations 


TABLE 2 
Schedule of Operations for Three Factor Analysis of Variance 
Deviates 
Applied to used for 
Line | Oper- | values in Will form totals or deviates analysis of 
ator lines variance 
components 
1 Input Leis 
2 1 X 
3 1D, 1 and.2 — Xs; 
4 2 
5 3 — 
6 D; 2 and 4 1X3; 
7 D; 3 and 5 — 1X 65 — TXe.5 
8 X 
9 5 7’ 
10 6 I 
ll 2; — — TX:.. + X 
12 D; 4 and 8 J 
13 D; 5 and 9 TIX 4.3 — + X 
14 D; 6and10 | IJX.4; — JX. — 1X4. +X IxJ 
15 D; Tand1l | — [JX 4; — |TXIXS 
— + 1X4. + 7TX:.. —X 
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TABLE 3b 
Analysis of Variance of Data, Table 3a 


Components TIJ (S.o.squares) 
6 
I 78 
J 4 
186 
182 
182 
TREKS 46 


with the help of a simple example in which T = 3, J = 3 and J = 2. 
It is these figures, i.e. the number of levels in each factor, that may 
vary from experiment to experiment and need to be conveyed to the 
machine. If the number of factors, k, exceeds 3 the sequence of operators 
> D --- is extended as part of the program. A decision may have to 
be made on the maximum number of factors (k) to be covered by the 
standard program. In this connection it will be observed from Table 2 
that when the operations of lines 1 to 15 are completed, the deviates 
formed in lines 1 to 7 are no longer required; if other operations, cor- 
responding to further factors, require additional storage space, the stores 
containing lines 1 to 7 could be cleared and used for subsequent deviates. 
If this saving of storage is carried to its full consequences the minimum 
storage space required is that for the final set of deviates to which the 
mean square operations (_ )* are applied. To give an example, take a 
rather large experiment of k = 4 factors each at p = 6 levels, ie. an 
experiment consisting of p* = 6* = 1296 observations. By the time 
the pair of © and D operators have been applied for the last (4th) 
factor the number of deviates to be stored has risen to 


k 


(4) = (p + = 7' = 2401. 


For the particular case of a 2‘ experiment in which all factors are at 
2 levels (T = I = J = 2) the operators 2, D, , --- , correspond to the 
well known procedure of taking sums and differences with regard to 
each factor in turn and the deviates formed in the 4th column of Table 2 
are in fact the well known ‘Effects’ (i.e. main effects, 1st order and 2nd 
order interaction effects) of the 2° analysis. The present procedure 
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may therefore be regarded as a generalization of this well known 2" | 
analysis to any number of factors with any number of levels, and it 
makes provision for a complete computation of sets of ‘effect values’ in 
terms of sets of deviates. In many cases it will be worthwhile to arrange 
for the machine to print out these deviates as they usually provide a 
better insight into the nature of the data than do the mere sums of 
their squares. In particular this program facilitates studies on the 
fundamental assumptions of Analysis of Variance such as additivity and 
normality as well as simultaneous studies of various variate transforma- 
tions. 

The detailed method of programming the operators 2, , D, , --- , will 
depend on the facilities of the particular machine. 


3. The reduction of the analysis of other designs to the factorial program 


3.1 Analysis by standard factorial program. The k-factor experiment 
as described in the preceding section will result in an analysis of variance 
as shown (for k = 3) in Table 1, and in the general case consisting of 
2* — 1 components of variance. Whilst such designs are not uncommon 
they are certainly not predominant and in the majority of cases the 
design employed will not be of this kind. For example, the simplest of 
all analysis of variance situations is an experiment in which T' experi- 
mental groups are each replicated J times, and x,; denotes the result 
from the 7th replicate in the 7th treatment group. Since the J replicates 
in each treatment group are completely unrelated they do not constitute 
a factor. Nevertheless we may analyze the data z,; formally as if 
they came from a 2-factor experiment with ‘factors’ T and I and thus 
obtain from the standard computing procedure the sums of squares 
(multiplied by TJ) for the components T, J and T X I. To complete 
the proper analysis of variance for this experiment the program will 
contain a code for a ‘Summary of Components’ which would be appro- 
priate for this particular design and would read as follows:— 


Summary:— I + (T X I) = Error (within treatments). 


Again if we have a ‘Split Plot’ design with T main treatments and I 
subtreatments arranged in J blocks we can obtain the appropriate 
analysis of variance by performing first the standard factorial analysis 
and then adding the summary instruction:— 


Summary:— Error (b) = (J X J) + (T XI X J) 


The general principle of this procedure is therefore to perform first a 
formal factorial analysis and then pool certain components in accordance 
with summary instructions which specifically apply to the particular 
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_ design. This would certainly be a wasteful procedure for a desk com- 
puter but is convenient for a high speed computer. In Table 4 (at end) 
we list examples of such experiments which can be dealt with by this 
method. We give for each design 


(a) The factors which would be associated with the standard factors 
of the factorial analysis, 

(b) The summary instructions for the design, 

(c) The final analysis of variance in terms of the main effects and 
interactions of the factorial analysis. 


It will be seen that a considerable diversity of designs can be covered 
in this manner. 


3.2 Analysis involving ‘rearrangements’. Designs not covered by the 
principle of reduction used in the preceding section are those involving 
a higher degree of balance, i.e. a higher degree of restriction in the 
randomization such as the Latin Square, Youden Square, Lattices, 
Incomplete Randomized Blocks and Confounded Designs. In order to 
deal with such situations only one intermediate operation must be 
introduced. This operation, called ‘rearrangement’ will be described 
in terms of the Latin Square below:— Let us denote by z;;,,) the record 
for the ‘plot’ in the 7th ‘Row’ and the jth ‘Column’ to which the éth 
treatment has been applied. There are only 7’ records and these are 
first conveyed to the machine store arranged by the two factors ‘Rows’ 
t= 1,2--- JZ = T and ‘Columns’ j = 1,2, --- , J = T. Along with 
Zija) We convey the treatment code ¢ and this is carried through all 
the D operations, but ignored in the = operations so that the results of 
the sequence 2; D; 2; D; are the quantities shown below:— 


Operation Yields values 
2; 
D; 2; T X;. 


At this stage the latter quantities D; D;x;;:,, will be ‘rearranged’ and 
conveyed to a store arranged by (say) rows, z, and treatments, ¢. To 
these quantities u;, = 7? 2;;..)5 — TX.; — TX;. + X.. will then be 
applied the operators 2/ and D{ in which the’ signifies that the operation 
occurs after rearrangement. The results of these operations are as 
follows: 
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Operation Yields values 
= — TX — TX + TX = T(TX, — X) 
= — TX, — TX;, — TX, + 2X) . 


In terms of these operations the final analysis of variance can clearly 
be written as follows: 


Rows (D; =I 
Columns (2, =J 
(7) Treatments (2! D; Dj)’ 
7 (Dt D, 


Here the symbol T~'(D; D; D;)* means:— To the x;; arranged by the 
1st factors (Rows Columns) the operators D; and then D, were applied, 
the deviates thereby resulting were rearranged by the 2nd factors 
(Rows 7 and Treatments ¢) and to these were applied operation D; and 
the mean square operation (_ )’, finally the answer was divided by 
T.* With the help of this additional operation a wider variety of 
experimental designs can be analyzed in terms of the standard factorial 
program. With certain designs it may be necessary to convey a double 
subscript code required for a subsequent rearrangement. Some of 
these designs are set out in Table 5 (below) in which are shown for 
each design the ‘factors’ before and after rearrangement and the final 
analysis of variance formulas in terms of the operators 2, D and 2’, D’. 
For certain designs it is convenient to also employ a ‘Total Sum of 
Squares of Deviations Operation’ (D)’ which is applied to the total 
N =T X number of observations as soon as they are fed 


into the machine. This operator is defined by 


8) (Dy = (Wz XY. 


It should be noted that the ‘rearrangement operation’ may be 
applied at different stages. In the five examples shown in Table 5 it 


*It will be seen that the only purpose for the ‘rearrangement’ of the zi;(t) by the factors J and T is 
to facilitate the subsequent summation ©,’ of the zi;(¢) over i for constant values of t and subsequent 
operation D;’. On most machines, the lengthy operation involved here would be the internal scanning 
of the stored code ‘t’.. With rearrangement one scanning would be required whilst if operations 2,’ and 
Dj’ were carried out directly without rearrangement two scannings will normally be required. 
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‘s applied as follows:— 


Design:— Rearrangement applied to quantities:— 
Latin Square D; D; x; 
Incomplete Blocks Dy 
Youden Square D, Di 
Triple Lattice D; 
Factorial 2* (confounded) 


In the last design, therefore, the data are fed into the machine 
arranged by r (Replicates) and treatment combinations (abc). They 
will be immediately rearranged by the 3 new factors, i.e. by r (Replicute), 
g (Treatment Group) and ¢ (No. of Treatment Combination within 
Group). If (say) ABC is the confounded effect then this arrangement 
of treatment combinations is as follows: 


No. of Treatment in Group 


g=1 1 ab ac be 
Treatment Group 
g=2 abc c b a 


The subsequent operations are then completely standard. For the 
arrangement by the Ist factors we have an ordinary (non-confounded) 
analysis like that in the last example of Table 4, regarding Blocks as 
Replicates, for arrangement by the second set of factors we have a 
standard split plot analysis. The two analyses should check on the 
following components: 


First Factors Analysis:— Second Factors Analysis:— 
(A)+(B)+(C)+(AB)+(AC)+(BC) = (T) + (T X G) 
(ABC) = G 
(9) R = R 
‘Error’ = Error (a) + Error (6) 


Finally we should mention that the analysis of variance for the 
‘Balanced Incomplete Blocks’ as set out in Table 5 gives the one for 
treatment comparisons without recovery of inter block information 
whilst that for triple lattices is the one providing the weights for treat- 
ment-mean comparisons. There is no inherent difficulty in program- 
ming for the alternative analysis in either design. 
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4. A universal missing plot formula 


It is well known that an orthogonal design in which one or several 
results are missing can be essentially reduced to standard orthogonal 
analysis with the help of what are known as ‘missing plot formulas’. 
Most of these formulas are for a single missing plot and have arisen by 
minimizing the appropriate ‘Error Sum of Squares’ as a function of 
the unknown missing value. The resulting ‘missing plot formulas’ 
differ from design to design and the use of these formulas for high-speed 
computing would require a different missing plot program for every 
design. In order to avoid this we offer here a universal program appli- 
cable to all designs:— 

Assume then that the standard program for the orthogonal analysis 
of variance (i.e. without missing values) is already programmed and 
let Q denote the appropriate ‘Error Sum of Squares’. Denote by 
@y , a, and a, any three trial values for the missing plot, unity apart, so 
that a; = a;_, + 1, and by Q) , Q, and Q, the corresponding ‘Error Sum 
of Squares’ computed from the observed results supplemented by 
@) , a, and a, respectively. The universal formula for the missing 
value is then given by 


(10) a = a, + (Qo — Q2)/2(Qo — 2Q; + Qz2) 


which is in fact that for the a-scale value at which the parabola through 
the three points (a) , Qo), (a; , Q:) and (a, , Q.) attains its minimum. 
It will be noted that the only addition to the standard program for analy- 
sis of variance is formula (10), if one is satisfied with the ‘approximate 
analysis of variance tests’, i.e. with supplementing the data by ‘the 
missing value’ a, computed from (10), and using the standard analysis 
of variance program with the error degrees of freedom reduced by 1. 
By a similar method it is possible to program the ‘exact’ missing plot 
analysis resulting from the likelihood ratio principle. Likewise it is 
possible to cover the case of several missing plots using (10) and the 
customary iterative procedure. 


5. The limitations of the present program 


It will be noted that certain Analysis of Variance situations are not 
covered by the present program. Notably these are cases of unbalanced 
non-orthogonal designs. In such cases the complete analysis usually 
requires a basic Least Squares model in which all levels of all factors 
are represented by ‘unknown effect constants’ to be fitted to the observed 
data. Such a computation can, of course, be covered by a general 
‘Least Squares’ program comprising the formation of the ‘normal 
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equations’ and their solution. This part of the analysis (i.e. the estima- 
tion of effect constants) is comparatively straightforward and follows 
the principle set out recently by Tocher (1952). However, the final 
analysis of variance tests depend on a hierarchical sequence in which 
the factors are to be arranged. 

Finally we should like to stress that the basic principle of the present 
approach is to reduce the analysis of a given design, by minor additional 
program instructions, to the standard program. Indeed, for the more 
complex designs here discussed it may be argued that the ‘standard’ 
part of the program consists merely of standard ‘subroutines’ in the 
form of the operators D, Z,( )* and the rearrangement. Whilst these 
latter could be used for every design a different ‘summary’ or ‘steering 
program’ would be needed to organize the operations in the correct 
sequence. For the bulk of the simpler designs, however, most of the 
program would consist of the standard subroutines. In the case of a 
very complex experiment it would, therefore, defeat the principle of 
this approach if very elaborate program additions are found to be 
required to reduce the analysis to the standard. In such cases it will be 
wiser to revert to a known desk computer layout and analysis, partic- 
ularly if it is unlikely that the design will be used frequently. 
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NOTE ON FITTING THE MULTI-HIT SURVIVAL CURVE * 
Joan M. GuRIAN 


Division of Biological and Medical Research 
Argonne National Laboratory, Lemont, Ill. 


The expected proportions (S) of a population of micro-organi ms 
surviving irradiation dose x has been described by 


Kimball [1] proposes a graphical method for evaluating the consts ts 
k and n in this expression. He defines his dependent variable 


u; = n(1 — S,) 7=1,2,---p 


where S; is the observed proportion surviving, and minimizes he 
quantity 


P 


V = Dl; —e 

t=1 

with respect to k, by a graphical procedure. This method has the 
advantages of being simple, requiring no arbitrary approximations, and 
permitting the use of all data. However, application of the transforma- 
tion used in this procedure, without making corresponding adjustments 
of the weights, will give the larger S; unduly heavy weight. This 
results in a poor fit at the lower values of S; , if S; varies by several 
orders of magnitude. 

Survival experiments are usually designed to give equal percentage 
plating error for all S; or constant variance o” for all log S;. If u; = 
In (1 — S,) is used as the dependent variable the sampling varia :ce 
of the dependent variable is no longer constant, but is approximately 
equal to o[S,/(1 — S,)]’ [2]. The theoretical weight for wu; is then 
[(1 — S,)/S,}’. Errors in dose, dilution, ete., are often propagated as S; 
decreases [3], hence these effects must also be taken into consideration 
in the choice of the weighting function. Ideally the appropriate weights 
should be derived empirically from a properly designed experiment; but 
prior knowledge may be used to give a rough estimate of the weights, 
which may prove adequate for the evaluation of the parameters. 


*This work was performed under the auspices of the U. S. Atomic Energy Commission, 
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TABLE 1 
Formulae for Approximating n and k 
Quantity Kimball’s formulae Weighted Kimball’s formulae 
k Graphical estimate Graphical estimate 
UY; > WUY; 
n 
7 ww? 
i=1 i=1 
v; = In(1 — e**’) v; = n(l — 
V Yui uy, —n way, 
i=1 i=1 i=l 
AB — C’ AB — 
AB — AB — 
where: 
p—2 
A wy? 
i=1 
2 P 2 
C n ( ) h ) 
TABLE 2 
n, k and Their Errors Computed 
Quantity Kimball’s Weighted Kimball’s 
Method Method 
n 2.220 2.654 
k .525 .600 
8n .140 .199 
-033 -028 
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Table 1 contains in Column 1 Kimball’s formulae used in the approxi- 
mations of n and k, and the sampling errors of their final values, and 
in Column 2 the corresponding modified formulae. The parameters k 
and n are determined by the same method of iterative approximations as 
in the original procedure. The only increase in labor and complexity is 
that due to the presence of differing weights. 

Table 2 contains the constants and their errors as computed by the 
two procedures, using the data of Dr. S. Pomper on the survival of 
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Saccharomyces cerevisiae, which were employed by Kimball. The 
weighting function used is w; = (1 — S;)*. It is based on the assump- 
tion that in this experiment errors in dilution, etc., increase: with de- 
creasing S; in such a manner, that they offset the decreasing theoretical 
sampling error, so that the variance of S; is constant. 

Figure 1 shows the data as fitted by the two methods. 
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QUALITY EVALUATION BY NUMERICAL AND SUBJECTIVE 
METHODS, WITH APPLICATION TO DRIED VENEER 


W. G. Kauman, J. W. Gorrstern, anp D. LANTICAN 


Officers of the Division of Forest Products, C.S.I.R.0., 
Instructor in Wood Technology, University of the Philippines 


I. INTRODUCTION 


Very different schemes have often been proposed for the evaluation 
of perceptual entities. Advocates of physical measurement, either by 
instruments or by unambiguously defined scoring systems, claim that 
any method of subjective evaluation will suffer from variability between 
and within observers, and from difficulties in training observers other 
than by personal tuition. Those advocating subjective evaluation 
schemes retort that for many of the quantities in question, physical 
measurement would be so involved as to be impracticable, and that the 
results, though reproducible, would often be meaningless for the intended 
purpose. 

If a simple, straightforward method of physical measurement is 
available, it would clearly be a waste of time and effort to substitute 
subjective evaluation, as, for instance, for the evaluation of weight or of 
electric current intensity. On the other hand, quantities or properties 
whose physical background is very complex or poorly understood, but 
which form well-defined sensory concepts, such as taste flavors or the 
notion of personal comfort, would appear more amenable to subjective 
evaluation than physical measurement. The use of subjective methods 
for taste testing is so widespread as to require no comment. With 
regard to personal comfort Muncey (1954) showed, for instance, that 
the subjective comfort—rated on a linear scale—may be used as satis- 
factorily as foot temperature to assess the ‘‘thermal comfort”’ of different 
floor materials. 

In many instances, quality judgments are based on a number of 
more or less understood physical variables, the assessment depending 
to a large degree on the weighting and integration of the individual 
factors. Perceptive concepts of this type are considered by some workers 
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to be of the nature of a “‘Gestalt’’*; for example the “firmness’’ of rubber 
samples (Scott Blair and Coppen, 1942) or the “tackiness” of glues 
(Strasburger, 1953). The perception of faults in veneer would also 
appear to be of ‘“‘Gestalt’’ character, the observer forming his subjective 
judgment of the quality by a perceptively integrated impression rather 
than by conscious observation and summation of individual defects 
present. 

Objective numerical evaluation of ‘Gestalt’ properties requires 
analysis into dimensionally simple components which are amenable to 
physical measurement. Scott Blair and Coppen (1942) consider that 
such an analysis by the physicist is “apt to destroy much of the original 
concept’’, and, at best, ‘‘can offer but a blurred picture of the subjective 
original’. 

Any gain in consistency and reproducibility of results provided by 
a numerical method based on measurements may thus be vitiated by 
the experimenter losing sight of the original property, quite apart from 
the considerable complexity of procedure which is often entailed. 

On the other hand, if a subjective evaluation method is chosen, the 
experimenter must appreciate the importance of the mental attitude 
(Scott Blair and Coppen, 1940) and strain and fatigue of the observers 
(Harper, 1952; Hopkinson, 1952; Wyatt and Langdon, 1932), and of 
proper randomization of samples (Harrison and Elder, 1950; Hopkins, 
1950; Sheppard, 1953). 

It is also necessary to select carefully the actual evaluation technique 
to be used. Bradley (1953) advocates wider use of ranking methods, 
but Hopkins (1950) considers this procedure of advantage only when 
samples are few in number or can be grouped in small sets. For veneer 
quality evaluation, ranking methods would be highly impracticable 
because of the large number of samples usually involved, as would also 
be paired comparisons which are commonly used in tasting techniques 
(Bradley, 1953). Scoring techniques which allow independent assess- 
ment of individual specimens and facilitate standardization of quality 
ratings would appear clearly preferable in this case. 

When selecting an evaluation method for a particular application, 
it should be appreciated that the method need not yield 100 per cent 
accuracy, nor is it necessary for every single judgment to be correct. 
It is sufficient ‘‘that the mean of a convenient number of judgments 
should not drift, and that the variance of the judgments should be 
small in relation to the differences in which (the) experiment is inter- 
ested” (Hopkinson, 1953). 


*The term “Gestalt” was introduced into psychology by C. Ehrenfels in 1890 to indicate the char- 
acter of perception as a unity (Encycl. Britt. 1947). 
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It would therefore appear that, a priori, numerical systems have no 
particular advantage over subjective evaluation methods. The selection 
of a method is largely a matter of defining the desired maximum level 
of variance, and having decided the permissible variance, ef considering 
which evaluation method requires least time to produce results of the 
desired accuracy. The actual choice of a method will, of course, be 
influenced by additional arguments, such as the cost of apparatus, 
availability of trained personnel, relative difficulty of training observers, 
and so on. In most practical cases, it will be necessary to base the final 
choice on the statistical analysis of an exploration experiment such as 
the one described in the present paper. 


2. SCHEMES FOR EVALUATION OF THE QUALITY OF DRIED VENEER 


The success or failure of experimental veneer drying schedules is 
usually assessed in terms of the “dried quality” of the veneer produced. 
It is therefore of great importance to evaluate quality by a simple, 
unambiguous method which should give accurate and reproducible 
results. 

The “‘dried quality” of a sheet of veneer may be defined as the relative 
freedom from drying degrade*, with careful specification of the types of 
degrade considered. Evaluation of dried quality is thus largely a 
problem of assessing degrade. 

The assessment may be made either by a method of wn 
measurement of the amount and severity of the different degrade 
types present, or by a system involving subjective judgment (Ellwood 
1952). Quality evaluation schemes described in the present paper have 
been developed for application to “ash” type eucalypt 38 X 38 in. 
veneer sample sheets of 1/16 in. thickness (green dimensions). Although 
the great majority of species used for veneers can be dried with little 
degrade using quite a wide range of drying conditions, in the ‘‘ash” euca- 
lypts drying conditions are extremely critical and the remarkable amount 
of degrade which may occur makes assessment of effects difficult and 
necessitates the use of special evaluation methods. Application of the 
schemes to veneer of other species, dimensions or thickness offers, how- 
ever, no major obstacles in principle and will be discussed in a later 
section. 

The schemes were conceived as a compromise between the require- 
ments of conformity to commercial standards and ease of application in 
experimental research work. They are intended for research purposes 


*Defects such as checks, splits or buckling arising in timber during drying are commonly termed 
“‘degrade’’. 
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only, and not to supplant grading schemes commonly used in commerical 
practice. 

In “ash” eucalypt veneer, it is convenient to distinguish four principal 
types of drying degrade, comprising, in order of decreasing importance: 


(i) Through-checking 
(ii) End splitting 

(iii) Face checking 

(iv) Buckling 


Buckling in these species is largely due to differences in collapse* 
in early and late wood zones and can usually be reduced to a marked 
degree by a steaming treatment known as “reconditioning”; it is there- 
fore considered a minor form of degrade. 

Face checks (i.e. checks in the tight face not penetrating through the 
whole thickness of veneer) often close up during the manufactuve of 
plywood and become inconspicuous; they too are considered minor 
degrade. 

End splits are usually a moderately severe form of degrade, but may 
be rated more or less severely according to their constitution. End 
splits with flat and adjacent edges may become practically invisible 
when the veneer is made up into plywood, whereas end splits with curly, 
widely separated or overlapping edges cause severe and often irreparable 
defects. 

Through-checks (i.e. checks penetrating through the whole thickness 
of veneer) generally constitute a severe form of degrade, except when 
they are of the nature of extremely fine hairline checks. Through-checks 
which have an appreciable width cause the glue to well up and spread 
on the surface during pressing, thus spoiling the appearance of the 
plywood. If the through-checks have curly edges, they cause severe 
overlapping defects. 

Although through-checking, as well as face checking, is usually of 
minor importance in veneer from most species, both these degrade 
types can cause very severe degrade in the “ash” type eucalypts. 

The amount of degrade of each type present is classified into the 
following six ‘‘severity classes’’: 

None 

Very slight 
Slight 
Moderate 
Severe 
Very severe 


wnre © 


*Collapse is an abnormal form of shrinkage caused by caving in or ‘‘collapse’’ of the cell lumina. 
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The evaluation schemes are intended to assess only degrade caused 
by the drying treatment. Any defects present prior to drying, and 
drying degrade clearly caused by blemishes such as knots inherent in 
the timber are therefore disregarded. 

The final result of the evaluation is termed the “quality rating” of 
a veneer sheet. It is obtained by arithmetical or mental summation of 
the severity classes (according to the scheme used), and provides a 
quantitative measure of dried quality. 


2.1 Numerical Evaluation Scheme 


In the numerical scheme the allocation of severity classes for checks 
and splits is based on the measured total combined length of each type 
present, as shown below: 


Total combined length of 
all checks or splits of Severity 
the type being assessed* class 


0 0 
Under 2 in. 1 
2 in. to 6 in. 2 
6 in. to 12 in. 3 
12 in. to 24 in. 4 
Greater than 24 in. 5 


*For 1/16 in. thick 38 X 38 in. veneer sheets. 


No satisfactory method of measurement could be found for buckling, 
and classification of this degrade type is therefore based on an estimate 
of the amount and average amplitude of the corrugations in the veneer 
sheet. The amount of subjectivity introduced by this procedure is 
very small since buckling makes only a minor contribution to the quality 
rating. 

Each degrade type is associated with a “numerical weight”? which 
indicates its relative importance for the overall quality assessment. 

The “degrade rating” of a degrade type is obtained by multiplying 
the severity class by the numerical weight. Table 1 summarizes degrade 
types with their respective numerical weights, and also gives the maxi- 
mum degrade rating of each type, obtained by multiplying the numerical 
weight by the highest severity class, i.e. by 5. 

In samples where several sub-types of a particular degrade type are 
present, the total severity class must not exceed 5 and is allotted accord- 
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ing to the total combined lengths of the checks or splits of all sub-types, 
but giving full effect to the sub-type having the greatest numerical 
weight. 

If any type of degrade is localized (i.e. contained in not more than 
two regions whose total combined area does not exceed approximately 
} of the area of the veneer sheet), the corresponding degrade rating is 
multiplied by 3, rounding off to the nearest integer. 

The “quality rating” of a veneer sheet in the numerical scheme is 
the sum of the degrade ratings of all degrade types present. A perfect 
sheet completely free from degrade will thus have a rating of 0, whereas 
a very bad sheet containing a maximum of each degrade type will have 
a rating of 50. ; 


TABLE 1 
Numerical Weights and Maximum Degrade Ratings of Different 
Degrade Types 
Maximum 
Degrade Degrade sub-types Numerical| degrade 
types weight rating 
Buckling _— 1 5 
Face checking — 1 5 
End End splits with flat and adjacent edges. 2 
splitting 
End splits with slightly curled and/or 
widely separated edges. 3 
End splits with severely curled and/or 
overlapping edges. 4 20 
-Through- Hairline through-checks 2 
checking 
Through-checks of width less than 1 mm., 
except hairline checks. 3 
Through-checks of width greater than 1 
mm and/or with curled or overlapping 
edges. 4 20 
Maximum quality rating* 50 


*Quality rating 50—'‘‘very bad” veneer sheet with maximum degrade. 
Quality rating O—‘‘excellent’”’ sheet free from degrade. 
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2.2 Subjective Evaluation Scheme 


In the subjective scheme, degrade is classified into severity classes 
by visual estimation. In the case of buckling this is again based on an 
estimate of the amount and average amplitude of the corrugations 
present, whereas checks and splits are judged on the basis of the number 
present, the distribution over the veneer face and the average length 
and width of individual checks or splits. As there are no “numerical 
weights” to be taken into account, the “‘degrade rating” in the subjective 
scheme is synonymous with the “severity class’. 

The quality of the sheet is assigned to one of the following nine 
“quality ratings”: 


Excellent 
Very good 
Good 
Very fair 
Fair 

Poor 
Very poor 
Bad 

Very bad 


The quality rating may be obtained either by mental integration of 
the previously established severity classes or, if no record of severity 
classes is required, by direct subjective judgment based on visual inspec- 
tion of the sheet. The latter of these methods should, however, onl,’ 
be used after some experience has been gained with the former. 


CONOR © 


3. EXPERIMENTAL PROCEDURE 


Twenty sheets of 1/16 in. thick veneer (green size 38 X 38 in.) 
from five trees of alpine ash (Eucalyptus gigantea Hook f.) were selected 
from material included in a study of the mechanical drying of “ash” 
type eucalypt veneer (Kauman, Gottstein and Lantican, 1956). The 
sheets had been selected to include, as far as possible, all types and 
severities of degrade and quality ratings. 

Three observers (A, B and C) each carried out two quality evaluations 
on these twenty veneer sheets, using first the subjective and then the 
numerical scheme. Several days were allowed to pass between the first 
and second evaluation by each observer, and the order of the sheets was 
changed between tests. No observer was permitted to see the others’ 
results, or to be present during the others’ tests, before he had completed 
all his own tests. ; 
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The procedure for the subjective scheme was first to assess the 
quality rating, and then to estimate the severity classes of the various 
degrade types present. In this manner it was hoped to obtain quality 
ratings which were independent of the individual severity classes. 

Observer A was most experienced, having completed a quality 
evaluation of some 340 sheets of veneer by both schemes immediately 
prior to the present experiment; B had acted as recorder to A during 
this evaluation, but had not himself carried out evaluations, whereas 
C had no previous knowledge of the schemes and carried out his evalua- 
tion after a brief training period. Observations were recorded by an 
assistant, and the observers did not consciously remember previous 
results during their second evaluation tests. 


TABLE 2 
Results of Numerical Quality Evaluation 


Quality Ratings 
Sheet Observer A Observer B Observer C Mean 
no. 
Test Test Test Test Test Test 
No. 1 No. 2 No. 1 No. 2 No. 1 No. 2 
1 22 15 21 24 19 11 18.7 
2 31 23 24 21 29 33 26.8 
3 34 37 32 35 40 37 35.8 
4 44 44 41 44 42 40 42.5 
5 10 6 13 19 8 7 10.5 
6 31 32 34 39 34 25 32.5 
7 42 39 39 41 25 30 36.0 
8 28 19 23 27 13 24 22.3 
9 29 29 27 27 28 25 27.5 
10 39 37 29 35 34 36 35.0 
ll 21 20 17 25 16 23 20.3 
12 42 44 44 40 41 ~ 38 41.5 
13 40 39 34 37 34 31 35.8 
14 1l 9 9 12 7 8 9.3 
15 38 37 30 36 44 40 37.5 
16 9 13 14 25 ll 12 14.0 
17 20 16 19 19 6 15 15.8 
18 28 31 30 31 24 37 30.2 
19 28 24 25 30 25 22 25.7 
20 15 8 11 19 9 10 12.0 


= i 
= 
Mean 28.1 26.1 25.8 29.3 24.4 25.2 | 
4 
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4. RESULTS 


Quality ratings evaluated by the numerical and the subjective 
method are given in Tables 2 and 3. Table 4 summarizes the mean 
degrade ratings for all degrade types evaluated by both schemes. 
Detailed evaluation results of degrade types by the different observers 
are not reproduced here. 


TABLE 3 
Results of Subjective Quality Evaluation 


Quality Ratings 
Sheet Observer A Observer B Observer C Mean 
no. 
Test Test Test Test Test Test 
No. 1 No. 2 No. 1 No. 2 No. 1 No. 2 
1 3 2 3 4 3 4 3.2 
2 7 7 5 4.5 7 5 5.9 
3 6 5 5 5 5 6 5.3 
4 7 7 8 8 7 6 7.2 
5 1 1 2 2 3 2 1.8 
6 5 5 5 6 5 5 5.2 
7 5 6 6 6 5 5 5.5 
8 3 3 5 4 4 4 3.8 
9 4.5 4 5 5 5 5 4.8 
10 6 6 7 6 7 7 6.5 
11 5 4 + 4 4 5 4.3 
12 8 7 7 8 8 6 7.3 
13 5 6 7 7 5 5 5.8 
14 1 1 2 2 2 3 1.8 
15 z 8 7 7 7 8 7.3 
16 1 3 3 3 4 4 3.0 
17 4 3 3 3 3 4 3.3 
18 6 4 6 5 5 5 §.2 
19 5 4 5 5 7 5 5.2 
20 3 2 3 4 2 2 2.7 
Mean 4.62 4.40 4.90 4.92 4.90 4.80 
Subjective Quality Ratings: 
Excellent 0 Poor 5 
Very good 1 Very poor 6 
Good 2 Bad 7 
Very fair 3 Very bad 8 


Fair 4 
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TABLE 4 


Mean Degrade Ratings 


| 
Sheet | Numerical scheme | Subjective scheme 

no. | 

| FC. | ES. | | B. | F.C. E.S. rz. 
1 2.0 3.2 | 65 | 6.4 23 1 oA 2.0 2.0 
2/45 | 28 | 47/43 |] 48 | 22 | 
3 2.2 3.5 | 16.7 | 14.5 2.4 2.8 3.6 3.5 
4 1.5 4.0 | 18.3 18.7 1.8 3.3 4.0 4.4 
5 1.7 1.4 5.2 2.3 2.0 0.7 1.6 0.5 
6 1.8 3.2 15.8 11.7 | 2.8 3.3 2.8 
7 2.7 8:7 15.8 | 14.0 3.2 2.8 3.1 3.4 
8 2.0 2.7 9.2 8.5 2.3 2.1 2.5 2.3 
9 2.3 3.7 13.0 8.5 2.4 2.5 3.2 2.6 
10 1.5 3.5 | 16.5 | 13.5 1.6 3.3 4.2 4.2 
11 3.2 2.3 7.7 7.2 4.2 2.1 he 3 2.0 
12 29 4.0 | 16.5 18.2 3.2 3.7 3.5 4.3 
13 $3 4.0 | 17.0 12.3 2.2 3.4 3.6 3:3 
14 i 1.2 4.7 1.8 | 1.8 0.5 1.7 0.7 
15 i) 11.3 | 20.0 | 26 3.2 3.2 | 4.7 
16 | 2.8 2.5 1.8 2.0 1.6 
17 27 |. 30 4.3 6.8 2.8 1:3 1.8 2.0 
18 | 48.7 | |. 3.2 3.7 3.5 
19 $0 3.7 10.3 | i 2:8 3.2 2.8 
20 15 | 28 2.2 | 5.5 | 1.4 | 2.2 0.7 2.3 


Abbreviations: 
B. = Buckling 
F.C. = Face checking 
E.S. = End splitting 
7T.C. = Through-checking. 


Figures 1 and 2 illustrate typical veneer sheets of “good” and “very 
bad” quality, whereas Figure 3 shows the appearance of severe through- 
checking viewed by reflected and transmitted light. (Viewing by 
transmitted light was found to be most effective for detection of fine 
through-checks). The correlation between the numerical and subjective 
quality ratings is shown graphically in Figure 4. 

A comparison of the standard deviations with the mean quality 
ratings for all specimens is given in Figure 5, each rating being the 
average of all observations by the three observers. It is obvious that 
the scattering of results is essentially random in nature, although it 
would appear that in the numerical scheme, low quality sheets had 
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“Good” quality veneer sheet (Sheet No. 5) 
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FIGURE 2. A “Very bad” quality veneer sheet. (This sheet was not included in the present 
experiment.) 


Numerical Scheme Subjective 


Sev. Cl. Num. Wt. Degr. Rat. 


Buckling 1 1 1 Very slight 
Face Checking 4 1 4 Severe 
End Splitting {' 4 19 Severe 

1 3 
Through-Checking 5 4 20 Very severe 
Quality Rating — _— 44 Very bad 
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FIG.4. CORRELATION BETWEEN SUBJECTIVE AND NUMERICAL 
QUALITY RATINGS. 
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smaller average standard deviations than higher quality ones. This 
trend is contrary to what one might expect, since low quality sheets 
offer more scope for error due to the higher numerical weights they 
usually carry on end splitting and through-checking. It might indicate 
that observers were more certain and more nearly unanimous in their 
judging of a sheet with a large number of serious defects than one with 
a smaller number. It is, however, possible that had there been a number 
of very good to excellent sheets, they would have been judged with 
greater assurance (note position of the point for sheet 14 in Fig. 5a). 

Table 5 shows the mean absolute differences between the first and 
the second degrade ratings by each observer for all degrade types in 
both schemes, and also the mean absolute differences between the 
quality ratings. Even though the numerical degrade ratings for end 
splitting and through-checking are multiplied by a numerical weight 
factor of 2, 3, or 4, so that their differences are of a greater magnitude 
than those for buckling and face checking, it would appear that all 
observers had more difficulty in being consistent in their numerical 
assessment of end splitting and through-checking than in assessing 
buckling and face checking. The effect of experience seems to be 
evident in a comparison of the mean difference of the numerical quality 
ratings of observer A (3.3) with those of the other observers (4.2 and 
4.6); the mean difference is smallest for the former. This is, however, 
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TABLE 5 
Mean Absolute Differences between First and Second Degrade and 
Quality Rating Evaluation by Each Observer 


Observer A Observer B Observer C 
Degrade type 

Num. Sub. Num. Sub. Num. Sub. 
Buckling 0.4 0.3 0.4 0.2 0.7 0.5 
Face checking 0.4 0.6 I 0.5 0.6 0.4 
End splitting 1.8 0.4 3.6 0.4 2.8 1.0 
Through-checking 2.8 0.4 2.3 0.4 3.2 0.6 
Quality rating 3.3 0.7 4.2 0.4 4.6 0.7 


Num. = Numerical scheme; Sub. = Subjective scheme. 


not the case for the subjective ratings for which observer B obtained 
the greatest consistency. It seems rather surprising that the personal 
factor of experience should play a greater role in the numerical than in 
the subjective scheme, and this observation may suggest that the 
numerical evaluation is influenced by subjective factors. 

Further analysis to elicit the reason for the uncertainty in the 
numerical evaluation of end splitting and through-checking showed 
that all observers had marked difficulty in the determination of severity 
classes, even though these are nominally based on supposedly unam- 
biguous measurement of the total length of checks or splits. This is 
due to the fact that an observer who has to inspect a large number of 
sheets in a limited time cannot possibly accurately measure the total 
length of the checks or splits present, but has to rely on visual estima- 
tion to make his assessment. The same applies to the numerical weights, 
which, in practise, were also often allotted by subjective judgment 
rather than detailed quantitative evaluation of degrade sub-types. 

There were also considerable inconsistencies in observers’ judgments 
of whether degrade was localized or not. 


5. STATISTICAL ANALYSIS OF RESULTS 


Analyses of variance were carried out on the experimental results, 
giving an assessment of the magnitude of effects of the following factors:— 


(i) Sheets (S) 
(ii) Observers (0) 
(iii) Repetitions by an observer (I?) 
(iv) Sheet-observer interaction (SO) 
(v) Sheet-repetition interaction (Si) 
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A summary of the analyses for quality ratings in both schemes is 
given in Table 6. 

Of the variance values given in this table, the observer component 
in both schemes, and the repetition component in the subjective scheme, 
are non-significant, but all other components are significant at the 1 per 
cent level. 

It should be noted that the highly significant sheet-observer inter- 
action in both schemes is ‘entirely consistent with the non-significance 
of the observer component. 


5.1 Variance of Sheet Means (Quality Ratings) 


The variance of sheet means, V, , in the present experiment was 
found to be 


V,= 
The value of this expression for the numerical scheme is 
V,mum.) = 5.01 


6 


so that the standard error of sheet means is 1.9. In the case of the 
subjective scheme, equation (1) has the value 


0.7928 
6 


and the standard error of sheet means is then 0.36. 

To predict the variance and standard error for any method of exam- 
ination, we may write, by analogy, the variance of sheet means if m 
observers each take n readings on each sheet as 


(SR) + n(SO) 


mn 


V.(sub.) = = 0.1321 


(2) 


which becomes, for the numerical scheme 


_ 10.45 , 5.61 
m 


V,(@mum.) + 


and for the subjective scheme 


n m 


Although equation (2) can be used to predict the performance of a 
single observer for purposes of comparison with a greater number of 
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observers, the analysis of an actual test by a single observer would have 
to be different from that given. It is clear that if there is only one obser- 
ver, no components of variance can be estimated for observers and 
sheet-observer interaction. 

The values of standard errors of sheet means (in per cent of the total 
range of ratings) for different numbers of observers and readings are 
plotted in Figure 6. It is obviously preferable to increase the number 


NUMERICAL SCHEME SUBJECTIVE SCHEME 


AL. 


STANDARD ERROR IN PERCENT OF TOTAL RANGE OF RATINGS 


2 4 6 610 20 40 60 100 2 4 6 8610 20 40 60 100 
NUMBER OF READINGS PER OBSERVER ’ NUMBER OF READINGS PER OBSERVER 


FIG.6 STANDARD ERRORS OF SINGLE SHEET MEANS OF QUALITY 
RATINGS 


of observers rather than to increase the number of readings per observer. 
Even the extreme case of one observer taking 100 readings will lead to 
a greater standard error than two observers each taking two readings, or 
three observers each taking one reading. 

The accuracy of the above formulae may be judged from the results 
of the drying study mentioned above (Kauman, Gottstein and Lantican, 
1956) in which observer A twice subjectively evaluated the quality of 
340 sheets from which the 20 specimens of the present investigation 
were later selected. The SR component of variance was 0.4225, com- 
pared with 0.3884 in the present results (Table 6). Using 0.2022 as 
the value of the SO interaction (no value is available for the one observer 
in the drying experiment), the standard error of sheet means in the 
drying experiment was 0.643, a difference of only 2 per cent from the 
predicted value (0.630). 
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5.2 Variance of Observer Means (Quality Ratings) 


Although the variance of sheet means is of greatest importance 
when defining desirable levels of accuracy and assessing the usefulness 
of an evaluation scheme, it is also of some interest to consider the vari- 
ance of observer means. The latter measures the internal consistency 
of results for any observer and is appropriate for testing the significance 
of the differences between observers. It could be useful in deciding 
whether observers can be changed in the course of an experiment. 

In the present experiment, the variance of observer means, V, , is 
given by 


_ (SR) + 2(SO) + 20(R) 
40 (3) 


Upon substitution of the components of variance, equation (3) 
yields for the numerical scheme 


Vo 


V,(num.) = wees = 1.68 (Standard error 1.3) 
and for the subjective scheme 
V,(sub.) = os = 0.01982 (Standard error 0.14) 


The general expression for the variance of observer means, for 
observers taking n readings on each of p sheets, is 


(SR) + n(SO) + p(R) 


np 


(4) 


the formula for the numerical scheme being 


0.45 2.28 
1 45 5.61 
np n 


V,(num.) = 


and for the subjective scheme 


0.3884 , 0.2022 
V,(sub.) = + 


As stated earlier, the observer means in the present experiment did 
not differ significantly. 

It should be noted that observer variance had to be calculated from 
three mean squares in the analysis of variance, unlike sheet variance 
which could be found from one mean square. Observer variances are 
thus not very accurate. 
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Figure 7 gives the standard errors of observer means (in per cent 
of the total range of ratings) for different numbers of sheets and readings. 
As would be expected, the variance of observer means decreases with 
increase in the number of readings. An appreciable improvement in 


NUMERICAL SCHEME SUBJECTIVE SCHEME 


STANOARD ERROR IN PERCENT OF TOTAL RANGE OF RATINGS 


NUMBER OF SAMPLES EXAMINED NUMBER OF SAMPLES EXAMINED 40 60 100 


FIG.7 STANDARD ERRORS OF OBSERVER MEANS OF QUALITY RATINGS 


observer consistency may be obtained by increasing the number of 
readings from 1 to 2, but an increase in the number of readings beyond 2 
does not result in improvements of any great importance in either 
scheme, except for a very small number of sheets. 

The actual observer means of the quality ratings in both schemes 
may be found in Table 7. 


5.3 Results for Degrade Ratings 


The degrade ratings for the individual degrade types in both schemes 
were analyzed in the same manner as the quality ratings. Means and 
standard errors for each degrade type (and also for quality ratings) are 
given in Table 7. In the numerical scheme, observer differences were 
highly significant only for face checking. Significance at the 5 per cent 
level was found for buckling, but differences for end splitting and 
through-checking were not significant. In the subjective scheme, 
observer differences for degrade ratings were found highly significant 
for buckling, face checking and through-checking, but not for end 
splitting. The differences were such that observer differences for 
quality ratings were not significant, as noted earlier. 
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The standard errors are a measure of the variation between observers, 
but also of the internal consistency of the observations by each observer. 
Thus, for instance, although the observer means for end splitting and 
through-checking in the subjective scheme are almost identical, the 
standard errors differ considerably, indicating that the internal consist- 
ency of observers in the evaluation of through-checking was considerably 
greater than for end splitting. No such differences in the standard 
errors were observed in the numerical scheme. 

The high significance of observer differences found for buckling and 
face checking might be expected since there is usually less variation 
within these two degrade types than in end splitting and through- 
checking. In other words, the consistency within observers, and hence 
the significance of observer differences should be higher for the former 
than for the latter. With the exception of the subjective degrade rating 
for through-checking, this was indeed the case. An alternative explana- 
tion could be that there was less variation in through-checking than in 
end splitting, but this seems hardly borne out by the results of the 
numerical scheme. 

It is interesting to note the correlation between observer means in 
the two schemes. Observer C, for instance, was consistently high in his 
evaluation of buckling, but low for end splitting. B obtained high 
readings for end splitting and face checking, the latter particularly in 
the numerical scheme where he appears to have classified a proportion 
of through-checks as face checks. It might be noted here that differentia- 
tion between face checks and very narrow through-checks was reported 
particularly difficult by all observers. 


6. DISCUSSION 


The present experiment has shown that subjective evaluation can 
yield results of an accuracy approaching that of the numerical scheme, 
although the accuracy of the latter was slightly superior. 

If it is postulated that the maximum permissible standard error of 
single sheet means for quality ratings should not exceed 5 per cent of 
the total range of ratings, then the standard error should be less than 
2.5 in the numerical scheme, and less than 0.4 in the subjective scheme. 

Inspection of Figure 6 shows that this condition is fulfilled in the 
numerical scheme for three observers taking one reading, or for two 
observers taking two readings each. One single observer would have 
to take about 100 readings to obtain a similar accuracy, although 17 
readings will give a standard error below 5 per cent. 

In the subjective scheme, on the other hand, two observers would 
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have to take + readings, and three observers 2 readings each to obtain 
the desired accuracy. <A single observer could never realize the postu- 
lated accuracy even with an infinite number of readings. 

More specifically, a comparison of variance ratios shows that in 
general, the subjective scheme requires 40 per cent more readings than 
the numerical scheme to give a similar accuracy. 

However, numerical quality evaluation, necessitating detailed rating 
of all degrade types, occupies considerably more time than subjective 
quality rating. The approximate average times were as follows: 


(i) Numerical evaluation of quality rating (which in- 


cludes detailed evaluation of degrade ratings) 3 min 
(ii) Subjective evaluation of quality rating only 0.5 min 
(iii) Subjective evaluation of quality rating and degrade 

ratings 1.5 min 


Figure 8 gives a graphical comparison of the time required to realize 
a desired degree of accuracy by different numbers of observers and 
readings when using the various evaluation schemes. A logarithmic 
time scale has been used for convenience. The plots indicate that in 
applications where only quality ratings are desired, the subjective 
method, though slightly less accurate for a given number of observers 
and readings, is clearly preferable on account of the much greater 
rapidity of evaluation. 

If a detailed record of degrade ratings is desired, the difference in 
time required by the two schemes for quality evaluation by up to three 
observers is largely eliminated. If four or more observers take part in 
the evaluation, the subjective scheme has again the advantage of re- 
quiring considerably less time. 

Where it is inconvenient to use more than two observers, the numer- 
ical scheme has a margin over the subjective scheme as regards accuracy 
and time required. In addition, the numerical scheme records greater 
detail and would therefore permit finer analysis of the factors contribut- 
ing to degrade. 

Neither of the schemes appeared to have any advantage over the 
other regarding the relative difficulty and the time required for the 
training of observers by personal tuition, but it is considered that the 
numerical scheme would be more suitable for initial training by written 
instructions only, and for standardization purposes. 

As noted before, the quality evaluation schemes presented above 
were developed to deal with the drying of ‘‘ash” eucalypt veneers which 
are prone to considerable drying degrade and collapse. The schemes 
are, however, sufficiently flexible to be applied to other types of veneer 
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FIG. 8 TIME REQUIRED FOR QUALITY EVALUATION 


and possibly to sawn timber; criteria for evaluation can easily be re- 
defined to take into account different forms of degrade. 

It will, in general, be simpler first to re-define the numerical scheme 
for a given class of material, and to use the criteria thus obtained for the 
re-definition of the subjective scheme. In the former, it will usually be 
possible to place different degrade types in proper perspective by re- 
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allocation of appropriate numerical weights. It will be convenient to 
allocate the weights in such a way that the highest possible numerical 
quality rating is a round number (e.g. 50 or 100), although this is not 
absolutely essential. 

In some cases, it might appear desirable to change the criteria for 
severity classes. It would, for instance, be possible to use the total 
surface area instead of the total length of checks, but such a procedure 
would entail considerable difficulties of measurement. In cases where 
buckling is the major form of degrade, it might be deemed convenient 
to define more closely the criteria for severity ratings for this degrade 
type. 

Although observers found no apparent difficulty in dissociating 
defects present prior to drying from the actual drying degrade, it might 
be of value to evaluate quality ratings before and after drying, taking 
into account all defects present, and subtracting the former rating from 
the latter to obtain “dried quality’. 

It might also be proposed to change the rating of a degrade type 
abruptly above a suitably defined threshold level. In the authors’ 
opinion, however, the evaluation schemes presented here are not suitable 
for defining a level of acceptance or rejection of a sheet in this manner. 
Their function is rather to provide continuous and reasonably uniform 
scales for the evaluation of degrade types in order to compare the effect 
of different experimental conditions. 

Having re-defined the numerical scheme to deal with a given set of 
requirements, it is an easy matter to apply the new definitions to the 
subjective scheme. In the latter, it might under certain circumstances 
be of value to increase the number of quality ratings; by such a procedure 
it might be possible to decrease the standard error. 

Although the association of subjective ratings with numbers greatly 
facilitates subsequent analysis of the results, it is considered that the 
verbal designations of quality ratings (excellent, very good, etc.) and 
of severity classes (very slight, slight, etc.) should be used by observers 
until sufficient experience in using the scheme has been gained. The 
use of intermediate ratings, such as slight to moderate, or fair to poor, 
though a convenient way of escape from making a decision, should be 
discouraged. If observers find the number of ratings available in- 
sufficient, it would be better to increase their number, as suggested 
above, so that each rating can be associated with an integral index 
number. 

It may be considered desirable, under certain circumstances, to 
define quality rating so that the highest rating corresponds to the best 
quality. This can easily be achieved by re-defining the quality rating 
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as follows: 


Quality rating = 50 — (Degrade ratings). 
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THE ESTIMATION OF AGE-SPECIFIC INFECTION RATES 
FROM A CURVE OF RELATIVE INFECTION 


P. WuittLe 


Applied Mathematics Laboratory, New Zealand D.S.I.R. 
Summary 


Consider a sample of individuals (in which the age classes are not 
necessarily proportionally represented) from a population subject to 
some permanent infection. A procedure is described for the estimation 
of the infection rate as a function of age, and of the expectation of life 
after infection. The method is applied to unpublished data supplied 
by P. C. Bull on the incidence of liver lesions due to coccidiosis (Eimeria 
stiedae) in the New Zealand wild rabbit. 


(1) The model 


Disease and parasitism belong to those features of an animal popu- 
lation for which laboratory examination must be supplemented by 
sampling investigations in the field if a satisfactory picture of the 
natural state of the animal is to be obtained. The sampling approach 
is regrettably indirect, but can prove informative. 

Consider an animal population subject to some permanent infection 
with recognisable symptoms, and suppose the population divisible into 
suitable age categories, e.g. 0-1 weeks, 1-2 weeks, etc. Suppose now 
that a sample is taken from the population in such a way that for any 
age category healthy and infected animals will be represented in their 
correct proportions (apart from sampling fluctuations). The sample 
will presumably include animals of a range of ages, and we shall suppose 
that in the jth age category h; healthy and 7; infected animals are 
captured. 

If one could be sure that all age categories were proportionally 
represented, and if the expected mortality due to natural causes at all 
ages were known, then it would be a direct matter to estimate the rates 
at which healthy animals were being infected and infected animals 
were dying of disease. Failing this the relative proportions of infected 
animals in the different age groups 


P= h; + i; (1) 
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prove to be useful quantities. We shall call the graph of these pro- 
portions against age the curve of relative infection. 

We shall assume that an animal healthy at age ¢ has a probability 
d(t) dt of becoming infected in the age interval (t, ¢ + dt) so that A(é) 
is the infection rate for age ¢. We shall suppose that there is a similar 
natural death rate a(t) which applies equally to healthy and infected 
animals of equal ages. Thus, in the case of the rabbits which we shall 
consider later, there are the natural hazards to young rabbits of hawks, 
etc. Finally, we shall assume that animals which contract infection at 
age t-s have a death rate due to the infection of 8(s, t) (so that B(s, t) dt 
is the probability of a life between ¢ and ¢ + dé, or a sickness duration 
s and s + dt, given that infection was contracted at age ¢-s). 

A useful extension of the model would be to allow recovery after 
infection, but to attempt estimation of recovery rates as well as infec- 
tion and death rates is to ask too much of the data. It is only thanks 
to rather special circumstances and by introduction of specialising 
assumptions that we can extract as much information from the data as 
we do. Further, many animal infections are permanent in that once 
they reach a certain stage they can only end in death—the animal 
either dies of sickness or is so weakened as to fall an easy prey to its 
enemies. 

If we begin with a group of N animals born at time t = 0 (so that 
i can be used to represent age or time indiscriminately) then it follows 
by conventional reasoning (see, for example, Feller p. 364) that the 
expected number of healthy animals remaining at time ¢ is 


H(t) = N exp [ [a(u) + in} (2) 
The expected number becoming infected in the age interval (¢, ¢ + dt) is 
K(t) dt = H(t)X(6) dt (3) 


whence the total number of infected animals living at time ¢ is 


I(t) = [ K(u) exp ‘-f [a(v) + BW — u, v)] aw} du 


= N exp an} [ exp NOW (4) 


) ~ [ BY — u,v) du 


The ratio of the expected numbers of infected and healthy animals of 
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age ¢ is thus 


R(t) = = [ exp if [A(v) — Bw — u, du (5) 


This ratio is independent of the natural death rate a(t), as might have 
been expected. 

The theoretical curve of relative infection, with which the empirical 
curve (1) is to be compared, is 


PO = TH) + ~ T+ RO (6) 


Now, we cannot estimate two arbitrary functions \ and 8 on the 
basis of a single empirical function p;. Only by restricting one or both 
of \ and 8 can we hope to obtain determinate estimates. In the next 
two sections we shall consider two such specialisations: the extreme 
cases of a fixed sickness duration and of a random sickness duration. 

Let it be noted, that we are not completely without information on 
the functions \ and 6. By their nature both are positive and are unlikely 
to have more than one peak. This knowledge may be turned to good 
account. 


(2) The case of fixed sickness duration 


Suppose that an animal contracting infection at age ¢ dies after a 
further fixed period of time L (it is easy to generalise to the case when 
the sickness duration L is a function of age at the time of contraction). 
Expression (4) then becomes 


= Nexp ole) au\| exp Au) au) (7) 
— exp {- d(u) au 


where c(é) is the earliest age at which a living infected animal of age ¢ 
could have contracted infection, and is thus given by 


c(t) = 0 (t < L) 


(8) 
ct) = t-—L (t > L) 
The proportion of infection is then a simple function of 2: 
P(t) = 1 — exp {- A(u) au) (9) 
e(t) 
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Putting 
Q(t) = —log [lL — PO] (10) 
we can invert (9) to obtain an expression for \: 


At) =  (t2 L) 


In practice, if p; is known at (say) weekly intervals then equations (11) 
will be applied in difference form to provide a solution for \. Thus, if 


qi = —log (1 — p,) (12) 


and the expected sickness duration is an integral number of weeks, L, 
then solution of the equations 


i, = — (j < L) 


(13) 
i, = (j = L) 


will yield a set of constants 4, where i; can be approximately interpreted 
as the estimated probability that an animal healthy when 7 weeks old 
will contract infection in the coming week. (Rather, that an animal 
apparently healthy when j weeks old will show symptoms in the coming 
week, since infection will presumably not manifest itself immediately). 

If the expected sickness duration LZ is not known, then we must 
choose a value for which the solutions 4; of (13) are positive, and, pre- 
ferably, unimodal. In this way, one can fix an upper bound for the 
true L, since substitution of too large an Z in (13) will produce an 
oscillating solution. 


(3) The case of random sickness duration 


We shall go to the opposite extreme in this section, and suppose the 
sickness death rate to be independent of sickness duration, so that we 
can write 


Bs, ) = (14) 
It is now readily verified that 


R(t) = i { — BO) (15) 


is a solution of the differential equation 


R(t) + [BO — AMIR] = (16) 
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so that the solution for \(¢) in terms of observed quantities becomes in 
this case 


RO) + BOR) 
1+ RO) (17) 


In the previous section a knowledge of L (or L(t)) was necessary to 
enable us to complete our solution: in this case we require correspondingly 
a knowledge of 8(¢), the age-specific death intensity. We can, of course, 
suppose it equal to a constant and thus make possible an approximate 
solution for \. Just as the positivity of \ enabled us before to find an 


upper bound for L, it enables us now to find a lower bound for 6. We 
have from (17) 


A(t) = 


p> (18) 


In this way we obtain a lower bound for 8, but it seems that one cannot 
establish an upper bound (or a lower bound for the fixed sickness dura- 
tion LZ of the previous section) without making further restricting 
assumptions. We shall return to this question in the next section. 

We have framed our model in continuous time, as is natural. How- 
ever, since observations are taken at discrete intervals of time the 
differential relation (16) must be approximated by a difference equation. 
Thus, if we have observed values r; of the ratio R(t) taken at unit 
intervals of time, then the obvious adaption of equation (17) would be 


= — 73) + +173) 


This is the estimating equation we shall use in the next section. 


(18) 


(4) An application to observed data 


The data came from an extensive investigation by Bull (now being 
prepared for publication) on the incidence of parasites in the New 
Zealand wild rabbit. The particular parasite considered is Eimeria 
stiedae, (Protozoa), which destroys the epithelial cells of the bile ducts 
and leads to the appearance of severe lesions in the liver (Becker, 1934). 

The sample was one of 4143 rabbits, whose ages were estimated 
from the paunched weight. (The weight/age curve for diseased animals 
may well differ from that for healthy animals. However, this point has 
not been investigated and in this present analysis we can make no 
allowance for it.) Animals were counted as infected if more than half 
the liver showed severe lesions. Stephens (1952) referring to 1’. stiedae 
in young rabbits states that “it can be assumed that in severe attacks 
the animal generally dies before reaching maturity.” 
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TABLE 1 


Incidence of Parasites, New Zealand Wild Rabbit 
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Age Estimated infection rate 
in p r ¥ - 
wecks B = 0.45 B = 1.00 

3 0.090 

0.093 —0.012 0.027 0.074 
+ 0.080 

0.169 0.163 0.204 0.284 
5 0.200 

0.689 0.878 0.703 0.928 
6 0.530 

1.473 0.689 0.547 0.874 
7 0.645 

2.163 0.692 0.526 0.903 
8 0.715 

2.716 0.413 0.440 0.842 
9 0.745 

2.747 —0.351 0.236 0.639 
10 0.720 

2.175 —0.793 0.059 0.435 
il 0.640 

1.453 —0.650 0.002 0.327 
12 0.530 

0.919 —0.319 0.049 0.313 
13 0.415 

0. 6g2 —0.134 0.094 0.309 
14 0.365 

0.518 —0.115 0.078 0.265 
15 0.315 

0.425 —0.071 0.084 0.248 
16 0.280 

0.361 —0.056 0.078 0.224 
iy 0.250 

0.308 —0.051 0.067 0.196 
18 0.220 

0.259 —0.047 0.055 0.168 
19 0.190 

0.228 —0.015 0.072 0.173 
20 0.180 

0.213 —0.015 0.067 0.163 
21 0.170 

0.201 —0.007 0.069 0.162 
22 0.165 

0.194 —0.008 0.066 0.156 
23 0.160 | 

| 0.173 —0.034 0.038 0.118 


= 
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TABLE L- Concluded 


Age Istimated infection rate 
weeks B= 0.45 B = 1.00 

24 0.135 

0.143 —0.026 0.033 0.102 
25 0.115 ; 

0.124 —0.013 0.038 0.099 
26 0.105 

0.108 —0.018 0.028 0.081 
27 0.090 

0.093 —0.012 0.027 0.074 
28 0.080 

0.076 —0.023 0.010 0.049 
29 0.060 


A graphical interpolation by weeks of the age-grouped data yielded 
the curve of relative infection in the diagram (full line), the figures p; 
are given in the second column of Table 1. 

An application of the method of section 2, assuming a constant 
sickness duration, was not at all successful. One had to choose an 
absurdly low value of L in order to obtain a \ solution which would not 
oscillate violently. The reason for this becomes apparent on an inspec- 
tion of the diagram. The curve rises much more quickly than it falls. 
If sickness duration were constant, however, the large increase of infec- 
tion in the 5-6 week period would be followed by a roughly corresponding 
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decrease some time later when this age-group died off. The fact that 
there is no such decrease indicates that the hypothesis of a fixed sickness 
duration is untenable. In fact, the gradual decay of the curve indicates 
that sickness duration must be fairly broadly distributed. 

Let us consider the alternative hypothesis of a uniform death rate. 
In the third and fourth columns of the table we have recorded the 
values of 3(r;,, + 7;) and 3(r;,, — r;) which we shall take\as estimates 
of R(t) and R’(t) at the mid-week ¢ = j7 + 34. The maximum observed 
value of —R’/R occurs at 113 weeks: 0.447. The constant 8 must thus 
be greater than 0.447, but probably not much greater, since A(t) prob- 
ably assumes fairly small values once the animal has passed a critical age. 

For trial values of 8 we shall take 8 = 0.45 and 1.00, corresponding 
respectively to expected sickness durations of 1/0.45 = 2.2 weeks = 
15.4 days, and 1/1.00 = 1 week. (These figures, particularly the first, 
are quite consistent with laboratory investigations of EZ. stiedae infection, 
which give a sickness duration of between 21 and 30 days (Becker, 1934). 
Since lesions appear in the liver about 10 days after infection, the actual 
time from appearance of the symptoms until death is of the order 11-20 
days. However, in Bull’s sample, rabbits were only counted as infected 
if more than half the liver showed severe lesions, and expectation of 
life for animals in which infection is so far advanced would almost 
certainly be less than 11-20 days.) 

The values of \ estimated using formula (18) with 8 = 0.45 and 1.00 
are recorded in the table and on the diagram. Neither of the two 
estimated curves conforms completely to our ideal of a smooth unimodal 
curve, but there are at least three reasons for this: 


(a) The model includes several rather drastic assumptions (in particular 
the assumption that the death rate A(s, ¢) is constant can be valid 
only as a very rough approximation) and the fitting of an over-simplified 
model will produce systematic errors in the fitted curve. 


(b) The curve is subject to sampling fluctuations. These have been 
largely smoothed out in the present case, but may have some effect at 
the extremes of the curve, since there were few rabbits in the extreme 
age categories. 


(c) In smoothing a curve by eye one can unknowingly introduce spurious 
deviations, which will become more pronounced if the curve is differ- 
entiated. Further, if one substitutes difference quotients for differential 
coefficients, one introduces a bias which is negligible if the gradient is 
small, but which may be appreciable when the ordinate changes quickly 
as it does at ¢ = 5-6 weeks. 
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The curve for 8 = 1 is more regular than that for 6 = 0.45, but at 
least part of this gain in regularity is illusory. If we let 8 become large 
in formula (17), we find that \(é) becomes proportional to P(é). The 
same holds if we let I become small in (11). In other words, the larger 
8 or the smaller L, the nearer will the estimated i, curve approach to 
the original p; curve. In this case, as in most based on a large sample, 
the p, curve is itself regular (i.e. smooth and unimodal), so that as we 
choose a larger 8 or a smaller L we shall observe an increase in regularity 
of the estimated i; curve which is quite deceptive. 

One would naturally wish to have an objective method of reaching 
determinate estimates of the rate \ and of testing the final fit of the 
model. However, this will be possible only if one is justified in making 
a number of very restrictive assumptions (concerning both the functional 
forms of a, 8 and \, and the nature of the sampling fluctuations) or if 
one can take a sample which is inherently more informative (i.e. sample 
age groups proportionally). 

However, it may be asserted that the method as it stands gives one 
a positive indication of the upper limit of sickness duration; and a 
lower limit can often be set from one’s experience of the particular 
infection. In the case of the rabbit data considered in this section it 
was even possible to form some conclusions on the nature of the sickness 
duration (i.e. that it was broadly rather than sharply distributed). 
Finally if the distribution of sickness duration is known, one can evaluate 
the age/susceptibility function \(¢); while if little or no a priori informa- 
tion on sickness duration is available, one can still set plausible limits 
on the form of A(#). Thus, in the present case, the true age/suscepti- 
bility curve for rabbits almost certainly lies between the two estimated 
curves. 


The author’s thanks are due to Mr. P. C. Bull for valuable advice. 
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AN EVALUATION OF THE REMOVAL METHOD 
OF ESTIMATING ANIMAL POPULATIONS* 


CALVIN ZIPPIN 


Cancer Research Institute, School of Medicine, 
University of California, San Francisco, California 


I. INTRODUCTION 


Moran (1951) presented a method for obtaining maximum likelihood 
estimates of population size from the results of a series of trappings in 
which the trapped animals are removed from the population. This type 
of trapping program may be appropriate when auxiliary studies of the 
trapped animals are to be made which necessitate their eventual sacrifice, 
or when economic or health reasons make it inadvisable to return trapped 
animals to the population. In such situations it is not possible to use 
the tagging and recapture method which has been discussed by numerous 
authors in recent years (see e.g. Leslie and Chitty (1951), Bailey (1951)). 
Moran’s method, which will be called the removal method, is a special 
case of a more general procedure described by De Lury (1947). 

The present paper will cover the following topics: 

1. Development of a rapid graphical procedure for obtaining maxi- 
mum likelihood estimates of population size from removal method data. 

2. Determination of the asymptotic precision of the removal method. 

3. Report on the results of experimental sampling which was con- 
ducted in order to compare the maximum likelihood estimates with 
estimates by an alternative regression method (Hayne, 1949) and to 
study the performance of the removal method on small populations. 

4. Determination of the proportion of the total population which 
must be trapped in order to reduce the coefficient of variation of estimates 
of population size to specified levels. 

5. Tests of the assumptions underlying the removal we and their 
application to actual trapping data. 


*From the Department of Biostatistics, The Johns Hopkins University School of Hygiene and 
Public Health, Baltimore, Maryland. Department Paper No. 311. This work was done under Navy 
Contract Nonr-248 (16), 
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Il. THEORY OF THE REMOVAL METHOD 


The removal method assumes a stationary population during the 
trapping program and also that the probability of capture during a 
given trapping is the same for each animal and does not change from 
trapping to trapping. 

In this section the theory leading to maximum likelihood estimates of 
N, the population size, and p, the binomial probability of capture during 
a single trapping, will be reviewed for the purpose of introducing a rapid 
graphical method of obtaining estimates. New formulas for the standard 
errors of these estimates will also be presented. 


Conditional Binomial Approach 


Under the assumptions presented above, Moran applies the method of 
maximum likelihood to obtain conditional binomial estimates of the 
population size and of the probability of capture during a single trapping: 

Let 


N = population size, 
probability of capture during a single trapping, 
= number of animals captured during the 7th trapping, 


and 
zx; = >,)21 y; = total number captured prior to the ith trapping. 


The probability of capturing y; animals during the ith trapping, 
given that x; animals had previously been captured is therefore 


where g = 1 — p = probability of escaping capture during a single 
trapping. 

The joint probability or likelihood of the sample of catches actually 
observed in k trappings is 


N! 
yt? 


P(S) = I] = pig (2) 


where 


k k 
T= and Q= Dy. 
i=1 i=1 
This same result may be arrived at by looking at the problem slightly 
differently. _We may consider the entire population as falling into 
(k + 1) categories—those captured during each of the k trapping 
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periods and those escaping capture. ‘The probability that an individual 
is captured during the first trapping is p. Under the assumption that 
the probability of capture does not change from trapping to trapping, 
the probability of escaping capture during the first trapping period and 
being trapped during the second trapping interval is pg. In general the 
probability of being captured during the ith trapping interval is pq‘ 
and the probability of escaping capture during each of the k trappings is 
q'. According to the multinomial probability distribution, the prob- 
ability that y,; , y2, ... , y, animals are trapped during the k successive 
trappings and that (N — 7’) animals escape capture is given by 
Yilyo! --- — 7)! 


which is another form of equation (2). 
The logarithm of the likelihood of the sample of catches is then 
L = log P(S) = T log p + (kN — Q) log q 


(3) 
+ log N! — log (N — T)! — log 


i=1 


Maximizing L with respect to p, Moran obtains as the estimator 
of the probability of capture an expression which is equivalent to 


(4) 


For his estimate of N, Moran obtains the value N which satisfies 
the equation 


Moran suggests solving equations (4) and (5) by a process of suc- 
cessive approximations. However, since @ = 1 — 7, from equation (4), 


*In this equation log — T)/N is used to approximate 1/N + + -T 
+1). It has been pointed out by one of the referees that log (® -T + 1) /(N + }) is a better ap- 
proximation. 
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Substituting this expression for @ in equation (5) we obtain as the 
estimate of V that value V which satisfies the equation 


(6) 


. This equation allows an estimate of N without first approximating 
p, although the process of substituting values of N in both sides of 
equation (6) until it is balanced is frequently quite tedious. 

It is of interest to consider a partition of the joint probability of 
the sample P(S). For any specified animal, the probability that it 
escapes capture during all k trapping periods is g*. Since the chance 
of being captured is assumed to be independent from animal to animal, 
the total number, 7’, of animals captured follows the binomial (p;+ qr)”, 
where gr = q' and pr = 1 — qr. 

Hence the partition 


P(S) = P(S | T)P(T) 


is found from equation (2) to be 


(7) 


It may be noted that N appears only in the second factor P(T). This 
shows that the maximum likelihood estimate of N is obtained from the 
binomial distribution of the total catch 7’. If N is estimated by equating 
T to its expectation N(1 — gq‘), it may be shown that this estimate will 
differ from the correct integral maximum likelihood estimate by less 
than unity. Hence, as pointed out by Moran, we may estimate N by 


N= (8) 


which is an alternative form of equation (5). 

Further g and p are estimated from the first factor P(S | 7), the 
probability of the distribution of catches. For if we take the derivative of 
log P(S) with respect to p, the contribution from the second bracket in 
(7) is 

MN — 7) 


(1 — q 


It is easily verified that this contribution vanishes when V = 7/(1 — @'). 
These results form the basis for a quick graphical method of estimating 
p and hence N. 


N kN — x, 
j 
| 
E 


Va 


ANIMAL POPULATIONS 167 


Setting the derivative of log PCS | T) with respect to p equal to zero, 
we find 


(9) 


k 


_ kT —Q 


i=1 
T 
R may be calculated quickly from the data of a set of trappings. 
The value f which satisfies (9) is the maximum likelihood estimate of p, 


and @ = 1 — jp. Substitution of T and ¢‘ in equation (8) provides the 
maximum likelihood estimate of N. 


where R 


Graphical Method of Estimating p and N 


Figure | has been drawn up to provide a means for rapidly obtaining 
p corresponding to a value of R calculated from the data of 3, 4, 5, or 7 
trappings. Estimates of (1 — g*) may be obtained directly from Figure 2. 
Similar graphs may be constructed for values of k other than those given 
here. 


With the aid of these graphs, the procedure for obtaining an estimate 
of N is as follows: 


1. Calculate 


2. Find estimate of (1 — q') corresponding to R from appropriate 
graph in Figure 2. 
3. Calculate 


T 


N= 


Example: 
Assume the following catches in a three night trapping program: 


y, = 165, y. = 101, ys = 54. 


Then, 
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Precision of Maximum Likelihood Estimates 

Assuming, as does Moran in the example which he gives in his paper, 
that the asymptotic properties of maximum likelihood estimates hold in 
the cases of N and #, the asymptotic variances of these estimates are: 
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where .V/, the variance-covariance matrix of Vand p is 


aL 


aN? log N! — an? log (N — T)!, 


Since 
E(T) = NU 
E(N — T) = Nq‘. 


It may also be shown that 


Since E(N — T) = Ng‘, then asymptotically 


_ log N! log Nq'! (10) 


aN* aN* 


E{a°L/aN*} may be evaluated exactly from tables of the trigamma 
function or may be approximated using Stirling’s approximation for 
factorials. If Stirling’s second order approximation is used, 


With N large, i.e. 200 or more, and (1 — gq‘) equal to .9 or less, the 
last two terms within the brackets become negligible. Therefore, the 
large sample approximation is 


’ OL | 

| op ON op q 

| In the above matrix 

oN? 

and 
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Substituting the expected values for T and Q and reducing, we obtain 


= NO (13) 
P14 


k 


Thus the asymptotic variance of N is 


V(N) = 


If the approximation for F’ given in (12) is used, 


(1 — — 


Figure 3 gives graphs of the standard deviations of estimates of N, 
as calculated from equation (15), for populations of size 50, 100, 300 
and 500 subjected to specified probabilities of capture for 3, 5, and 7 
trappings. Values of F’ were taken from tables of the trigamma func: ion 
(Davis, 1935). For a given N and k (number of trappings), it may be 
seen that, in general, as p increases the calculated standard deviation 
of the estimate of N decreases. For N = 50 and N = 100, the standard 
deviation appears to increase slightly as p becomes very large. Also at 
very low levels of p the trend of the standard deviation is abnormal. At 
high or low levels of p, the expected catch during the last trapping, 
Npq‘", is small and the asymptotic theory may not apply. 


V(N) = 


(16) 


Calculation of Standard Errors from Observed Data 


The variance of an estimate of N obtained from a given set of data 
is computed by substituting sample estimates for the parameters in 
equation (15). A simple form for calculating the variance of N is 


qT 


(17) 
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ASYMPTOTIC STANDARD DEVIATIONS OF MAXIMUM LIKELIHOOD 


ESTIMATES OF N=50,/00,300, 500 
(3,5,7 TRAPPINGS) 


N= 500 


NO. TRAPPINGS 


STANDARD DEVIATION OF ESTIMATES OF N 


PROBABILITY OF CAPTURE DURING A SINGLE TRAPPING 


FIGURE 3 


where 


k = number of trappings, 

T = total catch, 

_ Flog N! log 
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and p and @ are the maximum likelihood estimates of p and g. The 
variance of an estimate of p may be determined by 


Example: 


Referring to the previous example, T = 320, k = 3, N = 400. 

From Figure 1A, R of .65 corresponds to p = .42; hence @ = .58. 

F may be obtained from tables of the trigamma function (Davis, 
1935) by evaluating 


p — 2 log N! log (NG)! 


or, using the large sample approximation of equation (12) 


Applying either of the above procedures we find 


F = —.010. 
Hence, applying (17), 
T 185.60 
Vim) = = 691. 
™) = (er — 2084 
and 
SE(N) = V V(N) = 26.3. 
Also 
_ 
V(p) = — (kp? 2084 -002209 
and 


SE(p) = = .047. 


The 95% confidence limits for N and p are approximately 347 to 
453 and .33 to .51 respectively, i.e. 400 + 2(26.3) and .42 + 2(.047). 
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Regression Method of Estimating Population Size from Removal Data 


Hayne (1949) described a regression method of estimating population 
size which is based on the same assumptions as the maximum likelihood 
method. 

If p is the probability of capture during a single period of trapping, 
the number expected to be captured in the first period is Np; the number 
expected to be captured in the second period is (NV — y,)p, where y;, is the 
number caught in the first period. Similarly, the number expected to be 
caught in the ith period is 


= (N — x,)p 


where N is the original population size and z; represents the number of 
animals caught prior to the ith period. This equation may be written 
in the form 


E(y;) = Np — pz; . (19) 


This is the equation of the regression of the number caught in the 7th 
period on the number caught prior to the 7th period. The absolute value 
of the slope of the line estimates the probability of capture and the x 
intercept estimates the population size. 

There are several ways in which the regression of y; on x; may be 
fitted. The simplest way, suggested by Hayne, is to fit a line by eye to 
the plotted data. This method is demonstrated in Figure 4, using the 
data of the previous example (successive catches of 165, 101, 54). Other 
ways are to calculate an unweighted least squares regression or a weighted 
least squares regression. Since the binomial variance of the catch in the 
ith period is (N — x;)pq, a weighted regression might assign weights to 
this catch inversely proportional to (NV — z;), using the intercept of a line 
drawn by eye as a first approximation to N. 

It is easy to verify that the regression and maximum likelihood 
methods give identical estimates for data collected on two trappings. 
For three or more trappings the estimates by the two methods need not 
be the same. The asymptotic variance of an estimate obtained by an 
unweighted least squares procedure, using data from three or more 
trappings, is not easily determined. Estimates from a weighted least 
squares regression may be seen to be related to estimates obtained by 
the method of minimum chi-square. Minimum chi-square estimates 
are those values of N and # which minimize the expression 
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ESTIMATION OF POPULATION SIZE 
BY REGRESSION METHOD 
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FIGURE 4 


Here the weight of each point is inversely proportional to its binomial 
variance. Since the method of minimum chi-square is asymptotically 
the same as the method of maximum likelihood, it appears that the 
variance of estimates obtained from a weighted least squares regression is 
asymptotically the same as the variance of the maximum likelihood 
estimates. 


Sampling Experiments 


Sampling experiments were carried out in order to determine how well 
the large sample theory holds at various population levels as well as to 
compare the maximum likelihood and regression methods of estimation. 

Population estimates were made from ‘“‘catches” which were draw. 
at random with the aid of a table of random numbers (Snedecor, 1946) 
and Tables of the Binomial Probability Distribution (1950). The 
method of sampling produced “catches” which varied binomially about 
the expected “catches”. ‘Populations’ of sizes 49, 98, 196, and 392 
were subjected to a probability of capture of .4 for three ‘trappings’. 

In the case of three trappings, it may be shown that the maximum 
likelihood method does not give a finite estimate of N when the third 
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observed catch equals or exceeds the first catch.* This happened in 
four of the 328 samples drawn from N = 49. Estimates of p and N 
were made for the 324 remaining samples by the maximum likelihood 
and unweighted least squares regression methods. The frequency 
distributions of the estimates are shown in Figure 5 and are seen to be 


MAXIMUM LIKELIHOOD ESTIMATES 
OF N=49, 98, 196, 392 


100; 50 
N=49 


--- REGRESSION 
ESTIMATES 30 

— MAXIMUM 
LIKELIHOOD 

ESTIMATES 


10 


70 10 615006600 95— 130 


FREQUENCY 
uo 


50, 30, 
N=196 N=392 
sol 20+ 
10+ 
10+ 
150 190 230 270, 330 365 400 435 470 
N 
FIGURE 5 


skewed to the right. Seven maximum likelihood and five regression 
estimates fell outside the range of the graph. The distributions of the 
logarithms and the reciprocals of the estimates also show considerable 
skewness. 

Table 1 gives the mean, standard deviation, and per cent of the 
estimates by each of the two methods which fell within N = 49 plus 
or minus one and two times the theoretical maximum likelihood 
standard deviation, 11.56. 

Both observed standard deviations far exceed the theoretical value of 
11.56. Sample 19, with successive catches of 15, 9, and 14, gave a 


*When yi = ys, p = 0 and the estimate of population size is infinite. If ys > y:, the ratioR = 
DG — 1)yi/T exceeds unity, which cannot be satisfied (equation (9)) by values of p between 0 and 1. 
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TABLE 1 
Comparison of Two Methods of Estimation 


Regression Maximum Normal 
Parameters estimates likelihood curve 
estimates 

Number of estimates 324 324 
Mean 56.68 55.59 
Observed standard deviation 31.95 27.56 
Per cent falling within 49 + 1 

theoretical S.D. 76.9 78.4 68.3 
Per cent falling within 49 + 2 

theoretical S.D4s 90.1 90.4 95.5 


regression estimate of 175 and a maximum likelihood estimate of 342 
and is in large part responsible for the difference between the standard 
deviations of the two sets of estimates. 

The last two lines of Table 1 show that, compared to normal theory, 
too large a percentage of the estimates by either method lies within one 
theoretical standard deviation of 49, and too large a percentage lies 
more than two standard deviations away. 

Figure 6, which gives the cumulative percentage of estimates within 
a given absolute deviation of 49, shows little difference between the two 
methods. Of the 324 samples, 172 of the maximum likelihood estimates 
were closer to the true value of 49 than the corresponding regression 
estimate, while the reverse was true for the remaining 152 samples. 

Additional sampling was done for the maximum likelihood method 
with p = .4,k = 3, and N = 98,196, and 392. For N = 98 samples 
were obtained by combining pairs of the 328 samples drawn from 
N = 49. The four samples which had been excluded previously were 
retained in this distribution. Although no case was observed in which 
the first catch did not exceed the third catch, Sample 72, comprised in 
part of one of the rejected samples from N = 49, was the poorest with 
catches of 23, 28, 22 and a maximum likelihood estimate of VY = 1108. 

Eighty-two of the samples from N = 196 were obtained by combining 
successive groups of four of the original samples from N = 49. Eighty 
additional samples from N = 196 and the entire 160 samples from N = 
392 were drawn by use of random normal deviates (Wold, 1948). 

Frequency distributions of maximum likelihood estimates of N = 
98, 196, and 392 are shown in Figure 5. Table 2 gives the parameters of 
these distributions. 
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CUMULATIVE PERCENTAGE OF 
ESTIMATES WITHIN SPECIFIED 
ABSOLUTE DEVIATION FROM N=49 
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FIGURE 6 


As expected, agreement between the observed and theoretical results 
improves as N increases. ‘The means of the maximum likelihood esti- 
mates slightly overestimate the true N. Also the maximum likelihood 
standard deviation appears to represent an upper bound to the precision 
which can be attained, even if the assumptions underlying the trapping 
procedure are realized in practice. 

Regression estimates of N = 196 were made for the 162 samples 
and agree closely with the maximum likelihood estimates. The mean 
and standard deviation of the regression estimates are 199.4 and 23.5, 


respectively. 
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TABLE 2 
Maximum Likelihood Estimates 


Parameters N = 98 N = 196 N = 392 

Number of estimates 164 163* 162 160 
Mean 109.4 103.3 199.9 395.4 
Observed 8.D. 81.5 22.4 24.0 29.3 
Theoretical S.D. 14.5 19.5 26.9 
Per cent within N + 1 theor. 

73 75 70 
Per cent within N + 2 theor. 

S.D.’s 92 93 94 


*Sample 72 omitted 


Summary of Sampling Results 


The sampling experiments in this section demonstrate moderately 
close agreement between maximum likelihood and unweighted least 
squares regression estimates of population size made from three trap- 
pings. With the aid of Figures 1 and 2, the maximum likelihood esti- 
mates may be obtained more rapidly than the least squares regression 
estimates. 

Comparisons were made between the standard errors of the sampling 
distributions of N’, at various levels of N, and the corresponding asymp- 
totic maximum likelihood standard errors. These comparisons show 
that for three trappings and a probability of capture of about .4, reason- 
able conformance between the asymptotic and observed standard errors 
takes place for N of about 200 or more. Underestimation of the true 
standard error occurs when the asymptotic formula is applied to data 
from smaller populations. For N of less than 200 an interval of two 
asymptotic standard errors (formula given and applied in example 
earlier in this section) about NY may be considered to represent approxi- 
mately a 90% confidence interval for N rather than a 95% confidence 
interval. 


Ill. EFFECT OF DIFFERENT TRAPPING SCHEDULES ON THE PRECISION 
OF MAXIMUM LIKELIHOOD ESTIMATES 


The biologist is interested in the proportion of the total population 
that he must trap in order to be fairly certain that his population esti- 
mate lies within some specified per cent of the true value. Also should 
he trap intensively over a short period of time or apply less effort over 
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a greater number of trapping intervals? In studying the latter question, 
let us assume that the trapping region is fixed and that the home ranges 
of the animals overlap within this trapping area. Thus the number of 
animals exposed to capture during the first trapping is the same regard- 
less of the number of traps set. 

Let po represent the probability that an animal will be captured 
by a given trap during a single period of trapping, and let gq = 1 — po. 
Assuming p> to be the same for each trap and also assuming that each 
trap operates independently of all other traps, then the probability 
that an animal escapes trapping by any of ¢ traps in a single period 
is go = g. Also g* = qo‘ is the probability of not being caught in any 
of k periods. Consequently the proportion of the population which is 
expected to be captured throughout the trapping program is theoretically 
a function only of go and of the total number, tk, of trap-nights. 

We may now see how the variance of the estimates of N is affected 
by altering the trapping schedule, i.e., the values of ¢t and k, when the 
value of go and the total expected capture N(1 — q*) are kept constant. 
Equation (15) may be written in the following form: 


1 
(CP) 


(20) 


V(N) = 


where 


aN? oN? 


For a given total expected capture, the quantity Ng‘ remains con- 
stant, and consequently F’ remains constant. Thus the effect on 
V(N) of different schedules will be due to the changes in the value of 
(kp)’?/q as k and p vary but (1 — q*) remains fixed. 

Table 3 shows the standard deviation of NY when N equals 200, 300, 
500 and 1000, and when given proportions of the population are expected 
to be captured over varying numbers of trappings. The standard 
deviation decreases only slightly as the number of trapping intervals 
increases beyond three. For practical purposes the standard deviation 
can be considered essentially constant for a given N and a given total 
proportion captured. Therefore, under the assumption that infiltration 
of animals from nearby areas is of negligible importance, the trapper 
would choose that schedule which would minimize his cost. In some 
circumstances he might find it more economical to set many traps over 
a short period; in others, a few traps for a considerable length of time. 
If infiltration were expected to disrupt the assumptions of the method 
of estimation, a short intensive program would be indicated. 
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Discussion of Precision 


With the preceding results in mind, i.e., that the precision of the 
estimates depends primarily upon the expected proportion of the total 
population captured, we may now study the removal method to see 
approximately what proportion should be captured over an entire 
trapping program in order to reduce the standara Jeviation to a specified 
per cent of the population size. Table 4 gives these values for population 


TABLE 4 


Proportion of Total Population Required to be Trapped for Specified Coefficient 
of Variation of N 


Coefficient of Variation 


N 30% 20% 10% 5% 


Proportion (to nearest .05) of population to 
be captured (in 100 or fewer trappings) 


200 55 60 90 
300 50 60 75 
1,000 40 45 60 
10,000 20 2 50 
100,000 10 15 .20 .30 


sizes from 200 to 100,000 and for coefficients of variation of 30, 20, 10, 
and 5 per cent. This table shows that relatively large proportions of 
the population must be trapped in order to obtain precise estimates. 
It must also be remembered that the variances used in these calculations 
are at best valid only when the assumptions of a constant binomial 
probability of capture and a stationary population are satisfied. Since 
these assumptions are unlikely to be completely realized under actual 
trapping conditions, we should consider the variances obtained above 
as a minimum that might be reached under ideal conditions. The poor 
precision of the removal method raises a serious question as to its 
utility to the biologist. 

Changes in population density that are large may be detected by 
the removal method. Let us assume that at times 7’; and T, the popula- 
tion is of sizes N and rN respectively, i.e. r is the factor by which the 
population size has been changed during the interval between 7’, and T, . 
Assume that at time 7’, , the intensity of trapping (i.e. the number of 
traps set) is such that the coefficient of variation of the estimate of 
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population size is some value c. If the same number of traps are set at 
time 7’, , and are sufficiently numerous so that the possible increase in 
the population size would not affect the assumption of independent 
action of each trap, then theoretically the probability of capture during 
a trapping at time 7’, is the same as at time 7',. From equation (16), it 
is seen that asymptotically for given values of p and k, the variance of an 
estimate of population size is proportional to the population size; hence 
the coefficient of variation of the population estimate at time T’, is c/ Vr. 
Values of r, the amount of change in population size, which have a 
specified probability of being detected on the basis of trapping results 
at times 7, and T, , can therefore be determined. The following table 
gives the size of population change that has an eighty per cent chance 
of being detected (at the 5% level of significance) from trapping results, 
in which the trapping effort in the two trials is constant and is of sufficient 
intensity so that the coefficient of variation of the first estimate is a 
specified per cent of the population size at that time: 


Downward 
Coefficient of Upward Change Change 
Variation (Per Cent) (Per Cent) 
1 44 36 
2 96 65 
3 159 89 
4 233 — 


1V. VALIDITY OF ASSUMPTIONS 


In addition to knowing the theoretical properties of the removal 
method, the degree to which the assumptions underlying the method 
are met in practice will help determine its usefulness as a tool for obtain- 
ing estimates of population size. 

If fairly intensive trapping is employed over a short period of time, 
the assumption of a stationary population may be considered reasonable 
in many situations. The removal method also assumes that the prob- 
ability of capture during a trapping is the same for all animals, and 
that this probability does not change from trapping to trapping. Evi- 
dence of trap-proneness and shyness is reported in the literature (Chitty 
and Kempson, 1952; Tanaka, 1951; Young eé al, 1952); hence it would be 
valuable to have a quantitative measure of the amount of variation in 
eatchability within an animal population. It should be pointed out 
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that any population estimate based on trapping results applies only 
to the catchable population and therefore in situations where for 
example the young remain relatively inactive and not subject to capture, 
the adult population alone will be estimated. 

In addition to individual differences during a given trapping it 
would be important to know the changes in catchability that occur from 
trapping to trapping. Weather changes for example will influence 
amount of activity, which, in turn, will affect the probability of capture 
(Brown, 1954; Calhoun, 1945). Also, how much of a role does learning 
play as a factor in the probability of capture? Ecologists may help 
answer questions such as these in the future. 

The removal method has been used in the North American Census 
of Small Mammals (Calhoun, 1948-1951). Using the data from these 
censuses, two different chi-square tests were employed in order to study 
the assumption that the probability of capture remains constant 
throughout the trapping program. 


Chi-Square Tests 


In Section II, it is shown that a multinomial approach may be used 
as an alternative to the conditional binomial approach in obtaining 
maximum likelihood estimates of N, the population size, and p, the 
binomial probability of capture assumed constant throughout the trap- 
ping program. The multinomial approach assumes that the observed 
catches follow a multinomial distribution with the probability of being 
captured during the ith trapping equal to pg**. Consequently the 
probability that an animal will not be captured during the entire pro- 
gram of k trappings is g‘. If this hypothesis is satisfied and N and p 
are known then, asymptotically, the quantity 


(yi — Y , W-T— Nd)’ 
+ 


will be distributed as an with k degrees of freedom. In expres- 
sion (21) y; represents the ith observed catch and 7’ equals the total 
of the individual catches. Hence to test the hypothesis that the prob- 
ability of capture has a constant value p, X} may be calculated for a 
given set of data and compared with the critical tabular value. The 
asymptotic properties of this multinomial or unconditional chi-square 
test, including its limiting power, are known (see e.g. Cochran (1952)). 

If either or both of the parameters N and p are estimated from a 
sample, then expression (21) with the estimates substituted is also 
distributed as chi-square with decrease in degrees of freedom correspond- 
ing to the number of parameters estimated. When both N and p are 
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estimated and 7'/(1 — 4) is substituted for N (equation 8), then the 
form of the chi-square test becomes 


with (k — 2) degrees of freedom. This expression, which compares the 
ith observed catch with the proportion of the total catch expected during 
the ith trapping, is independent of N and therefore is a test of the 
assumption that the proportion of free animals expected to be captured 
during each trapping is constant. 

A conditional test of a constant binomial probability of capture 
suggests itself if consideration is given to x; , the total number captured 
prior to the ith trapping. By analogy with a chi-square test for the 
sum of a number of independent binomial variables, the conditional 
binomial test is of the form 


k 2 

— (N — 
The properties of this test, such as its power function, have not been 
thoroughly studied. However it may be shown that when the assump- 
tion of a constant binomial probability of capture p is satisfied, expres- 


sions (21) and (22) yield asymptotically identical results for a given set of 
data. To do this let us define the ith catch as 


y, = (N —2)pt+eV(N — (23) 


where e; is a normal random variable with mean 0 and variance 1, and 2; 
represents the total capture before the ith trapping. 

Substituting the right-hand side of (23) for y; in expression (21) we 
find, after considerable algebra, that 


k 
X? = Dé? + terms of order N~””. 


The same procedure applied to expression (22) yields 


k 
Xi= Dé. 
t=1 
Thus, asymptotically, for a given set of data the multinomial and 
the conditional binomial tests for the constancy of p provide numerically 
identical results. Since the multinomial test is known to follow a chi- 
square distribution, this result shows that the conditional binomial test 
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also follows a chi-square distribution with the same number of degrees 
of freedom. It may further be assumed that the conditional binomial 
test using estimates of N and p substituted for the parameters will also 
be distributed as chi-square with (k — 2) degrees of freedom. 

It should be pointed out that the conditional binomial test is not 
expected to detect changes in N, since it may be shown that a linear 
regression of catch, y; , on previous catch, x; , is expected when the 
exposed population changes by a constant factor between trappings 
and p remains constant. 

It would be of considerable interest to determine the relative powers 
of the conditional and unconditional tests of what is essentially the same 
hypothesis. This problem was not taken up actively, however, since it 
represents a digression from the main purpose of this paper. 


Application of Chi-Square Tests to Trapping Data 


In the North American Census of Small Mammals effort was made 
to sample rodent populations of the genera Microtus, Sigmodon, 
Peromyscus, Synoptomys, and Dicrostonyx. Trapping data were 
collected in many locations throughout the United States and Canada. 

Conditional and unconditional chi-square tests of the assumption 
of a constant probability of capture were performed on 149 sets of these 
trapping data. Each set of data represented the catches in a given 
area on three successive nights of trapping, a constant number of traps 
being set each night. The data presented here include only those cases 
where the total of three catches in each set equals or exceeds ten animals, 
and also where the first night’s catch exceeds the final catch in order 
that maximum likelihood estimates of N and p might be made. Fifty- 
five other sets, in which the second condition was not met, were excluded 
from the calculations. The catches were not separated according to 
genera since this would have resulted in groups of catches too small to 
analyze. Therefore the hypothesis tested is in effect that the probability 
of capture is the same for all the genera and remains constant over the 
entire trapping program. The results of the chi-square calculations are 
given below: 


Result Conditional Unconditional 
chi-squares chi-squares 
Significant at 5% level 9.4% 7.4% 
Total of 149 chi-squares 192.20 187.57 


P of total chi-square with 149 d. f. .009 .016 
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Thus both total chi-squares show significance at the 2% level of 
significance. 

A trend in the direction of the deviations of observed from expected 
catches was studied as a possible contributor to the significant total 
chi-squares. A comparison between the observed third catch and the 
expected third catch, as predicted from the first two observations, 
provides the following results: 


Ys — E(ys) 
Number | Positive | Zero | Negative Total 
No. of samples | 73 | 5 | 71 149 


Although these 149 samples are selected (y, required to exceed y;), 
consistently positive differences would be expected if major infiltration 
during the trapping programs were a general phenomenon. Conversely, 
if the probability of capture were decreasing throughout the trapping 
program, consistently negative differences might be observed. The 
above results do not suggest any such trend in the direction of the 
differences. Consequently, the significant chi-squares appear to repre- 
sent unsystematic deviations that are on the whole larger than simple 
binomial variability would predict. This may be related to differences 
among the species comprising the catches or this added variability 
might be due to any of a number of additional factors (climate, food 
supply, ete.) which the theory assumes to remain constant throughout 
the trapping program, but which may change in practice. Some varia- 
tions in trapping effort regularly occur due to variations in the number of 
sprung but unfilled traps observed throughout a trapping program. 

The formulas for the standard deviation of N given earlier in this 
paper assume that the size of each catch is subject only to binomial 
variability. The data studied above suggest that additional variability 
is present, which will be underestimated by the theoretical standard 
deviations. As an approximation to the true standard deviation, one 
might adjust the theoretical value by the factor 


total chi-square 
degrees of freedom 


which is about 1.2 for the trapping data studied. This kind of rough 
adjustment for variability in excess of purely binomial variability has 
been employed in the field of bioassay (Bliss, 1952). 
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V. SUMMARY 


The theory of the maximum likelihood method of estimation of 
population size has been reviewed and developed for application to data 
collected in a trapping program employing constant effort and removal of 
captured animals from the population. The assumptions underlying this 
method are that the population is stationary, that the probability of 
capture during a given trapping is the same for all animals exposed to 
capture, and that this probability of capture does not change from 
trapping to trapping. 

A rapid graphical method for obtaining maximum likelihood esti- 
mates of N (population size) and p (probability of capture during a single 
trapping) has been presented as well as formulas for the standard 
errors of these estimates. 

The asymptotic precision of the removal method of estimation has 
been determined for specified levels of probability of capture and for 
given numbers of trappings. A study was made of the effect on the 
asymptotic variance of maintaining constant the total number of units 
of trapping effort, but varying the trapping schedule. It was seen 
(Table 3) that only a slight gain in precision is expected from increasing 
the number of trappings with fewer traps set per trapping. The propor- 
tion of the population that must be captured for a specified coefficient of 
variation has been calculated. The precision of the method was shown 
to be poor, e.g. for N of 1000 or less, at least 40% of the total population 
must be trapped in order to obtain a coefficient of variation of 30% of less. 

Experimental sampling was carried out in order to compare the maxi- 
mum likelihood and regression methods of estimation and to compare the 
distributions of observed estimates of N with the asymptotic distribu- 
tions. The distributions of estimates obtained by the two methods were 
very similar. For small N’s the observed distributions were skewed in 
the direction of overestimating N, and the estimated variances exceeded 
the asymptotic values. With a probability of capture of .4 during each 
of three ‘‘trappings”, fairly close approximation to normality and the 
asymptotic variance was attained by the distributions of estimates of 
N = 200 and larger. 

Chi-square tests of the assumption of a constant binomial probability 
of capture were performed on trapping data collected in the North 
American Census of Small Mammals. Both conditional and uncon- 
ditional chi-squares were significant at the two per cent level. A study 
of the direction of the differences between observed and expected catches 
revealed no trend, suggesting that the significant chi-squares might have 
resulted from larger unsystematic differences between observed and ex- 
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pected than would be explained by simple binomial variability. Differ- 
ences in catchability among the animals and/or other factors, assumed 
constant by the theory, might explain the additional variability observed. 
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THE CONCEPT OF PATH COEFFICIENT AND ITS 
IMPACT ON POPULATION GENETICS* 


Cc. C. li 
Graduate School of Public Health, University of Pittsburgh 


INTRODUCTION 


The method of path coefficients was first published by Professor 
Sewall Wright thirty-five years ago. In 1921 there appeared in the 
Journal of Agricultural Research (Vol. 20) a general account of the method 
and the relationship between correlation and path coefficients, together 
with some examples of application; and in Genetics (Vol. 6) a series of 
five papers dealing exclusively with the application of path coefficients 
to genetic problems. Previously known results of various mating 
systems, obtained by laborious arithmetical procedures, were confirmed 
by the more elegant method, and many new results were reached, some 
of which were later corroborated by the method of matrix algebra while 
others are still difficult to obtain by any other method today. These 
classical papers, together with the pioneer work of Fisher (1918), still 
constitute the basic readings for students of population genetics, although 
the method has since taken a more sophisticated form and the field of 
application has been widened. However, one must admit that the 
method of path coefficients, as powerful and flexible as it is, was not 
immediately very popular among geneticists, still less so among pro- 
fessional statisticians. It was much later that its usefulness became 
gradually and generally appreciated. 

Path coefficients can be treated at various mathematical levels. 
The most important properties, however, can be deduced and studied 
by standard statistical tools. To understand the method requires little 
more than a knowledge of multiple and partial regression and correlation. 
It is a special type of multi-variate analysis—a method of dealing with a 
“closed” system of variables that are linearly related. (For non- 
linearly related variables, an appropriate transformation of scales may 


*Presented at the A.I.B.S. symposium on Sewall Wright’s contributions to Population Genetics, 
sponsored by the Genetics Society of America, the American Society of Human Genetics, and the 
Biometric Society (ENAR), and held in East Lansing, Michigan, September 6, 1955, 
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be necessary for its applicability.) By a closed linear system is meant 
that each variable in the system is either a linear combination of some 
other variables in the system or is one of the basic factors, which may 
be correlated with or independent of other basic factors in the system. 
In other words, the system is formally complete, including all the basic 
factors (the “causes’”’) and their resultant variables (the “effects’’). 

On account of the nature of a closed system, the practical employ- 
ment of the method of path coefficients is greatly facilitated by the 
formulation of a diagram showing the interrelationships of the variables 
concerned. In constructing such a diagram the convention is to use a 
double-headed arrow to indicate correlations between basic factors and 
a uni-directional arrow to signify the direct path of influence from one 
variable to another. The diagram, or causal network, must be consistent 
with a particular viewpoint and contain no duplicate or superfluous 
paths of connections. It is important to realize that various different 
viewpoints may be taken in regard to a system of interrelated variables, 
and, consequently, the causal diagram may be constructed in various 
ways. At least this is so for a purely mathematical setup. The choice 
of a particular viewpoint is up to the investigator, but, once it is taken, 
it must be held consistently in the entire diagram. 

At this point it is well to clarify a possible misapprehension of the path 
method. It is not an endeavor to infer causal relations from observed 
correlations among a set of interrelated variables. Quite on the con- 
trary, the employment of this method must be preceded by the formu- 
lation of a causal scheme, either based upon an a priori knowledge of 
the causal relations or based upon an hypothesis which the investigator 
chooses to accept or test. Consequently, the more we know of the true 
relations among the variables, the more meaningful will be the results 
of path analysis. In the case of incomplete external knowledge con- 
cerning the causal relations, the interpretations derived from the analysis 
are of course subject to revision in the light of further knowledge and 
revised hypotheses. 


PATH COEFFICIENTS 


A comprehensive mathematical presentation of the method of path 
coefficients was given by Professor Wright himself in 1934. His two 
more recent summaries of the subject are to be found in Annals of 
Eugenics (1951, Appendix) and Statistics and Mathematics in Biology 
(1954, Chap. 2). The only “pure” statistician (known to the writer as 
of now) who made a thorough discussion of the biometric concept and 
limitations of the method, as well as some topics related to path co- 
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efficients, is Prof. John W. Tukey (1954, Chap. 3), whose work should 
be consulted by both statisticians and geneticists. A much more 
elementary exposition of the subject has been outlined by the present 
writer (1955, Chap. 12). It is not proposed to reproduce all the mathe- 
matical details here. A very sketchy account must suffice for the present. 
Some of the more specific points will be discussed along with concrete 
examples in later sections. 

' The first premise is that the degree or extent of influence of one 
variable upon another can be expressed in quantitative terms. Then 
it is a matter of devising a numerical measurement of such an influence. 
After the construction of an adequate causal diagram referred to in the 
previous section, the task is to find a device for assigning a value to 
each of the arrows serving as a symbol of influence along that path 
(which is directional). The value so assigned to the path is called a 
“path coefficient”’. 

Suppose that a dependent variable Y is linearly determined by two 
factors, X and Z. This is our causal scheme, the simplest of its kind. 
In words, the path coefficient for the arrow from X to Y is defined as 
the portion of the standard deviation of Y that is due to the variation in 
X. The meaning of this verbal statement may be made clear by the 
following consideration. Let the ordinary multiple regression equation 
of Y on X and Z (fitted by method of least squares) be Y = A + 
BX + CZ, or, more conveniently, 


= — ¥) %) (1) 


where B is the regression coefficient of Y on X. It is the number of 
units expected to change in Y for each unit change in X, such as the 
change in number of bushels (Y) of grain per acre for each inch of 
rainfall (X) in the growing season, or the number of pounds expected to 
change in weight Y for each inch of height X. Clearly, the value of B 
not only depends upon the degree of influence of X on Y but also depends 
upon the actual physical units employed in measuring the variables. In 
other words, the value of B will change if the Y is expressed in units of 
tons instead of pounds and/or X is expressed in units of centimeters 
instead of inches, although the influence of X on Y remains the same. 
It is therefore desirable to devise a regression equation expressing the 
same relationship but independent of physical units. This may be 
accomplished by using the so-called “standardized”’ variables. Thus, 
equation (1) may be rewritten 


Oy Oy oy oz 


(2) 


Ox 


{ 
a 
| 
tie 
| / 


PATIL COLFFICIENT 193 


To be sure, the equations (1) and (2) are saying the same thing, except 
that (1) is in terms of actual units and (2) in units of standard deviations. 

It is convenient here to use lower case letters to denote the corre- 
sponding standardized variables. Writing x for (X — X)/cx etc. the 
foregoing equation takes the form 


y = br + cz (2’) 

in which the new regression coefficients for the standardized variables are 
B 

b = Pyx = c= Pryz = (3) 


These regression coefficients are also without physical units. The value of 
b = pyx is named the path coefficient for the direct path from X to Y. 
Since a path coefficient is merely a regression coefficient in standardized 
form, it may also be referred to as a “‘standardized regression coefficient”. 
Note that the definition (3) gives a precise meaning of our previous verbal 
definition: the portion of the standard deviation of Y that is due to the 
variation in X,. 

The path coefficient, so defined, possesses many properties which 
make it useful in statistical analysis. Being a type of regression co- 
efficient, it is directional (e.g. from X to Y), may be positive or negative, 
and may be greater or less than unity. Being without a physical unit, 
it resembles a correlation coefficient. Indeed, it reduces to an ordinary 
correlation coefficient under certain simple conditions. Some of these 
properties are illustrated in the following examples. 


EXAMPLES OF ANALYSIS 


Perhaps there is no single numerical example which would give all 
the properties of path coefficients. The basic purpose of this section is 
to show the similarity and difference between path analysis and the 
ordinary multiple regression method, and to point out the advantages of 
the former over the latter in interpreting and measuring the causal 
relationships. To begin with, let us consider the two sets of data pre- 
sented in Tables 1 and 2, and assume that Y is a linear combination of 
X and Z in both cases. Fitting a regression equation of type (1) by the 
usual procedure of least squares amounts to the solution for A, B, and 
C of the following set of ‘normal’ equations: 
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where N = 8 = number of observations in each variable. The equations 
(4) may be found in almost every textbook of statistics and are repro- 
duced here for reference and later use. Thus, for the data of Table 1, 
we find that Y = —1.70 + .30X + .32Z, or 


Y—Y = — X) + 3X2 — (5) 


Proceeding the same way with the data of Table 2, we find that A = 
—1.7, B = .30, C = .32; that is, the regression equations for the two 
sets of data are exactly the same. Based upon the regression equation 
(5) alone, we might conclude that the influences of X and Z on Y are 
the same for the two cases. However, expressing the equation in the 
standardized form (2’), we have for the data of Tables 1 and 2 respec- 
tively: 
= 60c + .80z, (5.1) 
y = 4932x2 + .6576z. (5.2) 


Now we proceed to examine why (5.1) and (5.2) are different and what 
the standarized regression coefficients mean. The chief difference 
between the two sets of data, as may have been noticed, is that X and Z 
are uncorrelated in the first example and correlated in the second. 
The results of analysis up to this point may be best summarized in the 
form of diagrams. Thus, Figures 1 and 2 represent the causal relations 


FIGURE 1 FIGURE 2 
Uncorrelated causes Correlated causes 


embodied in Tables 1 and 2 respectively. The coefficients in (5.1) and 
(5.2) are the path coefficients of the respective arrows in the diagrams. 

The path coefficient is taken as a measure of the degree of influence 
along the path. We note that when the causal factors are uncorrelated 
(Fig. 1), the path coefficient is simply the ordinary correlation coefficient 
between the two variables concerned (as given at the bottom of Table 1). 
Thus, 


= -60 (6.1) 
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When the causal factors are correlated (Fig. 2), the following relations 
hold: 


fyx = Pyx + TxzPyz = .4932 + (.5).6576 = .8220 (6 2) 


fyz = + TxzPyx = .6576 + (.5).4932 = .9042 


The meaning of the relations (6.2) becomes immediately clear upon 
examining the corresponding diagram. For instance, the first expression 
shows that the correlation between X and Y is the sum of the values of 
the two paths connecting them. In other words, the total correlation 
Tyx = .8220 has been separated into two components: the component 
Pyx = .4932 measures the influence of the direct path from X to Y, 
; and the other component rxzpyz = (.5) (.6576) = .3288 measures the 

influence of the indirect path from X to Y, via Z. The route X + Z— Y 
is also known as a “compound path’’, whose value is the product of 
< those of the two individual steps. The separation of a correlation co- 
efficient into various components is one of the chief accomplishments of 
the method of path coefficients. Analogous to “the analysis of variance’, 
the path method may be called “the analysis of correlation’’. 

Since the relation (6.1) is a special case of (6.2), we need only to 
demonstrate the truth of the latter. On transforming the actual values 
Y, X, Z, into the standardized values y, x, z, we obtain the pleasingly 
simple results: = Dor = = 0, = = o2 = 1, and 


Cov (y, x) = Ty. = Trx; etc. 


From the last expression, we see that a correlation coefficient may be 
called a “‘standardized covariance”. To fit a regression equation of y 
on x and 2, we need only substitute these simplifying values in normal 
equations (4) in which we use b and ¢ to stand for the path coefficients. 
The first equation yields a = 0. The remaining two equations reduce to: 


tyx = b+erxz 


= brxz 


(6) 


which is our (6.2) with b = pyx andc = pyz. Whenryz = 0, we obtain 
5 (6.1). The solution for b and c from equations (6) will of course yield 
$ the standardized regression equations (5.1) or (5.2). 
. There remains one more point to be mentioned. If Y is completely 
(without “residual error’) determined by X and Z according to relation 
(1), then 


oy = Brox + + 2BCoxozrxz (7) 
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Dividing by the variance of Y throughout, 
Oy Oy oy oy 


l= Pyx + Dyz + 2pyxPyzixz - (8) 


The values of Y in our numerical examples are completely determined by 
X and Z. Hence, for Figures 1 and 2, we have respectively 


(.60)* + (.80)’ = 1, 
(4932)? + (.6576)? + 2(.4982)(.6576)(.50) = 1, 


The terms of (8) are called the ‘coefficients of determination”. Thus, 
for Figure 1, we say that (.60)? = 36% of Y’s variation is determined 
by the variation in X; this conforms with the usual statement that o7 
will be reduced to — when X is “held constant” in the ordinary 
linear regression analysis of Y on X. In Figure 2 we say that X directly 
determines (.4932)” = 24.32% of Y’s variation while the joint effect of 
X and Z determines 2(.4932) (.6576)(.50) = 32.43% of Y’s variation. 
When the causal factors are positively correlated, the expression (8) 
for determination shows that the higher the correlation, the larger the 
joint effect of X and Z, and consequently the smaller the direct effects 
of X and Z separately. When rxz is negative, there is a semantic 
difficulty in speaking of a negative percentage. The meaning is never- 
theless clear: the negative correlation diminishes the variance of Y, 
making it smaller than it would be if X and Z were uncorrelated. 

General formulas given by Wright (1955), which are of the type of our 
(6), may be obtained immediately by writing down the standardized 
form of “normal” equations of type (4). 


or 


(8’) 


THE CORRELATION COMPONENTS 


The most direct application of the path method is the deduction of 
the correlation between two variables (say X and Y), which are linear 
functions of some common variables. For example, suppose 


X = Bo + B,Z, + B.Z, + BZ, + E, 
Y = Co + + C.Z, + + E, 
Upon standardization, the two linear expressions above become 


x= dizi + deze + + pier 


Y = C2, + C22 + + 
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Now assume that the two E’s are neither correlated with each other nor 
with any of the Z’s. Then the correlation between X and Y is 
Txy = cov (x, y) = + doze + + C2%2 + C523) 
= be, + + Dirises 
+ + deta, + der 


Translating the expression into a diagrammatic form (Fig. 3), we see 
that rxy is the sum of the nine compound paths by which X and Y are 


(9) 


FIGURE 3 
The correlation between X and Y due to common causes Z1, Z:, 73. 


connected. In the language of path analysis, we say that the three 
Z’s are the “common causes” of X and Y. The latter are correlated 
because of the presence of the common influences of the.Z’s. Some of 
the terms of (9) may drop out if any two of the common causes are 
uncorrelated between themselves. In particular, when all of the common 
causes are “‘independent”’ of each other, we simply have 


rxy = bic, + + dye; . 


Writing px, for b, to denote the path coefficient from the common 
cause Z, to X, etc., the general form of (9) may be written 


= PxiPyi + Px iPiiPyi (10) 
i4i 


E; 
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where « = a common cause. This is a fundamental theorem in the 
theory of path coefficients. In words, briefly, the correlation between 
two variables is the sum of all the paths connecting them. 

When this same principle is applied to genetics, we will obtain the 
correlation between relatives, for the true or “‘blood”’ relatives are simply 
those sharing one or more common ancestors, who constitute the common 
“‘causes”’ for the heredities of the relatives. This is the reason why the 
- method of path coefficient has been so successful not only in calculating 
the genetic correlations but also in breaking it up into components, so 
that each component represents the contribution by a certain common 
ancestor through a particular line of descent. 


MENDELIAN VARIABLES 


It is fully realized that the above sketchy account does not do full 
justice to the method of path coefficients but it must suffice; and now we 
turn to its applications in population genetics. In the course of con- 
sidering genetic problems, some further properties may be elucidated. 
In the following discussion we shall ignore linkage, sex-linked, and 
polysomic inheritance, as well as environmental effects, and concentrate 
on the autosomal genes of bisexual diploid organisms. 

Consider a metrical character that is entirely controlled by a certain 
number of genes. We shall first assume that each gene has a specified 
amount of effect upon the trait, without dominance or interaction 
(epitasis) between loci. If the contributions of AA, Aa, aa to the trait 
are 2a, a, 0; and those of BB, Bb, bb are 28, 8, 0, respectively, etc. each 
independent of other loci, then the value of a Whole genotype is the sum 
of the effects of individual loci. Briefly, all genic effects are additive. 
We shall refer the value of the metrical character of the whole genotype 
as the ‘“‘genotypical value’”’ (of an individual). 

One has undoubtedly heard the statement something like this: a 
child is a half-and-half of the parents. Thus, we talk about one’s being 
half-Irish and half-German. In this sense the statement is understand- 
able. From the observation that a child often shows a feature intermedi- 
ate between those of his parents, the statement taken descriptively is 
acceptable. With respect to the fact that a child receives an equal 
number of chromosomes from each of his parents, the statement is even 
scientifically accurate. Despite all these, when we talk about a metrical 
trait, it may be surprising to many that the above statement is simply 
not true, or at least not very precise. The height, for instance (assuming 
height to be entirely controlled by additive genic effects), of an adult 
offspring is not necessarily equal to } father’s and } mother’s heights, 
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except in very special cases. Were every adult offspring to assume the 
mean height of his parents, then in a few generations all adults of a 
population would be of the same height! 

Then, what is the relationship between an offspring’s and his parents’ 
measurements? There are several statistical approaches to this problem; 
the most novel one is that employed by Wright. His method rests upon 
the introduction of gametic variables and the subsequent analysis by path 
coefficients. If we imagine that a (haploid) gamete assumes a value 
corresponding to its genic content, then the genotypical value of an 
individual is the sum of the values of the two gametes that united to 
produce the individual. Hence, we obtain the fundamental causal 
scheme: 


(11) 


where y is the genotypical value of an offspring and g, is the value of the 
gamete contributed by the father and g, that of the mother. The state- 
ment (11) is an exact one; it is not the same as saying that a child attains 
the mean value of his parents. Assuming additive genic effects and 
autosomal inheritance, y is completely, linearly and equally determined 
by g, and g,. Furthermore, we shall treat g, and g. and y as standardized 
variables. Then we see that (11) is a regression equation of the type 
(5.2) or (5.1), depending whether g, and g, are correlated or not. The 
path coefficients from the two g’s to y must be of the same magnitude, 
since the two gametes have equal influence on y. Following Wright’s 
original notation, let a denote the path coefficient from a uniting gamete 
to the resultant zygote. The theorem (8) on complete determination 
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FIGURE 4 FIGURE 5 FIGURE 6 
Panmixia Inbreeding Full sib correlation 


yields in this case: (Figs. 4 and 5, Left): 
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for uncorrelated g’s @+a@é=1, a=vVi 


(12) 
1 


where F is the correlation between the two g’s. Note that for uncorre- 
lated causes, a = V3 is also the correlation coefficient between y and 
one of the g’s. 

The truth of (12), like all others to be presented in this section, is 
independent of the trait under consideration, independent of the units 
of measurement, independent of the number of loci involved, and in- 
dependent of the gene frequencies of the population. 

Next, let us consider the path coefficient, b, from an individual to a 
gamete produced by him (Figs. 4 and 5, Right). The genic content of 
a gamete is of course limited by that of the individual, but within that 
limitation it is subject to chance of segregation in each meiosis. The 
correlation between an individual y and a gamete g produced by him is 
the same as the correlation between y and a gamete g’ of the previous 
generation that produced the individual. Hence (for details see Li, 
1955, p. 174), 


for correlated g’s 2a” + 2a°F 


for random individual, b=a= V} 


for an inbred individual, 6 = a’(1 + F’) = 9 


where the primes indicate the previous generation. 

It will be noticed that both a and b are expressed in terms of the 
correlation between uniting gametes. The statistic F plays a cardinal 
role in the path analysis of genetical problems of this sort. It is also 
known as the inbreeding coefficient of an individual (whose parental 
gametes have a correlation F). Of course, F cannot be directly observed 
since the gametes possess no metrical value in the real physical sense. 
It is purely conceptual, or, shall we say, it is a mathematical device to 
facilitate the analysis. This brings out an important feature of the 
general method of constructing causal diagrams and the subsequent path 
analysis: a variable pertinent to the causal scheme should be included 
in the diagram whether it is observable or not. 

The correlation between the two uniting gametes is due to the corre- 
lation between the two parents which actually could be observed and 
measured. In Fig. 4 the two parents are uncorrelated. In Fig. 5 the 
correlation between the parents (also called ‘“‘mates’’) is m. Hence: 


F = bmb = bm. (14) 
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Without spelling out all the details except to reiterate that a compound 
path value is equal to the product of single path values and that the 
correlation between two variables is the sum of all paths connecting 
them, we obtain the following results by examining Figs. 4, 5, 6: 


Correlation Panmixia Inbreeding 
Parent-parent: | rpp = QO, 'pp = ™ 
Parent-offspring:| rpg = Vi Vi =i. rpo = ab(1 + m) (15) 


Sib-sib: roo = + (V3)! = 3, lroo = + m) 


We started out with the simple causal scheme (11) between gametes 
and zygotes, and subsequently obtained four formulas, (12)-(15), 
which can be considered the four corner stones of population genetics. 
It should be clear, however, that all of these relations are consequences 
of the Mendelian mechanism in heredity, only expressed in a very 
concise manner. When every gene has a qualitative effect, the results 
are the various “‘ratios’’; when every gene has a quantitative effect, the 
results are the various correlations. Therefore, the set of parental, 
gametic, and offspring’s metrical measurements may well be called 
“Mendelian variables’. A Mendelian variable is a random variable 
conditioned by Mendelism. Although the method of path coefficients 
is applicable to a wide range of problems, it is particularly suitable for 
the analysis of Mendelian variables that have linear additive properties, 
owing to the clearcut ‘“‘cause-and-effect’’ relationship in heredity. 


SYSTEMS OF INBREEDING 


The subject of continued inbreeding of a Mendelian population has 
been a field of active research since World War I. In this section we 
shall not study the effects of inbreeding as such, but rather cite one or 
two examples to illustrate how the results of continued inbreeding may 
be obtained by using path coefficients in a way that is by far more 
economical than any other known method. 

Briefly, to obtain the results of continued inbreeding merely involves 
repeated application of the four fundamental formulas (12)-(15) estab- 
lished in the previous section. The gist of the whole trick is this: when 
a system (or pattern) of inbreeding is specified, there corresponds a 
causal diagram in which the correlation between mates (and thus 
parents of the next generation) may be expressed in terms of path and 
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correlation coefficients of the previous generation. This will finally 
lead to a recurrence formula for F between the successive generations. 
For example, continued brother-sister mating, the most extensively 
studied inbreeding system by statisticians, geneticists, and animal 
breeders, may be summed up by one expression: m = r%o (Fig. 7), where 


Too Too 
m In’ 
O™ © 
b b b b b b 
a’ a a a a a 
FIGURE 7 FIGURE 8 FIGURE 9 
Full sib mating Full sisters X half brother Double first cousin mating 
mating 
m = m = a’2b’2(1 + 2m’ + r’’o0) m = 2a’2b'2(m’ + r’’o0) 


roo = 2(a’'b’)?(1 + m’) according to (15). This single expression spells 
out the key fact of the system: the correlation between mates of one 
generation is the correlation between full sibs of the previous generation. 
Expressing everything in terms of F, F’, F’’, we immediately obtain 
the recurrence relation F = (})(1 + 2F’ + F’’), from which the cor- 
responding recurrence relation for heterozygosis may be deduced. Since 
the results of full sib mating are so well known, we shall examine another 
system of inbreeding in some detail. 

The mating system we propose to consider involves three mating 
individuals in each generation—one male (B) and two females (C, D). 
One mating (B X C, say) produces a son and the other mating (B x D) 
produces two daughters. The son (B,) is thus a half brother of the 
daughters (C, , D,) who are full sisters of each other (Fig. 8). In the 
next generation the matings will be B, X C, yielding a son and B, X D, 
yielding two daughters again, and so on. For brevity we shall call it 
the “full sisters X half brother’? mating system. The correlations 
between the three individuals of the previous generation are indicated 
at the top of Fig. 8. Now let us consider the correlation between mates 
of the following generation. The male and one of the females (e.g. 


concentrating on the one next to the male in the diagram) are connected 
by four chains: 
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(i) directly via their father, . . . . . . a’b’b’a’ 

(ii) via father and female's mother, . . . a’b’m’b’a’ 
(iii) via male’s mother and father,. . . . a’b’m’b’a’ 
(iv) via the two mothers, ...... . a’b’riib’a’ 


The correlation between mates, being the sum of these four compound 
path values, is thus 


m = + 2m’ + ri, (16) 


which is the key formula of the mating system. The rest is simply a 
matter of algebraic substitution. From (12), we have a” = 1/2(1 + F’); 
from (13), b? = (3)(1 + F’) and b” = (4)(1 + F’”’); combining (12) 
and (13), b’a” = ba’? = 3; from (14), b?m’ = F’ and b’?m” = F”; 
and from (15), rig = 2a’”b'"(m’’ + 1). Substituting these values in 
(16), we have: 


m = a” + 2b"? m’ + +- b’”)] 
F = b’'m = 31+ + + + (0+ 17) 
F = 3(3 + 8F’ + + F’”). 


The recurrence relation (17) represents the way the “full sisters 
xX half brother” system works on F. Each inbreeding system increases 
the values of F of successive generations in a peculiar way. The re- 
currence relation is an intrinsic property of the mating system, independ- 
ent of the initial condition of the population, number of loci involved, 
and the gene frequencies. This is another reason why F, the correlation 
between uniting gametes, plays such a dominant role in population 
genetics. Indeed, it is not too much to say that the major responsibility 
of path coefficients is to find the value of F in a population under various 
circumstances. 

There are many population properties closely related to F; but they 
are beyond the scope of this review. We must, however, mention the 
relation of F and H, the proportion of heterozygosis in a population. 
The meaning of H may be viewed in three different ways, with as many 
different physical meanings. First, with respect to one pair of loci, 
H denotes the proportion of heterozygotes (Aa) in the population. 
Second, n individuals each having k pairs of loci may be regarded as a 
population of nk pairs of genes, ignoring the entities of individuals. 
Then H denotes the proportion of heterozygous pairs of loci in the 
population. Third, considering the many pairs of loci present in an 
individual, H denotes the proportion of heterozygous pairs in the indi- 
vidual. All these three viewpoints are useful, and, fortunately, they 
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are mathematically equivalent. Whatever interpretation is adopted, 
H is proportional to 1 — F in any generation; thus, 


H=H(1-—F), =HA(1-—F’),  ete., 


where H, is a constant (the initial value of H before inbreeding). Sub- 
stituting in (17), we obtain the recurrence relation of H: 


H = 1H’ + 1H" + (18) 


which is independent of the initial condition H, . The relation shows 
that H decreases in each generation as inbreeding proceeds. For 
example, if the inbreeding system starts out with one male mated to 
two full sisters (unrelated to him) and the initial value of H is 3, we will 
have the following H series in successive generations: 


y. 2.4 8 13 23 40 139 242 421 


4’8’16’32’ 64’ 128’ 512’ 1042’ 2048’ 


H/H’: 1,1, .8125, .8846, .8696, .8687, .8705, .8698, --- 


After a sufficiently large number of generations, the ratio H/H’ will 
approach to a constant value, implying that the heterozygosis will 
thereafter decrease at a constant rate. To find this constant eventual 
rate of decrease, we put H/H' = H'/H" = Xin (18), which then becomes 

The largest positive root of this equation is \ = .86995 = .870. That 
is, the heterozygosis decreases by 13% per generation under this in- 
breeding system. 

The results of continued double first cousin mating may be obtained 
in exactly the same manner. The six correlations between the four 
individuals of the previous generation are indicated at the top of Fig. 9. 
The mates are also connected by four chains. The key formula is given - 
at the bottom of that diagram. The whole analysis takes but four steps 
analogous to (16), (17), (18), (19). These and previous results have 
been summarized in Table 3, in which the size N refers to the number of 
mating individuals involved in each generation. 

Finally, perhaps a few words should be said about the method of 
matrix algebra (brief account in Li, 1955, pp. 108-112, 116-118). Without 
doubt the matrix method gives a fuller account of the inbreeding process. 
For instance, with respect to one locus, it gives the frequencies of the 
various types of mating in each generation. However, it becomes 
difficult to handle when applied to a not-so-simple mating system, 
although Fisher (1949) has displayed great skill in such cases. On the 
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TABLE 3 
Some Results of Continued Inbreeding 


Mating system Size Recurrence relation Limiting value 

N of H of H/H’ =x 
Selfing 1 H = 3H’ » = .500 
Full sibs 2 H = 3H’ + 3H” » = .809 
Full sisters X half brother 3 H = 3H'+ 3H" + 7;H’” »A = .870 
Double first cousins 4 H = 3H' + 3H” + 3H" » = .920 


other hand, the path method yields the most pertinent information on 
an inbreeding system in very short order but little else. 


PANMICTIC FINITE POPULATIONS 


The general pattern of analysis with respect to a random mating 
group of N, males and N, females is the same as that described in the 
previous section, viz., first obtain an expression for m, then deduce the 
recurrence relations of F and H, and finally solve for the limiting value 
of H/H’ = i. The important conclusion from this analysis is that 
heterozygosis in a finite population decreases, although maybe slowly, 
from generation to generation even though random mating is the rule. 
We shall not give all the details here except to say that the m in this 
case is the weighted average correlation for full sib and half sib matings 
as well as those between more remote relatives. Writing ¢ for 1/No + 
1/N, = (No + N,)/NoN, , which in general is a small positive fraction, 
the recurrence relation of H turns out to be 


H = (1 (20) 


and 
Vl6+ e}. (21) 


The recurrence relation (20) and the limiting ratio (21) are exact and 
do not involve approximations. As far as the writer’s knowledge goes, 
these exact results have never been worked out by any other method 
although certain approximations for certain special cases have. 

The original expression for percentage decrease per generation given 
by Wright (1931, p. 108) is equal to 1 — X, and the expression in his 
1951 review (p. 347) is \ — 1, a negative quantity. 

When N, and N, are of moderate size, the fraction ¢ will be so small 
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that V16 + = 4+ With this approximation, the proportional 
decrease per generation in H becomes 


In the special case in which there are equal numbers of males and females 
(No = N, = 4N) ina group of N individuals, the fraction « = 4/N, and 


| , ” , 


This is the basis of the frequently encountered statement that hetero- 
zygosis in a finite population decreases by 1/2N per generation. This last 
result is equivalent to the findings of Fisher (1930, p. 87) through a 
different method. After ¢ generations the heterozygosis will be e~‘/?” 
of the original value. Thus, putting e~‘”” = 4, we obtain t = 2N log, 2 
= 1.386N as the number of generations required to halve the hetero- 
zygosis proportion in a population of size N. It follows that a random 
mating population of a limited size will eventually, in the absence of 
other forces, attain complete homozygosis. This conclusion had a 
great impact in forming the modern theory of evolution and contri- 
buted much toward understanding natural populations. 


FURTHER APPLICATIONS IN GENETICS 


It does not seem practical to go much further into the subject. In 
the following we shall merely mention a few examples to illustrate the 
scope of generalization obtained by path method rather than discuss 
specific implications. 


1. Non-additive genic effects. As far as the proportion of hetero- 
zygosis is concerned, the above results are true whether there is domi- 
nance or not. Dominance only affects the various correlation values. 
The path coefficients can be applied to the case with dominance for 
unilineal relatives, but not for bilineal relatives (such as full sibs and 
double first cousins). Where the method applies, a path from genotype 
to phenotype may be introduced. 


2. Phenotypic assortative mating also leads to correlation between 
mates and also can be analyzed by path coefficients (Wright, 1921, III). 
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It must be admitted, however, that the analysis is a complicated one 
even with the help of path coefficients. A more recent treatment of this 
and some related subjects may be found in Reeve (1953). 


3. Equilibrium with inbreeding. A,population does not necessarily 
reach complete homozygosis if the inbreeding is only of a moderate 
degree, or if the correlation between mates stays at a certain constant 
level in each generation. When an equilibrium condition is reached, 
F = 3m/(1 — 4m), a’ = 4(1 — 3m), b* = 3/(1 — 4m), so that the path 
from a parent to an offspring is ba = } whatever the mating correlation. 

4. Sex-linked inheritance. The method of path coefficients can be 
applied to the analysis of sex-linked loci just as easily as autosomal loci 
(Wright, 1933a). Only slight modifications are needed: the path from 
father (heterogametic sex) to his gamete is b = 1 because there is no 
segregation for a haploid; and the path from the male gamete to daughter 
is a as usual, but that to son isa = 0. The path from mother (homo- 
gametic) to her gamete is b as usual, but that from her gamete to son is 
a = 1 since the female gamete completely determines the genotype of 
the son (illustrations in Li, 1955, p. 182, 185). With these modifications 
analogous results for systems of inbreeding and panmictic finite popu- 
lations have been obtained. Further examples may be found in Crow 
and Roberts (1950). The separation of sex-linked and autosomal genetic 
variances of a quantitative trait can also be easily achieved by the 
method of path analysis (Reeve, 1953). 


5. Polysomic inheritance. In 1938 the method was extended by 
Wright to the analysis of the effects of inbreeding in the case of poly- 
somic loci. He carried the results much farther than had been reached 
by the matrix method, and obtained a general recurrence relation of H 
for 2k-somic loci for some simple systems of inbreeding. 


6. Linkage. The lack of independent assortment between loci 
introduces difficulty in analysis by any method. However, the inter- 
ference of inbreeding with the recombination of linked genes has also 
been studied by Wright (1933b). 


7. Isolation by distance. Naturalists and ecologists have long been 
aware of the isolation effects of mere distance, even though the species is 
continuously distributed over a large area. The quantification of the 
problem and the subsequent analysis by Wright (1943, 1946, 1951 
Appendix F) is also based upon path coefficients and recurrence relations 
of F. The mathematics involved is unfortunately not as simple as we 
would like it to be, but the writer knows of no other treatment of the 
subject through a simpler method. 
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8. Environmental influence on a Mendelian variable may be arbi- 
trarily defined as the influences of all causes other than genetic.’ In 
path analysis a separate variable, assumed to be independent of heredity, 
may be introduced into the causal scheme so that the variations in a 
dependent variable may be completely accounted for. The causal 
scheme is then formally complete. 


9. Animal breeding. The studies of Wright have exerted much 
influence on animal breeding programs and helped in understanding 
some of the experimental results. Breeding, of course, is by no means a 
purely genetical problem. Those interested in this subject may refer 
to Lush (1949) and Lerner (1950). 


10. Theory of Evolution. It is entirely out of the scope of this paper 
to deal with Wright’s theory of evolution (which is quoted and discussed 
at length by Dobzhansky, 1951). Here it need only be said that path 
and inbreeding coefficients played a basic role in the early stages of the 
formulation of the theory. 
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ERRATUM 


OUTHWAITE, ANNE D. and A. RUTHERWOOD: Covariance 
analysis as an alternative to stratification in the control of gradients. 
Vol. 11, No. 4, (1955), 431-452. 


On page 433, Table II, Total >> x.y should read 17,474.8 instead of 
17,414.8. 
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ESTIMATION OF THE NUMBER OF DIFFERENT CLASSES 
IN A POPULATION* 


R. C. LEwontTiIn AND TimotHy Prout 


Genetics Faculty, North Carolina State College, Raleigh, N. C. and 
Division of Life Sciences, University of California, Riverside 

It is frequently of interest in a biological investigation to estimate 
the number of different kinds of objects in some discrete population. 
Two examples are the number of species in a plant or animal community 
and the number of gene loci on a chromosome. 

Clearly any estimation procedure will depend upon the frequency 
distribution of the various classes in the population. There are two 
types of distributions which are of general interest and which are easily 
amenable to estimation procedures. The first of these is that in which 
the classes have some unimodal distribution which can be specified in 
terms of the modal frequency. Examples are the Poisson or any discrete 
distribution approximating a normal or log normal distribution. Preston 
(1948) has dealt with the problem of the number of species in a com- 
munity and has derived an estimate based upon an approximate log- 
normal distribution of the classes. This method applies for any of 
the class of distributions specified above. 

A second case of importance is that in which the classes are uniformly 
distributed. It is this case which will be dealt with here, the derived 


estimate being applied to the estimation of the number of genes on a 
chromosome. 


THE MAXIMUM LIKELIHOOD ESTIMATOR 


It is assumed that a population contains n classes, all of which are 
equally represented. If a sample of N objects is taken from such a 
population it will be found that r; objects fall into the ith class. Then 


the probability of obtaining a given sample is the multinomial prob- 
ability 


*Contribution from the North Carolina Agriculture Experiment Station. Published with the 
approval of the Director as Paper No. 647 of the Journal Series. 
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Not all of the classes in the population will be represented in the sampie, 
however. Let us say k classes are represented at least once while 
(n — k) are not represented. This means simply that (n — k) of the 
r,’s in (1) are zero. Expression (1) then becomes 


! N 
(2) P= (4)’, 


TI r:! 


Not all samples are distinguishable from each other however. It is not 
known which set of k classes of the n classes in the population are found 
in the sample. There are 


n! 


(3) 


different samples containing k classes. In addition among the k classes 

observed it is not known which class is represented r; times. There are 
k! 

(4) To 


II m,! 


J 


ways in which k classes may be partitioned so that m; classes appear 7 
times. The probability of any observed or distinguishable sample is 
then by (2), (3) and (4) 


which is by definition the Likelihood function. 
Taking the logarithm of both sides of (5) and maximizing* with 
respect to n we get 


(6) 


The value of *% for which equation (6) holds differs by at most one, 
from the ML estimator of n. This expression may be solved for # 
using a trial and error method from a table of the sums of the reciprocals 
of the integers. Since such a table is not generally available, use may 
be made of the approximation 


N 


(7) ry >In (Feller p. 175) 


*Strictly, (5) is not differentiable with respect to n. However, if the factorials are expressed as their 
equivalent I’-function, then differentiated, and the result evaluated for integral values of n, expression 
(6) results. 
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whose solution requires only a table of natural logarithms or exponentials. 

There are several points of interest to be noted about this estimator. 
First, expression (6) or (7) is identical in form to the solution of the 
coupon collector’s problem as given by Feller (1950). This problem 
concerns the expected value of N for a given value of k, when n is known. 
That is, how large a sample must be drawn to have k out of n possible 
classes represented in the sample. The close relation of this problem 
to our estimation problem is clear. 

Second, any departure from equal representation of classes in the 
population makes (6) an underestimate of n. This is because the more 
frequent classes will appear more often in the sample, thus lowering the 
value of k. 

Third, if the sample size N is smaller than or equal to the number of 
classes in the population, there exists a finite probability that all of the . 
classes in the sample will be different. That is, k may be equal to N. 
Moreover the smaller is the ratio N/n, the greater is the probability 
that k equals N. Now an inspection of expression (7) shows that when 
k equals N, the estimate 7 is infinite. 

Last, the estimator (6) is sufficient. This can be seen as follows. 
The likelihood expression (5) is written as the product of two expressions. 
The left hand member in brackets is a function only of the observations, 
while the right hand member in brackets is a function of k andn. Then 
by definition k is sufficient estimate of n. It follows directly from the 
properties of the MLE that it too is sufficient. This is simply another 
way of saying that all of the information relating to n in the sample is 
contained in k. 


THE DISTRIBUTION OF THE NUMBER OF CLASSES IN THE SAMPLE 


To find the variance of the MLE it is necessary to find the expected 
value of k. In addition to its importance in deriving the variance, the 
expected value of k is itself of some interest as it is in a sense the reverse 
of the coupon collector’s problem. 

Feller (p. 69) has derived an exact distribution for the number of 
classes not represented in a sample when all the classes are equally 
represented in the population. If m is the number of missing classes, then 


v=0 n 


where all the notation has been previously defined by us. This exact 
distribution is unmanageable as it stands. However, Feller shows that 
this expression converges to a Poisson distribution. In this form 


af) 
ile 
Md 
ae 
4 


214 BIOMETRICS, JUNE 1956 


(9) =e" 
m! 

where 

(10) Xr ne 


The demonstration of equation (9) depends upon the fact that is 
bounded. Should n grow very large for some fixed value of N, d will 
increase without bound. This limitation on the usefulness of the Poisson 
approximation is of great importance to the problem of finding confidence 
intervals for n and we shall return to it when confidence interval esti- 
mation is discussed. 

Before expression (9) corresponds exactly to the model in which we 
are interested, it must be noted that for our purposes m cannot be equal 
to n. That is at least one class is always represented. The desired 
probability distribution of m is then 


P(m) 
(11) fm) = 1 — P(m = n) 
Now the expected value of m is 
1 
(12) E(m) = x m! F — P(m = | 


Letting m — 1 =z 


1 
(13) E(m) = 1—P@=n— 


Note that the expression in brackets is equal to 1. Then 


(14) E(m) = 
but since 

(15) m=n—k 
(16) E(k) n(i — 


Thus, we have derived an expression for the expected value of k. Ex- 
pression (16) may also be written as 


N _ 
which is virtually identical with our MLE and the solution of the coupon 
collectors problem. (17) itself may serve as an estimator of n. 
It is interesting that the validity of (16) is not affected by any 
assumption about A. If we allow n to grow very large in relation to N 
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in (16) it appears that the expected value of k is N. This is precisely 
what is to be expected. That is, if the sample size is very much smaller 
than n, all of the members in the sample will be different. 

Equations (12), (13) and (14) are simply the proof that the expected 
value of a variable with a Poisson distribution is the parameter \ of that 
distribution. This proof has been shown in detail here since a similar 
method will be used in the next section. 


VARIANCE OF THE MAXIMUM LIKELIHOOD ESTIMATOR 


In general the variance* of a Likelihood Estimator may be written as 


(18) = (Cramér, 1951, p. 500) 


where L(u) is the Likelihood function. Then 


n N -1 
E 
Now 
4 = 


nin + 1) 


for even moderately large k and n. Using result (11) and (15), the left 
hand member in brackets of (19) becomes 


(20) 


Setting m + 1 = 2, this takes the form 

i¢ en ( 1 ) 
which by reference to (13), (14) and (15) is seen to be simply 

1 = -N/ny 1- 

(23) [n ne ] ne 


From (23) and (19) the variance of the MLE is given by 


(24) 
N/n 
(1 n ) 


*Rigorously, the asymptotic variance, although for some estimates this can be shown to be the 
small sample variance as well. 
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It should be pointed out here that although this variance will give a 
standard error of the MLE, the relation of this standard error to actual 
confidence limits for n is not simple. The question of confidence limits 
will be explored in a later section. 


APPLICATION TO THE NUMBER OF LOCI ON A CHROMOSOME 


An important parameter in population genetics is the number of 
loci on a chromosome. It has, in addition, bearing on the basic problem 
of the structure of the chromosome and its relation to the structure of 
the genes. Many methods are available for estimating this quantity 
depending upon the criteria of “locus” used (see Herskowitz, 1950). 
The particular estimate most used in population genetics is that of the 
number of loci capable of lethal mutation, the assumption being that 
the mutation rate and selection coefficient are identical for all loci and 
that therefore in an effectively infinite population all loci will be equally 
represented by lethal mutations. 

The method used to determine the number of loci capable of lethal 
mutation is as follows. A number of organisms are sampled from a 
population. By appropriate genetic tests the details of which may be 
found in Wallace (1950) it is determined how many of these carry on 
one of their chromosomes a lethal gene. Since organisms in general 
possess a duplicate set of chromosomes, an individual carrying one 
chromosome with a lethal mutation will survive, the homologous normal 
chromosome serving to cover the lethal. Separate stocks are now made 
each of which derives from one of the original lethal bearing organisms. 
Each stock is then mated with every other stock so that N stocks will 
give rise to N(N — 1)/2 crosses. If two stocks both contain a lethal 
at the same position on the chromosome, that is at the same locus, then 
when these stocks are mated some of the offspring will have both of their 
chromosomes carrying the lethal and thus will die. On the other hand 
if the lethal genes contained in the two stocks are at different loci, none 
of the offspring will contain a lethal gene in double dose so that all the 
offspring survive. If two lethal genes are 21 the same locus they are 
said to be allelic. : 

The procedure described above provides three items of information. 
First it serves to estimate the proportion of all chromosomes in the 
population which contain lethals. Second, among the N lethal bearing 
chromosomes it shows how many loci are represented once, twice, 
thrice, and so on. Finally, as a corollary of the second, it gives k, the 
number of different loci (classes) present among the N lethals tested. 
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ESTIMATION OF CLASSES 
THE FREQUENCY OF ALLELISM ESTIMATE 


The unique feature of the genetic estimation problem just described 
lies in the inability to identify two lethals as being identical, or at the 
same locus, without the process of mating. In a probabilistic sense, the 
mating procedure is simply a method for choosing pairs of lethals. If 
two members of the pair are identical the pair will be lethal. If the two 
members are not identical, the pair will be viable. Using this point of 
view Wright and Dobzhansky (1941) have suggested an estimate of n 
in the following way. If p, is the probability that the ith locus is repre- 
sented in the sample by a lethal, then p; is the probability that a given 
pair of lethal chromosomes will both be lethal at the 7th locus. Since 
there are n loci all assumed to be equally frequent in the universe from 
which the sample is taken, the total probability that a given pair will 
consist of two identical members is 


(25) 
Then 
1 
2 i= 
(26) 


is an estimate of n, where 4 is the observed frequency of allelism, that is 
the observed proportion of all pairs sampled which contained two 
identical members. 

There are several points to be noted in connection with this estimator. 
Like the MLE estimator it may give an infinite estimate of n when 
N <n. In this case @ may be zero which exactly corresponds to the 
probability that k = N in the MLE. 

No matter what the sample size, n is upwardly biased since it is 
clear that 


1 

(27) 
but since 
(28) 

E(x) 
unless all the 2’s are identical 

{1 

(29) I () >n 


In addition, the frequency of allelism estimate is not sufficient. We 
have shown that k is a sufficient estimate of n. There is however no 
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one-to-one correspondence between a and k. Thus all of the available 
information is not contained in a. 

There are two experimental designs to which the frequency of allelism 
estimate may be applied. The first of these has been described in the 
previous section and will be denoted as the complete cross design. The 
second method will be referred to as the random cross design. It consists 
in dividing m lethals into two sub-groups of size m/2 and making m/2 
crosses or pairs, one member of each pair coming from each of the two 
groups. From a probabilistic standpoint this simply consists in choosing 
m/2 random pairs of lethals. 

While the frequency of allelism estimate is applicable to both of these 
designs, the estimate will have a different variance in each case. For 
the random cross method, the variance of the estimate may be derived 
as follows. If an estimate, ¢, is asymptotically normally distributed, 
then in large samples the variance of any function of ¢ is given by 


2 
(30) = (Cramér p. 354) 
For the specific case in question 
2 


But @ is the proportion of successes in m/2 independent Bernoulli trials 
since it is the proportion of random pairs which show allelism. Then 


(32) a) 


Finally from (31) and (32) 
2 — 
(33) 


which for moderately large values of n is approximately 
2 3 
(34) o; = = 


When the complete crossing scheme is used, however, the variance 
cannot be arrived at in such a simple manner. This is because the 
N(N — 1)/2 crosses are not independent. The outcome of these crosses 
will be a unique function of the distribution of lethals in the original 
sample of size N. There will be a positive covariance among the crosses, 
the value of which will depend upon the distribution of lethals in the 
sample. We have not been able to derive a satisfactory expression for 
the variance of # for this model. However, for an equal number of 
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pairs the random cross scheme will have a lower variance than the 
complete cross design. By equal number of pairs we mean that m in 
the random cross is equal to N(N — 1) in the complete design. Ex- 
pression (33), then, represents the lower limit of the variance of 7. 
Setting m = N(N — 1), it is clear that the frequency of allelism 


estimate is inefficient as compared with the MLE since the efficiency 
of % is defined as 


3 
2 
(35) > 
N/n iv 


which is always less than unity as can be seen from the expansion of 
the exponential. 

Again it should be pointed out that (35) is the efficiency of Wright’s 
estimate under the random crossing scheme as compared with the MLE 
under the complete crossing scheme when equal numbers of crosses are 
made for both. Should Wright’s estimate be applied to the complete 
cross design, it would presumably be even less efficient. 

Fig. 1 shows the efficiency of % for various values of N/n. 

Before the estimates can be compared for the number of loci on a 
chromosome, a correction must be made for the presence of more than 
one lethal mutation on a chromosome. It is not possible to distinguish 
a chromosome with more than one lethal from those with only one. 
An additional datum required for this correction is the frequency of 
lethal bearing chromosomes in the population from which the sample of 
lethals was derived. On the assumption that the number of lethal 
genes in a lethal chromosome has a Poisson distribution, Wright showed 
that his estimate takes the form 


where &@ = frequency of allelism of lethal genes 
A = frequency of allelism of lethal chromosomes 
Q = frequency of lethal chromosomes in the population 


This same correction factor may be applied to the Maximum Likelihood 
estimate where N and k in eq. (7) or (9) are calculated from: 


N 


(37) 
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FIGURE 1. Relative efficiency of the reciprocal of allelism estimate of the number of loci. The 
efficiency, E, is plotted against the ratio of N, the sample size, to n, the number of loci. 


where k’ = no. of different lethal chromosomes observed and N’ = 
size of sample (no. of chromosomes tested). 


CONFIDENCE INTERVAL ESTIMATION 


Confidence intervals for n are simple to calculate when the random 
cross scheme is employed. Since a, the frequency of allelism in the 
design, is the observed proportion of successes in m/2 independent 
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Bernoulli trials, upper and lower confidence limits for a are found by 
use of tables of the binominal sums, or else binominal probability paper. 
The upper and lower confidence limits for n are then simply the recipro- 
cals of the values for a. When the complete cross design is used, how- 
ever, this simple method does not apply. Unfortunately several authors, 
(Pavan & Knapp, 1954; Wallace, 1950) have given confidence limits 
based upon the random cross model, while actually employing the 
complete cross design. Such confidence intervals are clearly too small. 
One of the present authors (R.C.L.) is responsible for this error in the 
paper of Pavan & Knapp. What the correct limits are can only be 
stated when the distribution function of G@ under the complete cross 
design is demonstrated. 

Exact confidence intervals for n using the MLE can be found in 
theory, although greater or less practical difficulties may ensue depending 
upon the size of the sample. Expression (8) is the exact distribution of 
k given the sample size and n. The upper confidence limit for n is then 
the value of this parameter which when substituted into expression 
(8) makes 


Pr {k < observed k} = : 


where 1 — a is the confidence coefficient, let us say .95. 

Unfortunately for N and n of even moderately large size, the sum- 
mation of (8) becomes prohibitive both in terms of time and accuracy 
unless special computing machines are available. 

It is possible, then, to use the approximate expression (9). This is 
satisfactory for the lower confidence limit, but unless N is quite large 
with respect to n, the upper confidence limit cannot be estimated by 
the use of this expression. This is because, as we have pointed out, the 
approximation only holds if \ is bounded. Now for a fixed value of N, 
X increases without bound as n grows very large. The result of this 
increase is that the upper confidence limit for n will be grossly over- 
estimated and may often be infinite. If exact confidence limits are 
desired, recourse must be had to the exact distribution (8). 

For a large sample, of course, the MLE is approximately normally 
distributed with a variance given by (24), so that an approximate 95% 
confidence interval will encompass 1.96 standard deviations on either 
side of the estimate 7. 


A SPECIFIC EXAMPLE 


With the above reservations in mind we may contrast the frequency 
of allelism estimate with the MLE for two sets of data. The pertinent 
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information is shown in Table 1. The data on the third chromosome 
of Drosophila pseudoobscura are taken from Wright & Dobzhansky (1941) 
while those on the second chromosome of D. melanogaster are from 
Wallace (1950). In both cases the frequency of allelism estimate is larger 


TABLE 1 
Comparison of Maximum Likelihood Estimate and Frequency of Allelism Estimate 
for Number of Loci on a Chromosome. For Symbols, see Text. 


N' | k' | Q A S.E. S.E. 
Wright | MLE | MLE | Wright 


3rd chromosome 
D. pseudoobscura 105 | 86 |.1529}.00407| 289 268 55.2 66.4 


2nd chromosome 
D. melanogaster 100 | 87 |.1209).0028 | 406 342 85.2 116.1 


than the Maximum Likelihood estimate as might be expected from the 
bias of the former. However, since nothing has been demonstrated about 
the bias of the MLE, too much weight should not be given this point. 
Likelihood estimates are not, in general, unbiased. The standard errors 
shown for Wright’s estimate are minimum values on the assumption of 
a random cross model. Actually the complete cross design was used in 
both so that the true standard errors are somewhat higher. 


THE CHOICE OF AN ESTIMATE 


Each of the estimation procedures described has certain disadvan- 
tages. If a random cross procedure is used, the MLE is not applicable 
because the requisite information, k, cannot be obtained from the data. 
The random cross procedure has the advantage that exact confidence 
limits can be placed upon the estimate of n. If the complete cross 
design is used, the MLE is superior since it is more efficient, as well 
as being sufficient. Moreover the frequency of allelism estimate is 
demonstrably biased and insufficient, while it shares with the MLE the 
difficulty of establishing exact confidence limits. The choice then lies 
between the random cross design using Wright’s estimate and the 
complete cross design using the MLE. When the sample size is small 
compared with n, the efficiencies of the two methods do not differ greatly 
as shown in Fig. 1. This does not take into account however, that for 
an equal number of crosses to be made in the two methods, the original 
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task of sampling the parent population is greater under the random 
cross design. That is, N(N — 1) lethal chromosomes are required to 
make N(N — 1)/2 random crosses, while only N lethal chromosomes 
are necessary for this many crosses in a complete design. The labor is 
then approximately N times as great in the first stage of the operation. 

Pavan & Knapp found 1,063 as the estimated n in the second 
chromosome of D. willistoni using a modified complete cross design. 
These authors contrast this with Wallace’s estimate of 406 for the second 
chromosome of D. melanogaster (see Table 1). As they poimt out, 
these two chromosomes have been found by other means to be homol- 
ogous, that is they contain more or less the same loci. However, 
since the confidence limits for these two estimates are much larger 
than these authors show, the results of Pavan & Raagpend <a 
are not really in conflict. 

What these data 
mated number of loci may have no significance whatever because of the 
general inefficiency of estimation procedures. 

In general, neither of the estimates is very efficient when sample size 
is small as compared with n. Moreover unless N > n both methods 
run the risk of yielding an infinite value for % or %. This objection 
applies to the random cross design, no matter what the size of N, 
although decreasingly so for large N. If estimates of the number of 
loci are to be made then, it would be best to perform a large scale experi- 
ment with a sample size in excess of the current approximate values of n. 
Such an experiment, using the MLE would provide a reasonably secure 
estimate of n. The estimates thus far obtained provide only orders of 
magnitude rather than precision. 
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QUERIES 


Grorce W. Snepecor, Editor 


QUERY: In Query No. 102, June 1953, you gave a method for 
121 supplying the yield of a missing plot in an unreplicated factorial 

experiment. You indicated that the test of treatments is biased 
upwards but that the bias is assumed to be small. I have a similar 
experiment in which one of the treatment means is significant at about 
the 3% level. I would like to know if the bias is sufficient to change the 
conclusion. How can I determine the amount of bias? 


The amount of bias can be found and subtracted from the 
ANSWER: approximate sum of squares for treatments. The deter- 

mination of the amount of bias is made in the following 
way. 

Insert the symbol z for the missing value and obtain the expressions 
for the sums of squares corresponding to the individual degrees of 
freedom supposing that z is a number. Each sum of squares is then 
either a number independent of x or a quadratic in x with known 
coefficients. Denote the error sum of squares by E(x) so that 


= ax’? + 2br +c 
and the sum of squares corresponding to the effect of interest by T(x) 
with 

T(x) = dx? + 2ex + f. 


Then the minimum of E(x) occurs when z = 2 = —b/aand is Eyy;, = 
c — b’/a. This is the value of the error sum of squares used in both 
the exact and approximate tests of significance. ‘The value of the effect 
sum of squares used in the approximate test of significance is 


a a 
On the other hand the minimum of E(x) + T(x) occurs for 
a+d 
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and is 
= b 


The effect sum of squares for the exact test of significance is then 
given by general linear hypothesis theory as 


(E + T) mia Exuin = f 
Hence the bias in the effect approximate sum of squares is 

Bias = T, — ((E + T) min — Enial 

a a a a+d 

_ _ | 

a+d a+d 


As an example consider the evaluation of the bias in the SK effect 
in query 102. 
The sum of squares for the effect is 


T(x) = py [x’ — 56.62 + 800.89]. | 


Hence 


i 
12 


56.6 


d= ~ 412% 2 


= 0.08, e= = —2.358. 


Also 


E(z) = 0.502? — 23.622 + 381.85. 


Hence x) = 23.62, a = 0.50 and T, = 1.82. 
The bias therefore is 


+ dre)? _ [—2.358 + yy (23.62)]? 
a+d 0.50 + 0.08 


The value of the SK sum of squares used in the exact test of signifi- 
cance is therefore 


= 0.26. 


T, = bias = 1.82 — 0.26 = 1.56. 
Gzorce ZyrsKIND 
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Papers presented at the Third Colloquium of the German Region of the Biometric Society 
at Bad Nauheim, January 27-29, 1956 


366 R. K. BAUER. (Munchen). General Theory of an Anthropo- 
metric Test of Paternity. 


A test of paternity is carried out in two steps: first the survey, then 
the evaluation of evidence on the observational unit. The observational 
unit is 3-valued usually and consists of child, its certain mother and 
uncertain father. The evidence consists of combinations of, say, M 
hereditary characters, which can be found on the observational unit. 
As a result of the survey we get a matrix X of 3 X M variates. The 
best criterion of evaluation is the ratio of the probability of X repre- 
senting a “‘true” family, and the probability of X representing a “‘false’’ 
family. All anthropometric test methods for proof of paternity are 
based on this criterion. Practical statements, however, lead perforce 
to LUDWIG’s proposal, i.e. to abbreviated discriminant functions. 
The results of a recently finished model experiment are shown. 


W. U. BEHRENS. (Hanover). Problems in Correlation Sta- 
367 tistics. 

Frequently too much confidence is placed in correlation analysis. 
Significant correlation is no proof of a direct causal relationship. Re- 
gression lines are not suitable for representing structural relationships. 
Particular difficulties arise with the interpretation of partial correlation 
coefficients. The usual treatment in textbooks leads to the mistaken 
conclusion that partial correlation is a key for detecting causal rela- 
tionships between single variates. The danger in such interpretation 
is outlined and demonstrated by means of models and agricultural 
experiments. There is no argument against the use of regression for 
fitting empirical data. 


368 W. U. BEHRENS. (Hanover). The Fitness of Different Field- 
Trial Designs for the Removal of Soil Differences. 


For the layout of field trials with 4 to 10 treatments Latin square 
and block designs are frequently used. Under the assumption that 
small differences between treatments allow for efficient screening of 
soil differences, various field trial designs are discussed and particularly 
suitable (‘gerechte’) designs demonstrated. Highest precision of these 
designs is obtained with square or near square plots in contrast to long 
rectangular plots, though generally the latter allow a better screening 
of soil differences. 
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369 H. GEIDEL. (Rethmar). On Symbols and Terminology for 
Use in Agricultural Biometry. (Report.) 


370 A. JANOSCHEK. (Giessen). Quantum Biology and Reaction 
Kinetics. 


The time-trend of the decline of a homogeneous population of 
_ organisms under the effect of a noxious principle agrees with a reaction 
which can be formularized by the GLOCKNER proposition. The 
dose effect relation on the other hand is more informative generally. 
The effect of time (t) and dose (D), combined in a single formula, yields 


N = N{1 — exp [—(k-D)”-t]}” 


This formulation of the law of mass action applies also to the rate of 
population decline and growth. Examples show general validity of 
this basic law of reaction kinetics. 


371 G. A. LIENERT. (Marburg). Quantitative Analysis of Test 
Methods in Clinical Medicine. 


As yet no statistics have been defined to measure the reliability 
of clinical testing methods, other than the usual description by a 2 X 2 
table of percentages (negative/positive reaction X healthy/ill persons). 
The use of a suitable correlation coefficient as a means of characterizing 
reliability and symptomatic value is suggested; a possible conventional 
interpretation of the coefficient is outlined. 


O. LUDWIG. (Bad Nauheim). Theory and Application of 
372 Runs. 


Several kinds of runs are defined, and their importance for tests 
of random order in sampling for attributes is explained by means 
of examples from genetics, plant breeding, and meteorology. Several 
conditional and unconditional tests are considered for two or more 
characters. The theory may be applied to observations on variables 
too, e.g. by considering runs above and below the median. Runs up 
and down and runs of consecutive elements are also defined for variables 
(measurements); but their distribution under the null hypothesis is 
entirely different from that of the above-mentioned runs. In a certain 
sense an inversion of common theory is the “Pascal problem’ (inverse 
sampling) for runs; the connection for this problem with that of direct 
sampling and with FELLER’s theory of recurrent events is dealt with. 
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373 K. H. MULLER. (Jena). On the Accounting for Soil Varia- 
bility in Field Trial Analysis of Variance. 


In field trials the designed layout, i.e. the number of blocks, may or 
may not correspond with the number of actual different soil qualities. 
In the first case the block degrees of freedom have to be used in the 
statistical analysis; in the latter case it seems appropriate to allocate 
the experimental material and to assign the degrees of freedom cor- 
respondingly to realized soil groups. In that case, this seems to be a 


logical procedure to ensure a separation of systematic and random 
fluctuations. . 


374 H. RUNDFELDT. (Hanover). Review of Methods Usually 
Applied in Field Plot Technique. 


In field plot experiments the existence of soil variation by which 
means and variances might be considerably biased must be always 
taken into consideration. It is therefore necessary to reduce its influ- 
ence by choice of a suitable layout and statistical analysis of the experi- 
ment. In respect of the layout long narrow plots are preferable. Re- 
garding the analysis various systems were compared by means of dummy 
experiments. Based on the theoretical results of the evaluation of an 
experiment without soil differences (s? = 100) and the information 
resulting from a given field one finds the following “rule-of-thumb 
numbers” for the tested systems: A. for few treatments, 1. without soil 
balance (completely randomized designs) 250, 2. randomized blocks 
125, 3. method MITSCHERLICH 125, 4. Latin squares 120, 5. method 
LINDHARD 120, 6. comparison by means of systematic controls 200. 
B. for a greater number of treatments, 1. systematic controls 200, 
2. random controls 200, 3. lattices 130-150. Besides these ‘‘rule-of- 
thumb numbers” one will notice various technical advantages and 
disadvantages in applying the different systems. 


375 B. SCHNEIDER. (Giessen). On the Usefulness of Percentiles 
in Biometry.- 


The estimation of standard deviation from large samples by approxi- 
mation formulas frequently introduces uncontrollable and principal 
errors into further considerations. As a measure of dispersion the 
percentiles or a function of the same (e.g. semi-interquartile-range) are 
to be preferred; the usefulness of this standard deviation estimator is 
demonstrated. 
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376 W. SCHNELL. (Scharnhorst). On the Sphere of Permissible 
Generalization of Field Trial Results. 


Field trials are considered in which several “yield” factors such as 
“treatments” and “blocks” have been modified, the modifications 
(or “levels’”’) of each factor being orthogonal to the modifications of the 
other factors. (EISENHART’s model I and II, mixed model.) The 

sphere of permissible generalization of the result of a comparison is 
' said to comprise, with respect to some orthogonal factor modified at 
random, the corresponding population of modifications of that factor; 
however, as regards some orthogonal factor with fixed modifications, 
the generalization sphere is confined to the modifications actually used. 
As a special case, if the factor in question is acting additively within 
the limits so described, the generalization sphere widens as far as 
additivity holds. Asa rule, there are sizable interactions of “treatments” 
with “soils” and ‘‘years’’, respectively. One-year trials, though giving 
an unbiased estimate of the long-time average of some interesting 
comparison, do not furnish an estimate of the appertaining interaction 
with years and hence will not yield a valid error, if recommendations 
for future years are to be made. The resulting need for trials replicated 
over several years suggests the applicability of sequential sampling 
procedures to field experimentation. 


377 K. v. SOLTH. (Marburg). On Studies of the Gestation Process 
Under Different Gynaecological Diseases. 


The question of whether a relation exists between the ratio of 
abortions to births and different genital diseases is studied in cases of 
myoma (935), collum (450) and corpus (143) cancer. A method for 
characterizing this existing relation is given. 


378 R. WARTMANN. (Dusseldorf). Analysis of Variance When 
the Population is Finite. 


Multistage sampling of a finite population with p’, q’, r’, --- , units 
at each stage from which p < p’,q < q’,7r < 1’, --> , are sampled. 
The precision of a single one of the p-q-r- --- samples is asked for. 


Formulas for the variance, and separation of its components of a single 
sample are given. 


a7q ™. WERMKE. (Bochum). Combined Analysis of Variance of 
Field Trials by Computing Machine Equipment. 
The technique of computational procedure with excessive field 
trial material by use of IBM punched card equipment is detailed 
for a two-factorial analysis of variance. Generalizations. 
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380 R. WETTE. (Heidelberg). On the Use of Regression Lines in 
Biology. 


Frequently the conditions for identifiability of a linear or allometric 
structural relationship between two variables subject to error are not 
met in biological material. Application of the usual regression procedures 
to material of this kind yields fallacious results when the causal relation- 
ship between the structural variables is asked for, a fact frequently 
neglected. 


Papers presented at the first meeting of the Brazilian Region of the Biometric Society, 
at Instituto Biologico, Sao Paulo, Brazil, January 3, 1956 


381 A.CONAGIN. New Tests for Comparison of Means. 


In this paper the author discusses new tests for comparison of two 
means or contrasts ‘‘a posteriori’? when the null hypothesis of a group 
of n means is rejected at an a level. 

The tests considered are the least significant difference test and also 
tests of Newman, Keuls, Tukey, Scheffe and Duncan. The essential 
differences between them are pointed out in detail and the rules applied 
to a particular situation. 


382 A. GROSZMANN AND J. DOBEREINER. Problems in the 
Statistical Analysis of Population Growth in Bacteria. 


Two strains of bacteria of the genus Beijerinckia were studied for 
different purposes. Strain ‘“E’”’ was very efficient in free nitrogen 
fixation, fixing 15 mg of N, strain “‘F’’ was low, fixing only 8.5 mg of 
N per gram of sucrose. The experiment on counting the number of 
bacteria in the population was run in a split-plot design, with two 
replications, “E’’, “F’’ and “O” (check) strains, countings made at 
1, 2 and 3 week intervals, in two different culture media. The data 
were analysed on number of bacteria found under the field of the 
microscope. Ten fields were counted from every treatment. The 
analysis showed a C. V. = 38.9% for countings and a C. V. = 10.3% 
for “plots”. In spite of the relatively large error term, highly significant 
differences were found between strains and among weeks. The difficulties 
in sampling the population and making bacteria counts under the 
microscope were discussed. 
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W. R. JARDIM, A. M. PEIXOTO, S. SILVEIRA FILHO AND 

383 F. PIMENTAL GOMES. Study on the Precision and Accuracy 
of a Few Practical Methods of Estimating Milk Production of 
Dairy Cows. 


This paper deals with the estimation of milk production by means 
of biweekly, monthly and bimonthly observations, without taking 
into account, as is usual practice, the date of calving. The data studied 
were 72 records of cows of the Holstein-Friesian breed: 6 calvings 
in each month of the year and also 12 first calvings, 12 second calvings 
and so on, up to the sixth. These cows belong to the herd of the Escola 
Superior de Agricultura “luiz de Queiroz” (Piracicaba, S. P., Brazil), 
which has been kept since 1914 under approximately the same condi- — 
tions. 

The authors criticize the use of “maximum error” and also the use 
of mean deviation, both to be found in papers dealing with this subject. 
The former is completely superseded and inadvisable, and the latter, 
although equivalent to a certain extent to the usual standard deviation, 
has only 87.6% of its efficiency. 

The data obtained were compared with the actual production, 
corresponding to daily control, and the deviations observed were 
studied. Their means and standard deviations were the following, in 
kilograms per lactation period. 


Method Mean of Deviations Standard Error of Mean 
Biweekly control + 7.59 7.9 
Monthly control + 8.92 13.5 
Bimonthly control +121.86 21.7 


Comparison of these means with the expected value (zero) under 
the null hypothesis, by the ¢ test, shows that biweekly and monthly 
controls may be taken as unbiased, while bimonthly control is biased, 
the bias being positive and around 5% of the actual production. 

An analysis of variance of the observed deviations was carried out, 
this being correct in view of recent research by G. E. P. Box (Ann 
Math. Stat. 25: 290-302, 1954). The analysis, completed by Tukey’s 
studentized range test, shows that, with respect to a possible bias, the 
biweekly and monthly controls may be accepted to be equal to each 


other, but that they are both significantly different from the bimonth! 
control. 
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RUY AGUIAR DA SILVA LEME. (Escola Politecnica da 
384 Universidade de Sao Paulo). The ‘‘Simplex” Method in Multiple 
Regression. 


A modification in the method of Pivotal Condensation (Rao, Ad- 
vanced Statistical Methods, 1952) allows a prediction to be made of the 
improvement afforded by a new independent variate on the sum of 
squares due to regression. The method of Pivotal Condensation thus 
modified becomes formally equivalent to the “Simplex’’ Method of 
Linear Programming. 

An actual example is entirely worked out for illustration. 


A. M. PENHA, G. SCHREIBER, A. R. HOGE AND H. E. 
385 BELLUOMINI. Application of the Discriminant Function to a 
Problem of Sex Differentiation in Snakes. 


In a species of snake (Bothrops insularis) inhabiting exclusively a 
small island off the coast of the State of Sao Paulo, a large proportion 
of individuals, tentatively considered as inter-sexual, externally resemble 
males (male copulatory organs present) but after dissection show fertile 
female internal organs. 

With the values of 3 external morphological characters—head, 
body and tail lengths—Fisher’s discriminant function was calculated 
for a group of 155 typical males and 65 typical females of the afore- 
mentioned species. As indicated by the analysis of variance, the 
function had a highly significant discriminant power (F = 78.7). Its 
probability of misclassifying a randomly taken individual was 0.13, 
according to the corresponding “‘t’’ value, whereas the observed variates 
were much less discriminating, showing the following probabilities of 
misclassification under the same conditions: head length, 0.35; body 
length, 0.44; tail length, 0.38. Furthermore, the isolated contribution of 
each measurement to the discriminant was highly significant. 

When the calculated function was applied to a group of 151 inter- 
sexual individuals, 93% were classified as females, which indicated a 
proportion of misclassification well below the theoretical expectation 
for typical females. 


M. ROCHA E SILVA AND A. M. ROTHSCHILD. (Instituto 

86 Biologico, Sao Paulo, Brazil). A 4-point Design for Bioassay of a 

3 Material Inducing Strong Tachyphylaxis (Anaphylatoxin), With 
a Reference to the Mechanism of Desensitization. 


It has been shown in this laboratory that anaphylatoxin prepared by 
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incubating normal rat plazma with agar owes its stimulating action 
upon the guinea pig ileum to a release of histamine. The contracting 
effect is, therefore, indirect, and the response becomes less and less as 
the additions of the same dose of anaphylatoxin are repeated (tachyphy- 
laxis or desensitization). The possibility of comparing two different 
preparations of anaphylatoxin, a “standard” and the “unknown” by 
applying them alternatively to the same piece of intestine would give 
highly inaccurate results on account of the enormous bias introduced 
by the previous additions of the agent. By using the known set-up of 
. a 4-point assay with randomized blocks as utilized for histamine, 
bradykinin, oxytocin and so forth, the bias introduced by the tachyphy- 
lactic effect is so strong that only exceptionally could any reasonable 
ratio of potency be derived from the experimental data. 

After trying several designs in order to test the possibility of setting 
a 4-point assay for a quantitative estimation of anaphylatoxin, we have 
found that a 4 X 4 latin square design is able to correct for the bias 
introduced by the tachyphylactic effect, allowing a reasonable estimate 
of the ratio of potency between “unknown” and “standard”. The 
variation due to the tachyphylactic effect could be confounded with the 
so-called ‘order of additions” (““Rows’’) and its sum of squares eliminated 
from the experimental error, thus reducing to about 10 to 12% the 
standard error of each response. In a series of 19 squares the variances 
were found homogeneous, and the value of \ = S/b around 0.10. 
Deviations from parallelism and linearity were X trivial and non- 
significant in nearly all experiments. The standard error of M (= sy) 
was, on the average, around 0.030. 

A coefficient of desensitization (A) is described and the method for 
its calculation set up. The negative regression of ‘Rows’ on “Bias” 
(i.e. the number of previous additions of the same dose of anaphylatoxin) 
indicating a linear downward trend of the responses to the same repeated 
dose of anaphylatoxin was explained as a linear regression of response 
on dose of “intrinsic histamine’. From these considerations, a function 
D = Ce’ for doses of “intrinsic histamine’”’ has been derived and the 
responses 


Y= k-log, D = k-log, C — ka, where —log, C= A 


should be identified with the linear equation of the regression line 
“Rows” on “Bias’’, indicative of the responses to the released intrinsic 
histamine at the different levels of bias (0, 1, 2, 3, ---). The coincidence 
between the theoretical equations and the experimental ones calculated 
by the least squares method is remarkable. 
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Paper presented at the Joint Meetings of the American Institute of Mathematical 
Statistics, American Statistical Association, and the Biometric Society 
(ENAR), New York, N.Y., U.S.A., December 27-29. 


387 GEORGE FERRIS. Three Useful Designs in Taste-Testing. 


There has been a tendency to use statistical designs intended for 
other fields of application in taste-testing, whether the model is appli- 
cable or not. Alternatively, non-parametric methods have been used. 
It is the intention of the paper to show how existing models can be 
suitably modified for organoleptic purposes. 

The first model is intended for use by analytic taste-panels for 
quality control purposes when a number of samples are judged serially 
for flavor. It takes existing Latin Squares, Youden Designs, and In- 
complete Blocks, and modifies them in order to control “resid 
or “‘carry-over’’ effects. 

The second model is intended for use by analytic panels judging 
color or another physical property in samples simultaneously set out 
alongsid® one another. It modifies existing designs and analyses to 
allow for the phenomenon of ‘‘adaptation”’ known to psychologists. 

The third model illustrates how a balanced design for incomplete blocks 
of two may by a slight variation be adopted for organoleptic purposes. 

Each model is illustrated with a numerical example. 

Consumer testing is briefly referred to. 


Paper presented at the meeting of the French Region of the Biometric Society, February 8. 


388 D. SCHWARTZ, J. ULMO ET A. VESSEREAU. Problémes 
Relatifs 4 |’Echantillonnage Stratifié. 


Il s’agira, au cours de cette communication, de problémes d’échantil- 
lonnage stratifié ot le nombre de prélévements est le méme dans les 
diverses unités 4 un méme niveau, par example: 
dans la production d’un centre de fabrication d’ampoules de B.C.G., on 
choisit au hasard / lots de fabrication, on préléve dans chacun d’eux 
un méme nombre a d’ampoules, donnant lieu chacune a ¢ tubes de cul- 
tures sur lesquelles on mesure le nombre de particules vivantes 
ou bien: 
dans une livraison de papier on choisit au hasard r rames, dans chacune 
on préléve au hasard f feuilles, chacune donnant naissance a un méme 
nombre de mesures de résistance a |’éclatement. 

Et en particulier, la catégorie des problémes 4 deux niveaux, ot on 
représente une population par e échantillons donnant lieu chacun & m 


. Mesures. 


On exposera en particulier quels sont les plans d’échantillonnage 
optimum selon qu’on s’intéresse plus particuliérement 4 la valeur moyenne 
de la caractéristique étudiée, ou 4 sa variabilité aux divers niveaux. 
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THE BIOMETRIC SOCIETY 


General 


The paid up membership of the Society at 3lst December, 1955 
stood at 1265, divided between regions as follows:— 


E. N. American 483 German 61 
W. N. American 95 Indian 12 
Australasian 50 Italian 57 
Belgian 62 Japanese 43 
Brazilian 52 Netherlands 29 
British 155 Swedish 16 
Danish 12 Swiss 22 
French 65 At large 51 


The office of the Secretary has been transferred from New Haven 
to Rothamsted during the summer. Members are asked to excuse any 
delays that have arisen during the transfer. . 


Brazilian Region 


At a meeting held on January 3rd at the Instituto Biologico, Sao 
Paulo, the following were elected for 1956: 


President, Dr. C. G. Fraga, Jr; 
Secretary, Dr. P. M. Freire; 
Treasurer, Prof. A. Groszmann; 
Committee, Dr. Breiger 

Dr. Bueno 

Prof. Bitancourt 

Prof. Memoria 

Dr. Penha 


The following papers were presented at the same meeting—A. M. 
Penha and collaborators, on a problem of sex differentiation in snakes; 
M. Roche e Silva, on a factorial assay for anaphylotoxin; R. Silva Leme, 
on a new procedure for multiple regression; A. Groszmann, on bacterial 
growth curves; A. Conagin, on the comparison of several means; and 


F. Pimental Gomes and collaborators, on the evaluation of milk pro- 
duction. 


Belgian Region 


En novembre 1955, Prof. Lousse, Recteur de |’Ecole de Médicine 
Vétérinaire de Cureghiem a donné une conférence entitulée “Les re- 
cherches pures et appliquées en physiologie animale.” 

La Société Adolphe Quetelet a tenu son assemblée générale & Brux- 
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elles le 22 mars. Le Conseil d’administration pour 1956 se compose 
comme suite: Président, Prof. D. de Meulemeester; 
Vice-Présidents, MM. de Naeyer, Laurent, Lebrun, Welsch, 
Lecrenier; 
Secrétaire, M. L. Martin; 
Secrétaire adjointe, Mlle. A. Lenger; 
Trésorier, M. A. Rotti; 
Membres, MM. Roussel, Reuse, Bontemps. 
Ensuite, Mile. Lenger a fait un exposé sur le sujet ‘“Quelques 
commentaires 4 propos d’un voyage d’étude aux Etats-Unis.” 


German Region 


A colloquium held at Bad Nauheim during January 27-29 was 
attended by 150 biometricians, including 35 from East Germany and 
guests from Austria and Chile. Topics included:— 


I Distribution-free methods. L. Schmetterer, A general survey; 
H. Miinzner, Permutation tests. 


II Biometry and Medicine. W. Oemisch, Graduation of growth 
data by the Normal curve; H. J. Heite, Statistical treatment of experi- 
mental mortality curves; J. Hartung, Sickness Insurance; G. Bertram 
and J. Hartung, Remarks on sequential designs; H. Hosemann, Para- 
bolic curves of neo-natal mortality; G. A. Lienert, Quantitative analysis 
in clinical research; K. V. Solth, Effects on gestation of various gyne- 
cological diseases. 


III Statistical Methods in Agriculture. H. Geidel, A report on 
symbols and notation for agricultural biometry; W. Behrens, Correlation 
problems; H. Rundfeldt, A Critique of current techniques in agricultural 
research; W. Schnell, The accuracy and validity of field experiments; 
W. Behrens, The suitability of various experiment designs in the study 
of soil variation. 


Région Francaise 


La Région a tenu son assemblée générale 4 Paris le 8 février. Ont 
été elus—Président, M. E. Morice; Membre du Conseil, M. D. Bargeton. 
M. D. Schwartz, Mile. J. Ulmo et M. A. Vessereau ont presenté une 
communication entitulée “Problémes relatifs a l’échantillonage stratifié.”’ 


Australasian Region 


Visit of Professor M. G. Kendall: Professor Kendall gave a series of 
lectures and seminars in Sydney, Canberra, Adelaide and Melbourne to 


4 
: ‘ 
4 
a 


THE BIOMETRIC SOCIETY 237 


statisticians, economists, public servants and the general public on a 
wide variety of topics. The Region owes a great deal to C.S.I.R.O. for 
allowing non-C.S.I.R.O. members to see so much of Professor Kendall 
during his tour. A conference in Melbourne, 9th-13th April, was 
attended by 50 delegates from the Division of Mathematical Statistics 
C.S.1.R.0., the Universities, other Divisions of C.S.I.R.O. and Govern- 
ment departments, and from Commercial Organizations. It was easily 
the largest gathering ever held of Australian mathematical statisticians. 
On the two free afternoons in the week, visits were arranged to the 
Olympic Site, the Electronic Computer at the University of Melbourne 
and the Division of Meteorological Physics, C.S.I.R.O. 

The first three days of this Conference, with Dr. G. S. Watson in 
the Chair, were devoted to Time Series. Professor Kendall gave two 
lectures in which he sketched an outline of the subject as it stands 
today. In the first lecture the distribution theory of serial correlation 
coefficients was discussed. The second lecture dealt with questions 
of inference. In general the discussion was limited to the case of auto- 
regressive processes in discrete time. Contributed papers were read 
by: Dr. D. G. Lampard, “A method of estimation of correlation func- 
tions’; Mr. E. K. Webb, “Spectrum Analysis of continuous time 
series”; Mr. R. T. Leslie and Dr. F. E. Binet, ‘‘Relaxed runs’; Dr. G. S. 
Watson, “The joint distribution of circular serial correlation coefficients” 
and Dr. E. J. Hannan, ‘Testing for serial correlation in least squares 
regression’. 

The subject of the last two days of the Conference, with Professor 
E. J. G. Pitman in the Chair, was Distribution-Free Inference. In the 
first of his two lectures on this topic, Professor Kendall gave a general 
discussion and attempted to define and classify some of the concepts 
involved. In the second, he classified many distribution-free tests 
and compared the relative asymptotic efficiencies of alternative tests 
of the same hypothesis. Contributed papers were read by: Dr. H. S. 
Konijn, “Distribution-free procedures for testing treatment differences”’; 
Dr. H. O. Lancaster, “Some reconciliations of x*”; Mr. G. A. MacIntyre, 
“Distribution-free comparison of two sets of data” and Dr. G. S. Watson, 
“Distribution-free tests, similar regions and sufficient statistics’’. 
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NEWS AND ANNOUNCEMENTS 


The degree of Doctor of Laws was conferred on Harold Hotelling 
by the University of Chicago on 11 November, 1955, at the convocation 
in celebration of the twenty-fifth anniversary of the University’s Social 
Science Research Building. Dr. Hotelling was cited as the “foremost 
contemporary contributor of quantitative methods to the social sciences, 
who by mathematical analysis has notably advanced our understanding 
of fundamental problems in economics and statistics’. 

Dr. J. H. Bennett, who has been with Sir Ronald Fisher at the 
University of Cambridge, has gone to South Australia, to become 
head of the Department of Genetics, Adelaide University. 


Summer Session at the Massachusetts Institute of Technology 


The Statistical Summer Session formerly held at the University of 
Connecticut will this summer be held at Endicott House of the Massa- 
chusetts Institute of Technology during the weeks August 13 to 25 
inclusive. Two one-week programs will be offered. The first week 
will be under the chairmanship of Professor Leo Tick of New York 
University on the general topic ‘“Time Series’; the second week under 
the chairmanship of Professor Max Woodbury of George Washington 
University with the topic “The Impact of Computers on Statistics’’. 
Requests for information and reservations should be sent to Dr. M. E. 
Terry, Bell Telephone Laboratories, Murray Hill, N. J., U.S.A. 


Annual Meeting of the Biometric Society (WNAR) 


The annual meeting of the Bometric Society (WNAR) will be held 
at the University of Washington in Seattle during the period of August 
22-24, in conjunction with national meetings of the Institute of Mathe- 
matical Statistics, the American Mathematical Society, the Mathe- 
matical Association and the Econometric Society. The Special Invited 
Address will be given by Professor N. Rashevsky, Chairman of the 
committee of Mathematical Biology at the University of Chicago. 
Another address will be given jointly to the Institute and the Society 
by Professor Harold Hotelling, University of North Carolina, on “Some 
new light on the multiple correlation coefficient’’. There will be sessions 
on statistical problems in forestry, in medicine, and on prediction 
problems and possibly others. There will also be a session for contri- 
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buted papers. Abstracts of contributed papers and requests for accom- 
modations should be sent to the Program Chairman, Douglas G. Chap- 
man, Department of Mathematics, University of Washington, Seattle. 
A feature of the entertainment will be a salmon barbecue on Puget Sound. 


Meeting of the Association for Computing Machinery 


The annual meeting of the Association for Computing Machinery 
will be held on the University of California Westwood Campus, Los 
Angeles, August 27-29, 1956. (See January issue of the Journal of the 
Association for Computing Machinery). For information write G. W. 
King, Box 3251, Olympic Station, Beverly Hills, California, U.S.A. . 


Australian National University Research Scholarships 


Applications are invited from graduates for enrolment as research 
students in various Schools, including the Research Schools of Social 
Sciences and Pacific Studies, in which research may be done in Statistics, 
including Mathematical Statistics. Particulars and application forms 
may be obtained from: (1) The Australian Embassy, 2941 Massachusetts 
Avenue, Washington, D. C., U.S.A.; (2) The Australian Consulate- 
General, 206 Sansome Street, San Francisco 4, California, U.S.A.; 
(3) The Australian Consulate-General, Room 426, International Build- 
ing, 636 Fifth Avenue, New York 20, N. Y., U.S.A.; (4) The Office of the 
Australian High Commissioner, 5th Floor, Royal Bank Chambers, 
100 Sparks Street, Ottawa, Canada. 


Australasian Region of the Biometric Society 


The Society will hold a general meeting at 8 p.m. 16th August, in 
the Arts Building, University of Melbourne, at which Dr. E. A. Cornish 
will give his Presidential Address. It will form a part of the Australian 
Mathematical Society Meeting. 

Internationales Biometrisches Seminar, 24. September — 3 Oktober 1956 
u. Internationales Biometrisches Symposium, 1. Oktober — 3. Oktober 1956 
fiir die deutschsprachigen Gebiete. 

Veranstaltungsort: Linz a. d. Donau (Osterreich). Organisation: 
Prof. Dr. Arthur Linder, Genf (Schweiz), Avenue de Champel 24. 
Geschiftsstelle: Dr. Adolf Adam, Osterreichische Stickstoffwerke 
Aktiengesellschaft, Abt. UAW, Linz a. d. Donau, St. Peter 224 (Oster- 
reich)—Telefon Nr. 29181/Klappe 1508. 


Die Internationale Biometrische Gesellschaft, deren Zweck die 


Entwicklung, Anwendung und Verbreitung der quantitativen Methoden 
in der Biologie ist, veranstaltet in Linz a. d. Donau (Osterreich) ein 
Seminar tiber Planungsforschung fiir Biologen, Biochemiker, Agrikul- 
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turchemiker, Mediziner, Pharmakologen, Pflanzenschutztechniker, Land- 
und Forstwirte. Das Seminar gliedert sich in eine Grundausbildung 
im Planen und Auswerten von Versuchen und Beobachtungen durch 
Vorlesungen, Ubungen, Aussprachen und Besichtigungen, sowie in die 
eingehende Behandlung aktueller praktischer Probleme. Die Ausbil 
dung wird durch internationale Fachkrafte in deutscher Sprache 
geboten. 

Den Seminarbesuchern wird die Gelegenheit gegeben, an den 
Aussprachen im Rahmen des [nternationalen Biometrischen Symposiums 
iiber Wachstums- und Ertragsgesetze (1. Oktober 1956), Transforma- 
tionen bei der Auswertung von Haufigkeiten und Wirkungskurven 
(2. Oktober 1956), Beurteilung der Wirkung von Heilmitteln bei chro- 
nischen Erkrankungen (3. Oktober 1956), als Hérer teilzunehmen. 
Die teilnehmerzahl fiir das Internationale Biometrische Seminar wird 
auf 100 beschrankt. 

Die Gebiihr fiir die Teilnehmerkarte am Seminar wird S 150-, 
sfr. 25-, DM 25- betragen. Interessenten werden gebeten, ihre 
Anschrift zwecks Zusendung des Seminar programmes, von Anmelde- 
formularen usw. an die Geschiaftsstelle: Dr. Adolf Adam, Osterreichische 
Stickstoffwerke Aktiengesellschaft, Abt. UAW, Linz a. d. Donau, 
St. Peter 224 (Osterreich) bekanntzugeben. 
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